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Summary 



This memoir contains an overview of the proof of the bounded curvature conjecture. 
More precisely we show that the time of existence of a classical solution to the Einstein- 
vacuum equations depends only on the L^-norm of the curvature and a lower bound 
of the volume radius of the corresponding initial data set. We note that though the 
result is not optimal with respect to the standard scaling of the Einstein equations, it 
is nevertheless critical with respect to another, more subtle, scaling tied to its causal 
geometry. Indeed, bounds on the curvature is the minimum requirement necessary 
to obtain lower bounds on the radius of injectivity of causal boundaries. We note also 
that, while the first nontrivial improvements for well posedness for quasilincar hyperbolic 
systems in spacetime dimensions greater than 1 + 1 (based on Strichartz estimates) were 
obtained in [2], [3], [49], [50], [19] and optimized in [20], [36], the result we present here 
is the first in which the full structure of the quasilinear hyperbolic system, not just its 
principal part, plays a crucial role. 

The entire proof of the conjecture is contained in the following sequence of papers 

S. Klainerman, I. Rodnianski, J. Szeftel, The bounded L'^ curvature conjecture. arXiv: 1204 -1767, 
91 pp. This is the main part of the series in which the proof is completed based on the 
results of the papers below. 

J. Szeftel, Parametrix for wave equations on a rough background I: regularity of the phase 
at initial time. arXiv: 1204. 1768, PP- 

J. Szeftel, Parametrix for wave equations on a rough background II: control of the parametrix 
at initial time. arXiv: 1204. 1769, ^4 PP- 

J. Szeftel, Parametrix for wave equations on a rough background III: space-time regularity 
of the phase. arXiv:1204.1770, 276 pp. 

J. Szeftel, Parametrix for wave equations on a rough background IV: Control of the error 
term. arXiv: 1204. 1771, 284 PP- 

J. Szeftel, Sharp Strichartz estimates for the wave equation on a rough background. arXiv:1301.0112, 
30 pp. 
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CHAPTER 1 



Introduction 

1.1. General Introduction 

We present a summary of our proof of the bounded L^-curvature conjecture in General 
Relativity. According to the conjecture the time of existence of a classical solution to 
the Einstein-vacuum equations depends only on the L^-norm of the curvature and a 
lower bound of the volume radius of the corresponding initial data set. At a deep level 
the curvature conjecture concerns the relationship between the curvature tensor and 
the causal geometry of an Einstein vacuum space-time. Thus, though the result is not 
optimal with respect to the standard scaling of the Einstein equations, it is nevertheless 
critical with respect to a different scaling, which we call null scaling, tied to its causal 
properties. More precisely, curvature bounds are strictly necessary to obtain lower 
bounds on the radius of injectivity of causal boundaries. These lower bounds turn out 
to be crucial for the construction of parametrices and derivation of bilinear and trilinear 
spacetime estimates for solutions to scalar wave equations. We note also that, while 
the first nontrivial improvements for well posedness for quasilinear hyperbolic systems in 
spacetime dimensions greater than 1 + 1 (based on Strichartz estimates) were obtained 
in [2], [3], [49], [50], [19] and optimized in [20], [36], the result we present here is the 
first in which the full structure of the quasilinear hyperbolic system, not just its principal 
part, plays a crucial role. 

1.1.1. Initial value problem. We consider the Einstein vacuum equations (EVE), 

Ric,;3 = (1.1) 

where Rica/j denotes the Ricci curvature tensor of a four dimensional Lorentzian space 
time {Ai, g). An initial data set for (1.1) consists of a three dimensional 3-surface Eq 
together with a Riemannian metric g and a symmetric 2-tensor k verifying the constraint 
equations, 

/ - V,tr/c = 0, . . 

\ Rscai-\k\' + {trky = 0, ^ ^ ^ 

where the covariant derivative V is defined with respect to the metric g, Rscai is the 
scalar curvature of g, and trk is the trace of k with respect to the metric g. In this work 
we restrict ourselves to asymptotically flat initial data sets with one end. For a given 
initial data set the Cauchy problem consists in finding a metric g satisfying (1.1) and an 
embedding of Eg in A4 such that the metric induced by g on Eq coincides with g and the 
2-tensor k is the second fundamental form of the hypersurface C M.. The first local 
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existence and uniqueness result for (EVE) was established by Y.C. Bruhat, see [5], with 
the help of wave coordinates which allowed her to cast the Einstein vacuum equations in 
the form of a system of nonlinear wave equations to which one can apply^ the standard 
theory of nonlinear hyperbolic systems. The optimal, classical^ result states the following, 

Theorem 1.1 (Classical local existence [12] [14]). Let {T.o,g,k) be an initial data 
set for the Einstein vacuum equations (1.1). Assume that Sq can be covered by a locally 
finite system of coordinate charts, related to each other by diffeomorphisms, such that 
{g, k) E Hi^^{T.q) X Hi^^(T,q) with s > |. Then there exists a unique^ (up to an isom- 
etry) globally hyperbolic development {Ai,g), verifying (1.1); for which Eq is a Cauchy 
hypersurface^ . 

1.1.2. Bounded curvature conjecture. The classical exponents s > 5/2 are 

clearly not optimal. By straightforward scaling considerations one might expect to make 
sense of the initial value problem for s > Sc = 3/2, with the natural scaling exponent for 

based Sobolev norms. Note that for s = Sc = 3/2 a local in time existence result, for 
sufficiently small data, would be equivalent to a global result. More precisely any smooth 
initial data, small in the corresponding critical norm, would be globally smooth. Such 
a well-posedness (WP) result would be thus comparable with the so called e- regularity 
results for nonlinear elliptic and parabolic problems, which play such a fundamental role in 
the global regularity properties of general solutions. For quasilinear hyperbolic problems 
critical WP results have only been established in the case of 1 + 1 dimensional systems, 
or spherically symmetric solutions of higher dimensional problems, in which case the L^- 
Sobolev norms can be replaced by bounded variation (BV) type norms^. A particularly 
important example of this type is the critical BV well-posedness result estabhshed by 
Christodoulou for spherically symmetric solutions of the Einstein equations coupled with a 
scalar field, see [7] . The result played a crucial role in his celebrated work on Weak Cosmic 
Censorship for the same model, see [8]. As well known, unfortunately, the BV-norms are 
completely inadequate in higher dimensions; the only norms which can propagate the 
regularity properties of the data are necessarily based. 

The quest for optimal well-posedness in higher dimensions has been one of the major 
themes in non-linear hyperbolic PDE's in the last twenty years. Major advances have 
been made in the particular case of semi-linear wave equations. In the case of geometric 

^The original proof in [5] relied on representation formulas, following an approach pioneered by 

Sobolev, SCO [37]. 

■^Based only on energy estimates and classical Sobolev inequalities. 

^The original proof in [12], [14] actually requires one more derivative for the uniqueness. The fact 
that uniqueness holds at the same level of regularity than the existence has been obtained in [33] 
^That is any past directed, in-extendable causal curve in M. intersects Eq. 

^Recall that the entire theory of shock waves for 1-1-1 systems of conservation laws is based on BV 

norms, which arc critical with respect to the scaling of the equations. Note also that these BV norms 
are not, typically, conserved and that Glimm's famous existence result [13] can be interpreted as a global 
well posedness result for initial data with small BV norms. 
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wave equations such as Wave Maps and Yang- Mills, which possess a well understood null 
structure, well-posedness holds true for all exponents larger than the corresponding crit- 
ical exponent. For example, in the case of Wave Maps defined from the Minkowski space 
M^^^ to a complete Riemannian manifold, the critical scaling exponents is Sc = n/2 and 
well-posedness is known to hold all the way down to Sc for all dimensions n > 2. This 
critical well-posedness result, for s — n/2, plays a fundamental role in the recent, large 
data, global results of [47], [40], [41] and [28] for 2 -|- 1 dimensional wave maps. 

The role played by critical exponents for quasi-linear equations is much less under- 
stood. The first well posedness results, on any (higher dimensional) quasilinear hyperbohc 
system, which go beyond the classical Sobolev exponents, obtained in [2], [3], and [49], 
[50] and [19], do not take into account the specific (null) structure of the equations. Yet 
the presence of such structure was crucial in the derivation of the optimal results men- 
tioned above, for geometric semilinear equations. In the case of the Einstein equations it 
is not at all clear what such structure should be, if there is one at all. Indeed, the only 
specific structural condition, known for (EVE), discovered in [30] under the name of the 
weak null condition^ is not at all adequate for improved well posedness results, sec remark 
1.3. It is known however, see [29], that without such a structure one cannot have well 
posedeness for exponents^ s < 2. Yet (EVE) are of fundamental importance and as such 
it is not unreasonable to expect that such a structure must exist. 

Even assuming such a structure, a result of wcll-poscdness for the Einstein equations 
at, or near, the critical regularity Sc = 3/2 is not only completely out of reach but may in 
fact be wrong. This is due to the presence of a different scaling connected to the geometry 
of boundaries of causal domains. It is because of this more subtle scaling that we need 
at least L^-bounds for the curvature to derive a lower bound on the radius of injectivity 
of null hypersurfaces and thus control their local regularity properties. This imposes a 
crucial obstacle to well posedness below s = 2. Indeed, as we will show in the next 
subsection, any such result would require, crucially, bilinear and even trillinear estimates 
for solutions to wave equations of the form Dg^ = F. Such estimates, however, depend 
on Fourier integral representations, with a phase function u which solves the eikonal 
equation g^^daudpu = 0. Thus the much needed bilinear estimates depend, ultimately, 
on the regularity properties of the level hypersurfaces of the phase u which are, of course, 
null. The catastrophic breakdown of the regularity of these null hypersurfaces, in the 
absence of a lower bound for the injectivity radius, would make these Fourier integral 
representations entirely useless. 

These considerations lead one to conclude that, the following conjecture, proposed 
in [18], is most probably sharp in so far as the minimal number of derivatives in is 
concerned: 

Conjecture [Bounded Curvature Conjecture (BCC)] The Einstein- vacuum equa- 
tions admit local Cauchy developments for initial data sets {120, g,k) with locally finite 



'Note that the dimension here is n = 3. 
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curvature and locally finite norm of the first covariant derivatives of k'^. 



Remark 1.2. It is important to emphasize here that the conjecture should be primarily 
interpreted as a continuation argument for the Einstein equations; that is the space-time 
constructed by evolution from smooth data can be smoothly continued, together with a time 
foliation, as long as the curvature of the foliation and the first covariant derivatives of its 
second fundamental form remain L^- bounded on the leaves of the foliation. In particular 
the conjecture implies the break-down criterion previously obtained in [26] and improved 
in [SI], [52]. According to that criterion a vacuum space-time, endowed with a constant 
mean curvature (CMC) foliation Ht, can be extended, together with the foliation, as long 
as the LlL°°{'Et) norm of the deformation tensor of the future unit normal to the foliation 
remains bounded. It is straightforward to see, by standard energy estimates, that this 
condition implies bounds for the Lf^L'^iTit) norm of the space-time curvature from which 
one can derive bounds for the induced curvature tensor R and the first derivatives of the 
second fundamental form k. Thus, if we can ensure that the time of existence of a space- 
time foliated by depends only on the norms of R and first covariant derivatives of 
k, we can extend the space-time indefinitely. 

1.1.3. Brief history. The conjecture has its roots in the remarkable developments 
of the last twenty years centered around the issue of optimal well-posedness for semilinear 
wave equations. The case of the Einstein equations turns out to be a lot more complicated 
due to the quasilincar character of the equations. To make the discussion more tangible it 
is worthwhile to recall the form of the Einstein vacuum equations in the wave gauge. As- 
suming given coordinates x", verifying DgO;" = 0, the metric coefficients ga/3 = g{da,dj3), 
with respect to these coordinates, verify the system of quasilinear wave equations, 



where are quadratic functions of dg, i.e. the derivatives of g with respect to the 
coordinates x". In a first approximation we may compare (1.3) with the semilinear wave 
equation. 



with F quadratic in d(f). Using standard energy estimates, one can prove an estimate, 
roughly, of the form: 

\\m\s < 110(0) II. exp (^CsJ^ \mT)\\L^dr^ . 

The classical exponent s > 3/2 + 1 arises simply from the Sobolev embedding of H^, 
r > 3/2 into L°°. To go beyond the classical exponent, see [34], one has to replace 



g^'^df.d^^ga^ = F^p{g,dg) 




□0 = ^(0,90) 



(1.4) 



we shall see, from the precise theorem stated below, other weaker conditions, such as a lower 
bound on the volume radius, are needed. 
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Sobolev inequalities with Strichartz estimates of, roughly, the following type, 

(^J\mr)\\l^dT^ ' <c(^||a0(O)b.+e + ^*||n0(r)bi+.) 

where £ > can be chosen arbitrarily small. This leads to a gain of 1/2 derivatives, i.e. 
we can prove well-posedncss for equations of type (1.4) for any exponent s > 2. 

The same type of improvement in the case of quasilinear equations requires a highly 
non-trivial extension of such estimates for wave operators with non-smooth coefficients. 
The first improved regularity results for quasilinear wave equations of the type, 

g^\4>)d^d,4> = F{<^M) (1-5) 

with g'^'^{(f)) a non-linear perturbation of the Minkowski metric m'^^, are due to [2], [3], 
and [49], [50] and [19]. The best known results for equations of type (1.3) were obtained 
in [20] and [36]. According to them one can lower the Sobolev exponent s > 5/2 in 
Theorem 1.1 to s > 2. It turns out, see [29], that these results are sharp in the general 
class of quasilinear wave equations of type (1.3). To do better one needs to take into 
account the special structure of the Einstein equations and rely on a class of estimates 
which go beyond Strichartz, namely the so called bilinear estimates^. 

In the case of semilinear wave equations, such as Wave Maps, Maxwell-Klcin-Gordon 
and Yang-Mills, the first results which make use of bilinear estimates go back to [15], 
[16], [17]. In the particular case of the Maxwell- Klein- Gordon and Yang-Mills equation 
the main observation was that, after the choice of a special gauge (Coulomb gauge), the 
most dangerous nonlinear terms exhibit a special, null structure for which one can apply 
the bilinear estimates derived in [15]. With the help of these estimates one was able 
to derive a well posedness result, in the fiat Minkowski space ]R^+'^, for the exponent 
s = Sc + 1/2 = 1, where Sc = 1/2 is the critical Sobolev exponent in that case^. 

To carry out a similar program in the case of the Einstein equations one would need, 
at the very least, the following crucial ingredients: 



A. Provide a coordinate condition, relative to which the Einstein vacuum equations 
verifies an appropriate version of the null condition. 

B. Provide an appropriate geometric framework for deriving bilinear estimates for 
the null quadratic terms appearing in the previous step. 

C. Construct an effective progressive wave representation (parametrix) for solu- 
tions to the scalar linear wave equation Dg^ = F , derive appropriate hounds for 
both the parametrix and the corresponding error term E = F — Dg^F and use 
them to derive the desired bilinear estimates. 



"Note that no such result, i.e. weh-posedness for s = 2, is presently known for either scalar equations 
of the form (1.5) or systems of the form (1.3). 

^This corresponds precisely to the s = 2 exponent in the case of the Einstein-vacuum equations 
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As it turns out, the proof of several bilinear estimates of Step B reduces to the proof 
of sharp L'^{M.) Strichartz estimates for a localized version of the parametrix of step C. 
Thus we will also need the following fourth ingredient. 

D. Prove sharp L'^{A4) Strichartz estimates for a localized version of the parametrix 
of step C. 

Note that the last three steps need to be implemented using only hypothetical 
bounds for the space-time curvature tensor, consistent with the conjectured result. To 
start with, it is not at all clear what should be the correct coordinate condition, or even 

if there is one for that matter. 

Remark 1.3. As mentioned above, the only known structural condition related to 
the classical null condition, called the weak null condition, tied to wave coordinates, fails 
the test. Indeed, the following simple system in Minkowski space verifies the weak null 
condition and yet, according to [29], it is ill posed for s — 2. 

□ = 0, □'0 = 0- A0. 

Coordinate conditions, such as spatial harmonic^^, also do not seem to work. 

We rely instead on a Coulomb type condition, for orthonormal frames, adapted to a 
maximal foliation. Such a gauge condition appears naturally if we adopt a Yang-Mills 
description of the Einstein field equations using Cartan's formalism of moving frames^^, 
see [6]. It is important to note nevertheless that it is not at all a priori clear that such 
a choice would do the job. Indeed, the null form nature of the Yang-Mills equations in 
the Coulomb gauge is only revealed once we commute the resulting equations with the 
projection operator V on the divergence free vectorfields. Such an operation is natural 
in that case, since V commutes with the flat d'Alembertian. In the case of the Einstein 
equations, however, the corresponding commutator term [Dg,?^] generates^^ a whole host 
of new terms and it is quite a miracle that they can all be treated by an extended version 
of bihnear estimates. At an even more fundamental level, the flat Yang-Mills equations 
possess natural energy estimates based on the time symmetry of the Minkowski space. 
There are no such timelike Killing vectorfleld in curved space. We have to rely instead 
on the future unit normal to the maximal foliation whose deformation tensor is non- 
trivial. This leads to another class of nonhnear terms which have to be treated by a novel 
trilinear estimate. 

We will make more comments concerning the implementations of all four ingredients 
later on, in the section 1.2.4. 

^"^Maximal foliation together with spatial harmonic coordinates on the leaves of the foliation would 
be the coordinate condition closest in spirit to the Coulomb gauge. 

^^We would like to thank L. Anderson for pointing out to us the possibility of using such a formalism 
as a potential bridge to [16] . 

''^^Note also that additional error terms are generated by projecting the equations on the components 
of the frame. 
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Remark 1.4. In addition to the ingredients mentioned above, we also need a mecha- 
nism of reducing the proof of the conjecture to small data, in an appropriate sense. Indeed, 
even in the flat case, the Coulomb gauge condition cannot be globally imposed for large 
data. In fact [17] relied on a cumbersome technical device based on local Coulomb gauges, 
defined on domain of dependence of small balls. Here we rely instead on a variant of the 
gluing construction o/[10], [11], see section 1.2.3. 

1.2. Statement of the main results 

1.2.1. Maximal foliations. In this section, we recall some well-known facts about 
maximal foliations (see for example the introduction in [9]). We assume the space-time 

g) to be foliated by the level surfaces of a time function t. Let T denote the unit 
normal to E^, and let k the the second fundamental form of E^, i.e. kah — — g(Do2^) Cf,), 
where e^, a = 1, 2, 3 denotes an arbitrary frame on E^ and = Dea^- We assume that 
the Ej foliation is maximal, i.e. we have: 

tr^A; = (1.6) 

where g is the induced metric on Ej. The constraint equations on E^ for a maximal 
foliation are given by: 

V'kab = 0, (1.7) 

where V denotes the induced covariant derivative on Ej, and 

Rscai = (1.8) 

Also, we denote by n the lapse of the i-fohation, i.e. = — g(Dt, Dt). n satisfies the 
following elliptic equation on E^: 

An = n|A;|l (1.9) 
Finally, we recall the structure equations of the maximal foliation: 

Vokab = RaObO -n^VaVfen- kackb", (1.10) 

Va/Cbc - Vfc/Cac = RcOa& (l-H) 

and: 

Rab — kack'^ b — ^aObO- (1-12) 

1.2.2. Main Theorem. We recall below the definition of the volume radius on a 
general Riemannian manifold M. 

Definition 1.5. Let Brip) denote the geodesic ball of center p and radius r. The 
volume radius ryoi{p, r) at a point p & M and scales <r is defined by 

\Mp)\ 



rvoi{p,r) = inf 



r'<r r"^ 
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with \Br\ the volume of Br relative to the metric on M. The volume radius ryoi{M,r) of 
M on scales < r is the infimum ofr^oiip, r) over all points p e M. 

Our main result is the following: 

Theorem 1.6 (Main theorem). Let (A^,g) an asymptotically fiat solution to the Ein- 
stein vacuum equations (1.1) together with a maximal foliation by space-like hypersur- 
faces Ht defined as level hypersurfaces of a time function t. Assume that the initial slice 
(T,Q,g,k) is such that the Ricci curvature Ric G L^(So), VA; G L^(So), and Eq has a 
strictly positive volume radius on scales < 1, i.e. r^o/(So, 1) > 0. Then, 

(1) regularity. There exists a time 

T = T(||/?zc||i2(s„), ||V/c|U2(So),r,„KSo, 1)) > 
and a constant 

C = C(||i?zc|U2(s„), II VA:|U2(s„), r,,z(So, 1)) > 
such that the following control holds on <t <T: 

||R||l-^,l2(so < C, ||VA;||Lcx,^jL2(Et) < C and ^mf^r^o,(Et, 1) > ^. 

(2) Higher regularity. Within the same time interval as in part (1) we also have 
the higher derivative estimates^^, 



|a|<m Kl<"* 

where Cm depends only on the previous C and m. 



:i.i3) 



Remark 1.7. Since the core of the main theorem is local in nature we do not need to 
he very precise here with our asymptotic flatness assumption. We may thus assume the 
existence of a coordinate system at infinity, relative to which the metric has two derivatives 
hounded in , with appropriate asymptotic decay. Note that such hounds could he deduced 
from weighted bounds assumptions for Ric and Vk. 

Remark 1.8. Note that the dependence on ||i?zc||/^2(j](,), || VA';||i2|-S(,) in the main theo- 
rem can be replaced by dependence on ||R||l2(So) where R denotes the space-time curvature 
tensor^'^. Indeed this follows from the following well known estimate (see section 8 in 

[26];. 

/ \Vk\'' ^-M" < [ |Rp. (1.14) 

and the Gauss equation relating Ric to R. 

^■^ Assuming that the initial has more regularity so that the right-hand side of (1.13) makes sense. 
"'^^Here and in what follows the notations R, R will stand for the Riemann curvature tensors of 
and A4, while Ric, Ric and Rscah^scai will denote the corresponding Ricci and scalar curvatures. 
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1.2.3. Reduction to small initial data. We first need an appropriate covering of 
Eo by harmonic coordinates. This is obtained using the following general result based on 
Cheeger-Gromov convergence of Riemannian manifolds. 

Theorem 1.9 ([1] or Theorem 5.4 in [32]). Given ci > 0, C2 > 0, cs > 0, there 

exists Tq > such that any 3- dimensional, complete, Riemannian manifold {M,g) with 
\\R'i'(^\\L'^(M) ^ ci and volume radius at scales < 1 bounded from below by C2, i.e. ryoi{M, 1) > 
02, verifies the following property: 

Every geodesic ball Br{p) with p & M and r < r^ admits a system of harmonic 
coordinates x — {xi,X2,xs) relative to which we have 

(1 + cs)-'5ij < gij < (1 + C3)5ij, (1.15) 

and 

r [ \^^gij\^^Mdx<C3■ (1-16) 
Jbm 



We consider e > which will be chosen as a small universal constant. We apply 
theorem 1.9 to the Riemannian manifold Sq. Then, there exists a constant: 

ro = ro(||Ric||L2(2Q), \\Vk\\L2(^^^),r^oii^o, 1),^) > 

such that every geodesic ball Br{p) with p G Sq and r < ro admits a system of harmonic 
coordinates x = {xi,X2,Xs) relative to which we have: 

{l + e)-^6ij <gij < {l + e)Sij, 

and 



Ib, 



'Br{p) 

Now, by the asymptotic flatness of Sq, the complement of its end can be covered 
by the union of a finite number of geodesic balls of radius tq, where the number Nq of 
geodesic balls required only depends on ro- In particular, it is therefore enough to obtain 
the control of R, k and r^oii^t, 1) of Theorem 1.6 when one restricts to the domain of 
dependence of one such ball. Let us denote this ball by B^.^. Next, we rescale the metric 
of this geodesic ball by: 

( \ 

9x{t,x) = g{Xt,Xx), A = min j-^ , ^ , roe > 0. 

\\\^\\L^Br,) W^f'WL^Br,) J 

Let^^ Rx, k\ and B^^ be the rescaled versions of R, k and B^.^. Then, in view of our choice 
for A, we have: 

\\R\\\l^{B>^^) = ^\\R\\L-^{Br^) < £, 



""■^Since in what follows there is no danger to confuse the Ricci curvature Ric with the scalar curvature 
R we use the short hand R to denote the full curvature tensor Ric. 



10 1. INTRODUCTION 

and 

Note that B^^ is the rescaled version of B,.^. Thus, it is a geodesic ball for gx of radius 

^ > ^ > 1. Now, considering on < t < 1 is equivalent to considering 5^ on < t < A. 
Thus, since Tq, Nq and A depend only on ||-R||l2(£p), || VA;||i:2(S(,), r^o;(Eo, 1) and £, Theorem 
1.6 is equivalent to the following theorem: 

Theorem 1.10 (Main theorem, version 2). Let (A^,g) an asymptotically flat solution 
to the Einstein vacuum equations (1.1) together with a maximal foliation by space-like 
hypersurfaces Tit defined as level hypersurfaces of a time function t. Let B a geodesic ball 
of radius one in Sq, and let D its domain of dependence. Assume that the initial slice 
(Eo, g, k) is such that: 

11^11x2(5) < \\^k\\L'HB) < £ andryoi{B,l) > ^. 

Let Bf = D nHf the slice of D at time t. Then: 

(1) regularity. There exists a small universal constant Eq > Q such that if 
Q < £ < £q, then the following control holds on Q <t <1: 

||R|U-j,L2(B,) < £, ||V/i;||L-^,L2(B0 < e and M^^r^oiiBt, 1) > |. 

(2) Higher regularity. The following bounds hold on <t < 1: 

J2 I|D(")R||ljs=,,,l2(bo < l|V«i^^c|U2(B) + ||V»V^|U2(s). (1.17) 

Notation: In the statement of Theorem 1.10, and in the rest of the paper, the notation 
/i ^ h for two real positive scalars /i,/2 means that there exists a universal constant 
C > such that: 

fi < C/2. 



Theorem 1.10 is not yet in a suitable form for our proof since some of our constructions 
will be global in space and may not be carried out on a subregion B of Sq. Thus, we 
glue a smooth asymptotically fiat solution of the constraint equations (1.2) outside of B, 
where the gluing takes place in an annulus just outside B. This can be achieved using the 
construction in [10], [11]. We finally get an asymptotically flat solution to the constraint 
equations, deflned everywhere on Sq, which agrees with our original data set {T,Q,g,k) 
inside B. We still denote this data set by (Eq, g, k). It satisfies the bounds: 

||i?||L2(So) < 2£, ||VA;||l2(So) < 2e and r,<,^(Eo, 1) > ^. 
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Remark 1.11. Notice that the gluing process in [10] -[11] requires the kernel of a 
certain linearized operator to be trivial. This is achieved by conveniently choosing the 
asymptotically flat solution to (1.2) that is glued outside of B to our original data set. 
This choice is always possible since the metrics for which the kernel is nontrivial are non 
generic (see [4]j. 

Remark 1.12. Assuming only bounds on R and VA; is not enough to carry out the 
construction in the above mentioned results. However, the problem solved there remains 
subcritical at our desired level of regularity and thus we believe that a closer look at the 
construction in [10] -[ll], or an alternative construction, should be able to provide the 
desired result. This is an open problem. 

Remark 1.13. Since ||fc|||4(So) < ||-R«c||i2 we deduce that ||A;||l2(b) < e^/^ on the 
geodesic ball B of radius one. Furthermore, asymptotic flatness is compatible with a decay 
of |a;|~^ at infinity, and in particular with k in L^(Eo). So we may assume that the gluing 
process is such that the resulting k satisfies: 

\\k\\L'\ji,,) < e. 

Finally, we have reduced Theorem 1.6 to the case of a small initial data set: 

Theorem 1.14 (Main theorem, version 3). Let (A^,g) an asymptotically flat solution 
to the Einstein vacuum equations (1.1) together with a maximal foliation by space-like 
hypersurfaces defined as level hypersurfaces of a time function t. Assume that the 
initial slice (So, 9', k) is such that: 

||-R||l2(So) < £, ||^||l2(So) + ||VA;||l2(S(,) < e and r^oii^o, 1) > ^■ 

Then: 

(1) regularity. There exists a small universal constant £o > such that if 
< e < So, the following control holds on < t < 1: 

||R|U^yL2(so < e, ||A;||L-yL2(s,) + ||VA;||l-^jL2(s,) < e and ^M^r^oii^tA) > \- 

(2) Higher regularity. The following control holds on <t <1: 

l|D(")R|Ufj_^,L2(E0 < E l|V^')i?^c|U2(So) + ||V»VA;|U2(s„). (1.18) 

|a|<m |i|<m 

The rest of this paper is devoted to the proof of Theorem 1.14. 

1.2.4. Strategy of the proof. The proof of Theorem 1.14 consists of four steps. 

Step A (Yang-Mills formalism) We first cast the Einstein-vacuum equations in a 
Yang-Mills form. This relies on the Cartan formalism of moving frames. The idea is to 
give up on a choice of coordinates and instead express the Einstein vacuum equations in 
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terms of the connection 1-forms associated to moving orthonormal frames, i.e. vectorfields 
Ca-i which verify, 

g(e„, ep) = m^is = diag(-l, 1, 1, 1). 

The connection 1-forms (they are to be interpreted as 1-forms with respect to the external 
index /i with values in the Lie algebra of so(3, 1)), defined by the formulas, 

(A^)a/3 = g(D/,e^,e„) (1.19) 

verify the equations, 

D'^F^, + [A^F^,] = (1.20) 
where, denoting (F^^)„/3 := Ra^^u, 

(F,,,)„,fl = (D^A„ - D,.A^ - [A^, A,])^^. (1.21) 

In other words we can interpret the curvature tensor as the curvature of the so(3, l)-valucd 
connection 1-form A. Note also that the covariant derivatives are taken only with respect 
to the external indices fj,, v and do not affect the internal indices a, 13. We can rewrite 
(1.20) in the form, 

□gA, - D,(D'^A^) = J,(A, DA) (1.22) 

where, 

J, = D'^([A^,A,])-[A^,F^,]. 

Observe that the equations (1.20)-(1.21) look just like the Yang-Mills equations on a 
fixed Lorentzian manifold (A1,g) except, of course, that in our case A and g are not 
independent but rather connected by (1.19), reflecting the quasilinear structure of the 
Einstein equations. Just as in the case of [15], which establishes the well-posedness of 
the Yang-Mills equation in Minkowski space in the energy norm (i.e. s = 1), wc rely in 
an essential manner on a Coulomb type gauge condition. More precisely, wc take Cq to 
be the future unit normal to the foliation and choose 61,62,63 an orthonormal basis 
to Et, in such a way that we have, essentially (sec precise discussion in section 2.1.2), 
div A — VMj = 0, where A is the spatial component of A. It turns out that Aq satisfies 
an elliptic equation while each component Ai — g(A,ej), i — 1,2,3 verifies an equation 
of the form, 

□gA = -di{doAo)+A^djAi + A^diAj+ l.o.t. (1.23) 

with l.o.t. denoting nonlinear terms which can be treated by more elementary techniques 
(including non sharp Strichartz estimates). 

Step B (Bilineeir and trilineeir estimates) To ehminate 5j(9oAo) in (1.23), we need 
to project (1.23) onto divergence free vectorfields with the help of a non-local operator 
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which we denote by V. In the case of the flat Yang- Mills equations, treated in [15], this 
leads to an equation of the form, 

UAi = V{AWjAi) +r{A^diAj) + l.o.t. 

where both terms on the right can be handled by bilinear estimates. In our case we 
encounter however three fundamental differences with the flat situation of [15]. 

• To start with the operator V does not commute with Dg. It turns out, fortunately, 
that the terms generated by commutation can still be estimated by an extended 
class of bilinear estimates which includes contractions with the curvature tensor, 
see section 2.2.5. 

• All energy estimates used in [15] are based on the standard timelike Killing 
vectorfield dt- In our case the corresponding vectorfield cq — T { the future unit 
normal to E^) is not Killing. This leads to another class of trilinear error terms 
which we discuss in sections 2.5 and 2.2.5. 

• The main difference with [15] is that we now need bilinear and trilinear estimates 
for solutions of wave equations on background metrics which possess only limited 
regularity. 

This last item is a major problem, both conceptually and technically. On the conceptual 
side we need to rely on a more geometric proof of bilinear estimates based on a plane 
wave representation formula^^ for solutions of scalar wave equations, 

□g0 = 0. 

The proof of the bilinear estimates rests on the representation formula^^ 

p poo 

(t)f{t,x)= / / e'^""(*'^)/(Aa;)A'dAda; (1.24) 

i§2 Jo 

where / represents schematically the initial data^^, and where '^u is a solution of the 
eikonal equation^^, 

g^^^a^ '^'ua^ "^M- 0, (1.25) 

with appropriate initial conditions on Eq and dui the area element of the standard sphere 
in M?. 



We follow the proof of the bilinear estimates outlined in [21] which differs substantially from that 
of [15] and is reminiscent of the null frame space strategy used by Tataru in his fundamental paper [48]. 

^^(1.24) actually corresponds to the representation formula for a half-wave. The full representation 
formula corresponds to the smn of two half- waves (see section 2.7) 

^^Hcrc / is in fact at the level of the Fourier transform of the initial data and the norm ||A/||j^2(g3') 
corresponds, roughly, to the i?^ norm of the data . 

■"^^In the flat Minkowski space "u{t, x) =t±x ■ lj. 
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Remark 1.15. Note that (1-24) is a parametrix for a scalar wave equation. The lack 
of a good parametrix for a covariant wave equation forces us to develop a strategy based on 
writing the main equation in components relative to a frame, i. e. instead of dealing with 
the tensorial wave equation (1.22) directly, we consider the system of scalar wave equations 
(1.23). Unlike the fiat case, this "scalarization" procedure produces several terms which 
are potentially dangerous, and it is fortunate, as in yet another manifestation of a hidden 
null structure of the Einstein equations, that they can still be controlled by the use of an 
extende(f^ class of bilinear estimates. 

Step C (Control of the parametrix) To prove the bihnear and trihnear estimates of 
Step B, we need in particular to control the parametrix at initial time (i.e. restricted to 
the initial slice Eq) 

0^(0, a;)=/ / e*^''"('^'^)/(Aa;)A'ciAcicc; (1.26) 
Jo 

and the error term corresponding to (1.24) 

Ef{t,x) = □g0/(t,a;) =i / e^^""(*'^) (Dg ^u)f{\uj)\^dXduj (1.27) 

i.e. 0j is an exact solution of Dg^ = only in flat space in which case Dg '^u = 0. This 
requires the following four sub steps 

CI Make an appropriate choice for the equation satisfied by '^u{Q,x) on Eg, and 

control the geometry of the foliation ofEo by the level surfaces of '^u{0,x). 
C2 Prove that the parametrix at t = given by (1.26) is bounded in C{L'^ {E?) , LP' (Tiq)) 

using the estimates for ^u{Q,x) obtained in CI. 
C3 Control the geometry of the foliation of Ai given by the level hypersurfaces of 

C4 Prove that the error term (1.27) satisfies the estimate ||-E'/||i2(_y\4) < C||A/||i2(]R3) 
using the estimates for and Dg '^u proved in C3. 

To achieve Step C3 and Step C4, we need, at the very least, to control Dg '^u in 
This issue was first addressed in the sequence of papers [22]-[24] where an L°° bound 
for Dg '^li was established, depending only on the norm of the curvature flux along 
null hypersurfaces. The proof required an interplay between both geometric and analytic 
techniques and had all the appearances of being sharp, i.e. we don't expect an bound 
for Dg '^u which requires bounds on less than two derivatives in for the metric^^. 

To obtain the bound for the Fourier integral operator E deflned in (1.27), we need, 
of course, to go beyond uniform estimates for Dg '^u. The classical bounds for Fourier 
integral operators of the form (1.27) are not at all economical in terms of the number of 
integration by parts which are needed. In our case the total number of such integration 

^'^involving contractions between the Riemann curvature tensor and derivatives of solutions of scalar 
wave equations. 

^-"^classically, this requires, at the very least, the control of R in 
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by parts is limited by the regularity properties of the function Dg '^u. To get an L bound 
for the parametrix at initial time (1-26) and the error term (1.27) within such restrictive 
regularity properties we need, in particular: 

• In Step CI and Step C3, a precise control of derivatives of '^u and Dg '^u with 
respect to both u as well as with respect to various directional derivatives^^. To 
get optimal control we need, in particular, a very careful construction of the initial 
condition for '^u on Eq and then sharp space-time estimates of Ricci coefficients, 
and their derivatives, associated to the foliation induced by '^u. 

• In Step C2 and Step C4, a careful decompositions of the Fourier integral oper- 
ators (1.26) and (1.27) in both A and uj, similar to the first and second dyadic 
decomposition in harmonic analysis, see [39], as well as a third decomposition, 
which in the case of (1.27) is done with respect to the space-time variables relying 
on the geometric Littlewood-Paley theory developed in [24]. 

Below, we make further comments on Steps C1-C4: 

(1) The choice of u{0,x,u) on Eg in Step CI. Let us note that the typical choice 
u{0, x,uj) = X ■ CO in a given coordinate system would not work for us, since we 
don't have enough control on the regularity of a given coordinate system within 
our framework. Instead, we need to find a geometric definition of u{0,x,uj). A 
natural choice would be 

HgU = on Eo 

which by a simple computation turns out to be the following simple variant of 
the minimal surface equation^^ 

div ( I = k I , Y , , , Y , 1 on Er 



Unfortunately, this choice docs not allow us to have enough control of the deriva- 
tives of u in the normal direction to the level surfaces of u. This forces us to look 
for an alternate equation for u: 

^. f Vu\ ^ 1 , / Vm Vm \ 
div —— = 1 - + k ——, —— on Eq. 
\|V-u|/ \\/u\ \\yu\ \Vu\J 

This equation turns out to be parabolic in the normal direction to the level 
surfaces of u, and allows us to obtain the desired regularity in Step CI. A closer 
inspection reveals its relation to the mean curvature flow on Eq. 
(2) How to achieve Step C3. The regularity obtained in Step CI, together with 
null transport equations tied to the eikonal equation, elliptic systems of Hodge 
type, the geometric Littlewood-Paley theory of [24], sharp trace theorems, and 

^^Taking into account the different behavior in tangential and transversal directions with respect to 
the level surfaces of '^u. 

^^In the time symmetric case k = 0, this is exactly the minimal surface equation 
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an extensive use of the structure of the Einstein equations, allows us to propagate 
the regularity on Eq to the space-time, thus achieving Step C3. 

(3) The regularity with respect to ui in Steps CI and C3. The regularity with respect 
to X for u is clearly limited as a consequence of the fact that we only assume 
bounds on R. On the other hand, R is independent of the parameter cj, and one 
might infer that u is smooth with respect to ui. Surprisingly, this is not at all 
the case. Indeed, the regularity in x obtained for u in Steps CI and C3 is better 
in directions tangent to the level hypersurfaces of u. Now, the u derivatives of 
the tangential directions have non zero normal components. Thus, when differ- 
entiating the structure equations with respect to a;, tangential derivatives to the 
level surfaces of u are transformed to non tangential derivatives which in turn 
severely limits the regularity in oj obtained in Steps CI and C3. 

(4) How to achieve Steps C2 and C4- The classical arguments for proving bounds 
for Fourier operators are based either on a TT* argument, or a T*T argument, 
which requires several integration by parts either with respect to x for T*T, 
of with respect to (A.cj) for TT*. Both methods would fail by far within the 
regularity for u obtained in Step CI and Step C3. This forces us to design a 
method which allows to take advantage both of the regularity in x and u. This 
is achieved using in particular the following ingredients: 

• geometric integrations by parts taking full advantage of the better regularity 
properties in directions tangent to the level hypersurfaces of u, 

• the standard first and second dyadic decomposition in frequency space, with 
respect to both size and angle (see [39]), an additional decomposition in 
physical space relying on the geometric Littlewood-Paley projections of [24] 
for Step C4, as well as another decomposition involving frequency and angle 
for Step C2. 

Even with these precautions, at several places in the proof, one encounters log- 
divergences which have to be tackled by ad-hoc techniques, taking full advantage 
of the structure of the Einstein equations. 

Step D (Sharp L'^{Ai) Strichartz estimates) Recall that the parametrix constructed 
in Step C needs also to be used to prove sharp L^(A^) Strichartz estimates. Indeed 
the proof of several bilinear estimates of Step B reduces to the proof of sharp L^{A4) 
Strichartz estimates for the parametrix (1.24) with A locahzed in a dyadic shell. 
More precisely, let j > 0, and let t/j a smooth function on supported in 



Let (pfj the parametrix (1.24) with a additional frequency localization A ~ 2^ 




(1.28) 
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We will need the sharp^^ L^{Ai) Strichartz estimate 

UfjLHM)<2'-W2-'X)f\\LHRSy (1.29) 

The standard procedure for proving^^ (1-29) is based on a TT* argument which reduces 
it to an L°° estimate for an oscillatory integral with a phase involving '^u. This is then 
achieved by the method of stationary phase which requires quite a few integrations by 
parts. In fact the standard argument would require, at the very least^^, that the phase 
function u — '^u verifies, 

dt,.u e L^, dt,,dlu e L^. (1.30) 

This level of regularity is, unfortunately, incompatible with the regularity properties of 
solutions to our eikonal equation (1.25). In fact, based on the estimates for '^u derived 
in step C3, we are only allowed to assume 

dt,,u e L°°, dt,,,d^u e L°°. (1.31) 

We are thus forced to follow an alternative approach^^ to the stationary phase method 
inspired by [35] and [36] . 

Remark 1.16. Note that apart from the results of Chapter 2 which require the pro- 
jection of various tensors on a frame, the computations and estimates in all the other 

chapters are covariant. 

1.2.5. Structure of the paper. In Chapter 2, we perform Step A and Step B, 
i.e. we recast the Einstein equations as a quasilinear Yang-Mills type system, we prove 
bilinear estimates, and we reduce the proof of Theorem 1.14 to Step C and Step D. Next, 
wc perform Step C on the control of the plane wave parametrix (1.24). More precisely, in 
Chapter 3, we perform Step C4 on the control of the error term (1.27). Next, in Chapter 
4, we perform Step C3 on the space-time control of the optical function ^u. Then, in 
Chapter 5, we perform Step C2 on the control of the parametrix at initial time (1.26). In 
Chapter 6, we perform Step CI on the control of the optical function on the initial 
slice Eq. Finally, in Chapter 7, wc prove sharp L^(A4) Strichartz estimates localized in 
frequency which corresponds to Step D. 

Remark 1.17. Chapter 2 summarizes the results obtained in [27]. Chapter 3 sum- 
marizes the results obtained in [45]. Chapter 4 summarizes the results obtained in [44]. 
Chapter 5 summarizes the results obtained in [43]. Chapter 6 summarizes the results 
obtained in [42]. Finally, Chapter 7, summarizes the results obtained in [46]. 

^^Note in particular tliat tlie corresponding estimate in the flat case is sharp. 

^^Note that the procedure we describe would prove not only (1.29) but the full range of mixed 
Strichartz estimates. 

^^The regularity (1.30) is necessary to make sense of the change of variables involved in the stationary 
phase method. 

^'^Wc refer to the approach based on the overlap estimates for wave packets derived in [35] and [36] 
in the context of Strichartz estimates respectively for C^'^ and metrics. Note however that our 

approach does not require a wave packet decomposition. 



18 



1. INTRODUCTION 



Remark 1.18. The structure of this overview is such that each part motivates the next 
one. In particular, Chapter 2 relies on the control of the parametrix (1.24) (see Theorem 
2.27 in Chapter 2), and thus motivates Chapters 3, 4, 5 and 6 which precisely deal with 
the control of that parametrix. Next, in order to control the error term (1.27) in Chapter 
3, we rely on estimates for the optical function '^u, which motivates Chapter 4 where 
these estimates are proved. In turn, the space-time estimates for '^u in Chapter 4 o,re 
obtained in particular using transport equations, and we need the corresponding control 
for '^u on the initial slice Sq, which motivates Chapter 6. Finally, in order to control the 
parametrix at initial time (1-26) in Chapter 5, we rely on estimates for the function '^u 
on Eo; which motivates again Chapter 6. Finally, Chapter 2 also relies on sharp L'^{M.) 
Strichartz estimates localized in frequency (see Proposition 2.32 in Chapter 2), and thus 
motivates Chapter 7. 

1.2.6. Conclusion. Though this result does not achieve the crucial goal of finding 
a scale invariant well-posedness criterion in GR, it is clearly optimal in terms of all cur- 
rently available ideas and techniques. Indeed, within our current understanding, a better 
result would require enhanced bilinear estimates, which in turn would rely heavily on 
parametrices. On the other hand, parametrices are based on solutions to the eikonal 
equation whose control requires, at least, bounds for the curvature tensor, as can be 
seen in many instances in our work. Thus, if we are to ultimately find a scale invariant 
well-posedness criterion, it is clear that an entirely new circle of ideas is needed. Such 
a goal is clearly of fundamental importance not just to GR, but also to any physically 
relevant quasilinear hyperbolic system. 

1.2.7. Acknowledgements. This work would be inconceivable without the extraor- 
dinary advancements made on nonlinear wave equations in the last twenty years in which 
so many have participated. Wc would like to single out the contributions of those who 
have affected this work in a more direct fashion, either through their papers or through 
relevant discussions, in various stages of its long gestation. D. Ghristodoulou's seminal 
work [8] on the weak cosmic censorship conjecture had a direct motivating role on our 
program, starting with a series of papers of the first author and M. Machcdon, in which 
spacetime bilinear estimates were first introduced and used to take advantage of the null 
structure of geometric semilinear equations such as Wave Maps and Yang-Mills. The 
works of Bahouri- Ghemin [2]- [3] and D.Tataru [50] were the first to go below the clas- 
sical Sobolev exponent s = 5/2, for any quasilinear system in higher dimensions. This 
was, at the time, a major psychological and technical breakthrough which opened the 
way for future developments. Another major breakthrough of the period, with direct 
influence on our approach to bilinear estimates in curved spacetimes, is D. Tataru's work 
[48] on critical well posedness for Wave Maps, in which null frame spaces were first in- 
troduced. His joint work with H. Smith [36] which, together with [20] is the first to 
reach optimal well-posedness without bilinear estimates, has also infiuenced our approach 
on parametrices and Strichartz estimates. The authors would also like to acknowledge 
fruitful conversations with L. Anderson, and J. Sterbenz. 



CHAPTER 2 

Einstein vacuum equations as a Yang-Mills gauge theory 

Recall Steps A, B, C and D introduced in section 1.2.4. In this chapter, we perform 
Step A and Step B, i.e. we recast the Einstein equations as a quasilinear Yang-Mills type 
system and we prove bilinear estimates. This allows us to reduce the proof of Theorem 
1.14 to Step C and Step D. Here, we only outhne the main ideas, and we refer to [27] for 
the details. 

2.1. Yang-Mills formalism 

2.1.1. Cartan formalism. Consider an Einstein vacuum space-time (A1,g). We 
denote the covariant differentiation by D. Let be an orthonormal frame on A4, i.e. 

g(ea, e^) = T^ap = diag(-l, !,...,!). 

Consistent with the Cartan formalism we define the connection 1 form, 

{A)ais{X)=g(Dxe^,ea) (2.1) 

where X is an arbitrary vectorfield in T{M.). Observe that, 

i.e. the 1 -form A^dx^ takes values in the Lie algebra of 5*0(1,3). We separate the 
internal indices ct, ^ from the external indices ii according to the following notation. 

(A^),^ := {A)^p{d^) = g(D^e^, (2.2) 

The Riemann curvature tensor is defined by R(X, Y, U, V) = g(X, [Df/Dy — DyDj; — 
D[[/y]F] ) with X, y, [/, V arbitrary vectorfields in T{M). Thus, taking U — d^j^.V — d^, 
coordinate vector-fields, 

R(ea, e^, 9^, d^) = ^^(A^)^,^ - 9^(A^)„^ + (A^)„ ^(A^)a^ - (A^)„ ^{A^)xf}. (2.3) 

Defining the Lie bracket, 

([A^, A,]),^ = (A^), T (A,)^^ - (A,), T (A^)^^ (2.4) 

we obtain: 

R-a^/iz^ = dn{A„)afi — di,{A^)ap — ([A^, A„])ap, 

or, since 9^(A^) - 9j.(A^) = D^A^ - D^A^ 

(F^,u)a0 = '^afinu = (D^A^ - D^A^ - [A^, A^])^^. (2.5) 
where interpret F is the curvature of the connection A. 
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The usual covariant derivative of the Riemann curvature tensor can be expressed as 
follows: 

D^R„^^, = (^>D^F^, := D^F^, + [A,, F^,] (2.6) 

where we denote by '•^•t) the covariant derivative on the corresponding vector bundle. 
More precisely if U = U^^^^j ^-tensor on M. with values on the Lie algebra of 

SO{3,l), 

(^>D,U = D,U+[A,,U]. (2.7) 

Remark 2.1. Recall that in (A^)^,^, a,f3 are called the internal indices, while fi are 
called the external indices. Now, the internal indices are mostly irrelevant in our work. 
Thus, from now on, we will drop them, except for rare instances where we will need to 
distinguish between internal indices of the type ij and internal indices of the type Oi. 

The Bianchi identities for Ra/3//i/ take the form 

(^>D.F^, + (^>D^F,, + (^>D.F^^ = 0. (2.8) 

As it is well known the Einstein vacuum equations Ha/s = imply D'^Rq,^^,^ = 0. Thus, 
in view of equation (2.6), 

= (^to'^ F^, = D'^F^, + [A^ F^,] (2.9) 

or, in view of (2.5) and the vanishing of the Ricci curvature of g,. 

□A, - D,(D'^A^) = J, (2.10) 

where 

J, = D''([A^, A,]) - [A^, F^,]. (2.11) 

Using again the vanishing of the Ricci curvature it is easy to check, 

D^J^ = 0. (2.12) 

Finally we recall the general formula of transition between two different orthonormal 
frames and Cq on A4, related by, 

where niap — O^O^ m^^, i.e. O is a smooth map from M. to the Lorentz group 0(3, 1). 
In other words, raising and lowering indices with respect to m, 

O^aO^" = 5^ (2.13) 
Now, (A^)a^ = g(D^e/3,ea). Therefore, 

(A,U = OIOUA,),s + d,{Ol) (2.14) 
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2.1.2. Compatible frames. Recall that our space-time is assumed to be foliated 
by the level surfaces of a time function t, which are maximal, i.e. denoting by k the 
second fundamental form of we have, 

tr^A; = (2.15) 

where g is the induced metric on S^. Let us choose e(o) = T, the future unit normal to 
the Et foliation, and e(j), i = 1,2,3 an orthonormal frame tangent to E^. We call this a 
frame compatible with our Ej fohation. We consider the connection coefficients (2.2) with 
respect to this frame. Thus, in particular, denoting by Aq, respectively A^, the temporal 
and spatial components of 

{Ai)oj = (Aj)oi^-kij, = 1,2,3 (2.16) 

(Ao)oi = -n-^ViU 1,2,3 (2.17) 

where n denotes the lapse of the i-foliation, i.e. = — g(Dt, Dt). With this notation 
we note that, 

Vihj = Vi{ki)j + kin{Ai), " = V\Ai\j + kin{Ai)j'' 

where, as before, the notation 'Wi{ki)j or V'(Aj)oj, is meant to suggest that the covariant 
differentiation affects only the external index i. Recalling from (1.7) that k verifies the 
constraint equations, 

V%j = 0, 

we derive, 

V\Ai)oj^kr{Ai)mj. (2.18) 

Besides the choice of cq we are still free to make a choice for the spatial elements of the 
frame ci, 62, 63. In other words we consider frame transformations which keep Cq fixed, i.e 
transformations of the type, 

with O in the orthogonal group 0(3). We now have, according to (2.14), 

{AJ,j = OlO^jiAJki + d^{0'i)0'^5ki 

or, schematically. 

Am = OAraO'^ + {dmO)0-^ (2.19) 

formula in which we understand that only the spatial internal indices are involved. We 
shall use this freedom later to exhibit a frame ei, 62, 63 such that the corresponding con- 
nection A satisfies the coulomb gauge condition V'(^;)ij = (see Lemma 2.6). 
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2.1.3. Notations. We use greek indices to denote general indices on which do 
not refer to the particular frame (eo, ei, 62, 63). The letters a, b, c, d will be used to denote 
general indices on which do not refer to the particular frame (61,62,63). Finally, the 
letters i,j,l,m,n will only denote indices relative to the frame (61,62,63). Also, recall 
that D denotes the covariant derivative on Ai, while V denotes the induced covariant 
derivative on E^. Furthermore, d will always refer to the derivative of a scalar quantity 
relative to one component of the frame (60,61,62,63), while d will always refer to the 
derivative of a scalar quantity relative to one component of the the frame (61,62,63), so 
that d = (80,8). For example, 8 A may be any term of the form 8i{Aj), 8o{A) may be 
any term of the form 8Q{Aj), 9(^40) may be any term of the form 8j{Ao), and dA — 
{dA,d{Ao)) = (do{Ao),d{Ao),do{A),d{A)). Note that we use brackets such as (Aj) to 
emphasize that we are dealing with su{3,l) objects. Often, however, we will simply drop 
them. 

We introduce the curl operator curl defined for any su{3, l)-valued triplet {cui, 002,003) 
of functions on E^ as follows: 

icurlco)i =e, ''8,itoi), (2.20) 

where 6,^7 is fully antisymmetric and such that €123= 1- We also introduce the divergence 
operator div defined for any su{3, l)-valued tensor A on E^ as follows: 

div A = V'(A) = 8\Ai) + A^. (2.21) 

Remark 2.2. The term A? in (2.21) corresponds to a quadratic expression in compo- 
nents of A, where the particular indices do not matter. In the rest of this part, we will 
adopt this schematic notation for lower order terms (e.g. terms of the type and A^) 
where the particular indices do not matter. 

Finally, DAq and DAi will always be understood as □(Aq) and □(Aj), while (□A)^, 
refers to the tensorial wave equation. Also, AAq will always refer to A.{Ao). 

Remark 2.3. Since do and dj are not coordinate derivatives, note that the commuta- 
tors [dj,do\ and [dj,di\ do not vanish. In fact we have, schematically, 

[8i, 8j](f) = Ad(j) and [8j, 8o](j) = Adcj), (2.22) 

for any scalar function (f) on Ai. 

2.1.4. Main equations for (^40, A). Using the conventions above one can prove the 
following proposition. 

Proposition 2.4. Consider an orthonormal frame 6q, compatible with a maximal 
Ej foliation of the space-time A4 with connection coefficients defined by (2.2), their 
decomposition A — {Aq, A) relative to the same frame e^, and Coulomb- like condition on 
the frame, 

div A = A^. 
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In such a frame the Einstein-vacuum equations take the form, 

AAo = AdA + Ad{Ao) + A\ (2.23) 

DAi + di{doAo) = A^djAi + AWiAj + AodA + Ad{Ao) + A\ (2.24) 

Remark 2.5. It is extremely important to our strategy that we have reduced the co- 
variant wave equation (2.10) to the system of scalar equations (2.23) (2.24) (see remark 
1.15). 

2.2. Strategy of the proof of the main Theorem 

In this section, we discuss the strategy of the proof of the main theorem after reduction 
to smaU initial data, i.e. Theorem 1.14. 

2.2.1. The Uhlenbeck type lemma. In order to exhibit a frame ei, 62, 63 such that 
together with eo = T" we obtain a connection A satisfying our Coulomb type gauge on 
the shce E^, we will need the following result in the spirit of the Uhlenbeck lemma [51]. 

Lemma 2.6. Let (M, g) a 3 dimensional Riemannian asymptotically flat manifold. Let 

R denote its curvature tensor and r^,oi(M, 1) its volume radius on scales < 1. Let A a 
connection on M corresponding to an orthonormal frame. Assume the following bounds: 

\\Ml^{m) + II Vi||L2(M) + \\R\\L^iM) < S and r^oi{M, 1) > ^, 

where 5 > is a small enough constant. Assume also that A and VA belong to L'^{M). 
Then, there is another connection A on M satisfying he Coulomb like gauge condition, 
and such that 

WMWm + l|VA||i2(M) < s 

Furthermore, ifW'^A belongs to L^{M), then W^A belongs to L'^{M). 

The proof of Lemma 2.6 is a straightforward adaptation, in a simpler situation, of 
[51]. 

2.2.2. Classical local existence. We rely on the following standard well-posedness 
result for the Cauchy problem for the Einstein equations (1.1) in the maximal foliation. 

Theorem 2.7 (Well-posedness for the Einstein equation in the maximal foliation). 

Let {HQ,g,k) be asymptotically flat and satisfying the constraint equations (1.2), with 
Ric, VRic, k, Vk and V^A; in L^(So), and r^oi(So, 1) > 0. Then, there exists a unique 
asymptotically fiat solution (A^,g) to the Einstein vacuum equations (1.1) corresponding 
to this initial data set, together with a maximal foliation by space-like hypersurfaces 
defined as level hypersurfaces of a time function t. Furthermore, there exists a time 

T, = T;(||V«i?ic|U2(So),0 < I < 1, ||V(^U|U2(s„),0 < j < 2,r,oi(^o, 1)) > 



24 2. EINSTEIN VACUUM EQUATIONS AS A YANG-MILLS GAUGE THEORY 

such that the maximal foliation exists for on < t < with a corresponding control in 
L[J^ ]L2(Et) for Ric, VRic, k, Vk and V^k. 

Theorem 2.7 requires two more derivatives both for R and k with respect to the main 
Theorem 1.6. Its proof is standard and rehes solely on energy estimates (as opposed 
to Strichartz or bilinear estimates). We refer the reader to [9] chapter 10 for a related 
statement. 

Remark 2.8. In the proof of our main theorem the result above will be used only as 
an extension tool (see steps 1 and 3 below), only for very tiny values of the time interval. 

2.2.3. Weeikly regulcir null hypersur faces. We shall be working with null hyper 

surfaces in M. verifying a set of reasonable assumptions, described below. These as- 
sumptions will be easily verified by the level hyper surfaces l-Lu solutions u of the cikonal 
equation g'^'-'d^di, = discussed in section 2.7. The regularity of the eikonal equation is 
studied in detail in [44] (see also Chapter 4). 

Definition 2.9. Let % be a null hypersurface with future null normal L verifying 
g{L,T) — —1. Let also N — L — T. We denote by ^ the induced connection along the 
2-surfaces Tin Tit. We say that T-L is weakly regular provided that, 

||DL|U3(^) + ||DA^|U3(H) < 1, 

and the following Sobolev embedding holds for any scalar function f on %: 

wnwu) < wffWL^iu) + wmwL^iu) + wfWL^m- (2-26) 

2.2.4. Main bootstrap assumptions. Let M > 1 a large enough constant to be 
chosen later in terms only of universal constants. By choosing e > sufficiently small, we 
can also ensure Me is small enough. From now on, we assume the following bootstrap 
assumptions hold true on a fixed interval [0, T*], for some < T* < 1. Note that % 
denotes an arbitrary weakly regular null hypersurface with future normal L, normalized 
by the condition g{L,T) = —1. 

• Bootstrap curvature assumptions 

||R|U-L2(s,) < Me 

Also, 

\\R- L\\l2(^u) < Me, (2.28) 

where R-L denotes any component of R such that at least one index is contracted 
with L. 

• Bootstrap assumptions for the connection We also assume that there exist A = 
(^0,^) verifying our Coulomb type condition on [0,T*] , such that, 

||^IUrL2(so + ||9^||l-l2(eo + ||^||L?Lr(E,) < Me, (2.29) 



(2.25) 



(2.27) 
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and: 

||^o||L^L2(St) + ||^^o||L(°°L2(Et) + ||^o||L2L°°(St) + II || L^L3(Sj) 

+ ll^^^o||,.,3^,^^ < Me. (2.30) 

Remark 2.10. Together with the estimates in [44] (see section 4-4 that paper, and 
also section 4-2.3), the bootstrap assumption (2.27) yields: 

\\k\\LrL^{i:,) + ||VA;|Uc«^.(s,) < Me. (2.31) 

Furthermore, the bootstrap assumption (2.28) together with the estimates in [44] (see 
section 4-2 in that paper, and also section 4-2.2) yields: 

infr,„KSt,l) > \. (2.32) 

In addition we make the following bilinear estimates assumptions for A and R.: 
• Bilinear assumptions I. Assume, 

\\A^d,A\\vHM) < Mh\ (2.33) 

Also, let B = {—A)~^curl A (see (2.62) and the accompanying explanations). 
Then, we have: 

\\A^d,{dB)\\L.^M)<Mh^ (2.34) 

and: 

\\R..jod'B\\L2^M)<Mh^. (2.35) 

Finally, for any weakly regular null hypersurface T-L and any smooth scalar func- 
tion (f) on M., 



and 



\kj.d^<P\\LHM) < supWf <j)\\LHn), (2.36) 

H 



\\A^djcl>U.^M) < MhsupmUnu), (2.37) 

n 



where the supremum is taken over all null hypersurfaces T-L. 
• Bilinear assumptions II. We assume, 

||(-A)-i(g,,(^,^))|U.(^) < MV, (2.38) 

where the bilinear form Qij is given by Qij{(f), ip) — di(l)djip—dj(l)diip. Furthermore, 
we also have: 

\\{-A)-kd{A%A)U.^M)<Mh'. (2.39) 



26 



2. EINSTEIN VACUUM EQUATIONS AS A YANG-MILLS GAUGE THEORY 



• Non-sharp Strichartz assumption 





sharp. Nevertheless, this estimate will be sufficient for the proof as it will only be used to 
deal with lower order terms. 

Finally we also need a trilinear bootstrap assumption. For this we need to introduce 
the Bell Robinson tensor, 



Jm 

Let us conclude this section by remarking that the bootstrap assumptions are verified 
for a sufficiently small final value T*. 

Proposition 2.12. The above bootstrap assumptions are verified on < t < T* for 

a sufficiently small T* > 0. 

The only challenge in the proof of Proposition 2.12 is to show the existence of the 
desired connection A using in particular the Uhlenbeck type Lemma 2.6. All other esti- 
mates follow trivially from our initial bounds and the local existence theorem above, for 
sufficiently small T*. We refer to Proposition 4.6 in [27]. 

2.2.5. Proof of the bounded curvature conjecture. In the following two 
propositions, we state the improvement of our bootstrap assumptions. 

Proposition 2.13. Let us assume that all bootstrap assumptions of the previous sec- 
tion hold for < t < T* . If e > is sufficiently small, then the following improved 
estimates hold true on < t < T* : 



• Trilinear bootstrap assumption. We assume the following. 



(2.42) 




(2.43) 



ll^ll 



(2.44) 



(2.46) 



(2.45) 



+ ||^A||L?°L3(St) + \\ddAo 



0||L2L°°(Et) 



< e + M^si+MV, (2.47) 
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Proposition 2.14. Let us assume that all bootstrap assumptions of the previous sec- 
tion hold for < t < T* . If e > is sufficiently small, then the following improved 
estimates hold true on < t < T* : 

\\A^djA\\r^2^M) < (2.48) 

\\A^d,{dB)\\i^.^M)<M\\ (2.49) 

and 

\\R..j^d^B\\L2^M)<Mh\ (2.50) 
Also, for any scalar function on M., we have: 

\\k^.d^(t>\W^M) < MesMmWL^iu), (2.51) 

and 

WA^dML^M) < MesupmUHH), (2-52) 

H 

where the supremum is taken over all null hypersurfaces T-L. Finally, we have: 

||(^A)-i(Q,,(AA))|U.(^) <MV, (2.53) 

\\{-A)-kdA%A)\\^.^M) < MV. (2.54) 

PIIl?l^(e.) < Me, (2.55) 

ll'9i?|U?LT(E,) < Ms. (2.56) 

and 



M 



< (2.57) 



The proof of Proposition 2.13 is postponed to section 2.6, while the proof of Proposition 
2.14 is postponed to section 2.8. We also need a proposition on the propagation of higher 
regularity. 

Proposition 2.15. Let us assume that the estimates corresponding to all bootstrap 

assumptions of the previous section hold for < t < T* with a universal constant M. 
Then for any t G [0,T*) and for £ > sufficiently small, the following propagation of 

higher regularity holds: 

||DR||iooi2(St) < 2 (||i?ic||L2(So) + \\VRic\\L2(s„) + ||A-||l2(So) + II V/i;||l2(Eo) + || V^/c||l2(So)) . 

The proof of Proposition 2.15 follows along the same lines as the proof Proposition 
2.13 and Proposition 2.14, and we refer to [27] for its proof. Next, let us show how 
Propositions 2.12, 2.13, 2.14 and 2.15 imply our main theorem 1.14. We proceed, by the 
standard bootstrap method , along the following steps: 
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Step 1. We show that all bootstrap assumptions are verified for a sufficiently small final 
value T*. 

Step2. Assuming that all bootstrap assumptions hold for fixed values of < T* < 1 and 
M sufficiently large we show that, for £ > sufficiently small, we may improve 
on the constant M in our bootstrap assumptions. 
Step 3. Using the estimates derived in step 2 we can extend the time of existence T* to 
T* -\- 5 such that all the bootstrap assumptions remain true. 

Now, Step 1 follows from Proposition 2.12. Step 2 follows from Proposition 2.13 
and Proposition 2.14. In view of Step 2^ the estimates corresponding to all bootstrap 
assumptions of the previous section hold for < t < T* with a universal constant M. 
Thus the conclusion of Proposition 2.15 holds, and arguing as in the proof of Proposition 
2.12, we obtain Step 3. Thus, the bootstrap assumptions hold on < i < 1 for a universal 
constant M. In particular, this yields together with (2.31): 

llRlUj^LaCso < £ and ||A;||lool2(s,) < £ for all < i < 1. 

In view of (2.32), we also obtain the following control on the volume radius 

Furthermore, Proposition 2.15 yields the following propagation of higher regularity 

|a|<m |i|<m 

6i55[||V«Ric|U2(So) + ||V»VA;|U2(So) 

where Cm only depends on m. 

Remark 2.16. Note that Proposition 2.15 only yields the case m = 1 in (2.60). 
The fact that (2.60) also holds for higher derivatives m > 2 follows from the standard 
propagation of regularity for the classical local existence result of Theorem 2. 7 and the 
bound (2.60) with m — 1 coming from Proposition 2.15. 

Finally, (2.58), the control on the volume radius (2.59) and the propagation of higher 
regularity (2.60) yield the conclusion of Theorem 1.14. Together with the reduction to 
small initial data performed in section 1.2.3, this concludes the proof of the main Theorem 
1.6. 

The rest of the chapter deals with the proofs of propositions 2.13 and 2.14. The core 
of the proofs is to control A, the spatial part of the connection A. As explained in the 
introduction we need to project our equation for the spatial components A onto divergence 
free vectorfields. This is needed for two reasons, to eliminate the term c?i(c?oAo) on the left 
hand side of (2.24) and to obtain, on the right hand side, terms which exhibit the crucial 



(2.58) 
(2.59) 



(2.61) 
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null structure we need to implement our proof. Rather than work with the projection V, 
which is too complicated, we introduce instead the new variable, 

B = {-A)-^curl (A) (2.62) 

for which we derive a suitable wave equation. Since we have (see Lemma 2.20) A = 
curl (B) + l.o.t it suffices to obtain estimates for B which lead us to an improvement of 
the bootstrap assumption (2.29) on A. In section 2.4, we derive space-time estimates for 
OB and its derivatives. Proposition 2.13, which does not require a parametrix represen- 
tation, is proved in 2.6. Proposition 2.14 is proved in sections 2.8 and 2.9 based on the 
representation formula of theorem 2.30 derived in section 2.7. 

2.3. Simple consequences of the bootstrap assumptions 

2.3.1. Sobolev embeddings and elliptic estimates on E^. The bootstrap as- 
sumption (2.27) on R and the estimate for k (2.31) together with the estimates in [44] 
(see section 4.4 in that paper) yield the following lapse estimates: 

\\n - M\l'^{M) + ||Vn||Loo(_A4) + ||V^n||iooi2(St) + ||V^n||iooi3(s^) (2.63) 
+||V(9on)|U..3(,,) + ||V^n||^^^^3^^^^ + ||V^(9o(n))||^^^3^^^^ < Me, 

where V denotes the induced covariant derivative on Ej. 

Remark 2.17. Recall from (2.17) that (^o)oi = — n~^Vin. Thus, the estimates (2.63) 
forn could in principle be deduced from the bootstrap assumptions (2.30) for Aq. However, 
notice that Vn G L°°{A4) in view of (2.63), while Aq is only in L^L°°{T:t) according to 
(2.30). This improvement for the components (Ao)oi of Aq will turn out to be crucial (see 
remark 2.24)- Rs proof is given in section 4-4 of [44] (see also the discussion in section 
4.2.3). 

Next, we record the following Sobolev embeddings and elliptic estimates on derived 
under the assumptions (2.28) and (2.27) in [44] (see sections 3.5 and 4.2 in that paper). 

Lemma 2.18 (Calculus inequalities on S^). Assume that the assumptions (2.28) and 
(2.27) hold, and assume that the volume radius at scales < 1 on Sq is bounded from below 
by a universal constant. Then, the Sobolev embedding on E^ holds for any tensor F 

\\F\\lH^,) < ||VF|U2(s,), (2.64) 
Also, we define the operator (— A)~5 acting on tensors on Ej as: 

-1 /■+0O 

(-A)-^F=— ^ / T--^U{T)FdT, 

where F is the Gamma function, and where U{t)F is defined using the heat flow on E^; 

(dr - A)U(r)F = 0, U(0)F = F. 
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We have the following Bochner estimates: 

||V(-A)-5|U(i.(s,)) < 1 and ||V2(-A)-ir(i2(s,)) < 1, (2.65) 

where £(L^(St)) denotes the set of hounded linear operators on L^(Et). (2.65) together 
with the Sobolev embedding (2.64) yields: 

||(-A)-^F|U.(s,)<||F||^6^^^^. (2.66) 

2.3.2. Elliptic estimates for B. Here we record simple estimates for B, based the 
bootstrap assumptions (2.29) (2.30) for A and Aq and standard elliptic estimates such as 
the Bochner and Sobolev inequalities on E^, see (2.65) and (2.66). 

Proposition 2.19. Let Bi = (-A)~^(atr/ {A)i). Then, we have, for each component 
ofB: 

ll^^lUr^^CsO + ll^'^IUrL^Cs,) + ||9(9o5)|U,-l2(2,) < Me. (2.67) 

2.3.3. Decomposition for A. Recall that B = {—A)~^{curl (A)). We indicate 
below how to recover A from B: 

Lemma 2.20. We have the following estimate: 

A = curl {B) + E 

where E satisfies: 

l|9£;|U^c.^3(s,) + \\d'E\\^^^s^^^^ + ||i?|U2^o.(s^) < mV. 

Proof. We have, symbolically, 

A = {-A)-^curl {curl {A) + {-A)-\AdA + A^). 

from which, 

A = curl{B)-{-A)-^[A,curl]B + {-A)-\AdA + A^) 

The rest of the proof uses elliptic estimates on E^, the bootstrap assumptions for A and 
R, and the bootstrap assumption (2.40). We refer to [27] for the details. □ 

2.4. Estimates for DB 

We outline the proof of two important propositions concerning estimates for Ocurl A 
and OB, with B = A^^curl {A). The proofs makes use of the special structure of various 
bilinear expressions and thus is based not only on the bootstrap assumptions for tIq, ^, 
k and R but also some of our bilinear bootstrap assumptions. 

We record first the following straightforward commutation lemma, see [27]. 
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Lemma 2.21. Let (p a so{3, 1) scalar function on Ai. We have, schematically, 

dj{Dcl)) - = 2{A^)j dxd^cl) + do{Ao)dct> + A^d(t>. (2.68) 

We also have: 

[□, A]0 = -4r''V„V6(9o0) + 471-^571 V6(9o(9o0)) - 2VoA;«^V„V60 (2.69) 

F(2) = 9aAo + AaA + A^ 

where Va and denote induced covariant derivatives on Ej applied to the scalars 4>, docj) 
and do{do(p). 

The estimates for Ocurl A and OB are given by the following propositions. 
Proposition 2.22. We have 

3 

5^||(-A)-in(cwr/ A,)|U.(^,) < MV. (2.70) 

i=l 

Proposition 2.23 (Estimates for dB). The components Bi = {-Ay^{curl {A)i) 
verify the following estimate, 

3 

J2 iimu^M) + wdDBiU^^M)) < mv. (2.71) 

1=1 

We also have, 

3 

J2\\dodoBi\\L^^M) < Me. (2.72) 

i=l 

The proof of Proposition 2.22 and Proposition 2.23 are similar in spirit. We give below 
a short outhne of the proof of Proposition 2.23 which is slightly more difficult. 

Proof. In what follows we outline the main steps in the proof of space-time estimates 
(2.70), (2.71) for DB and d^B. We have: 

D{B,) = [D,{-A)-']{curl{A),) + {-A)-\D{curl{A),)) 

= -(-A)-^p, A]{-A)-\curl (A),) + {-A)-\n{curl {A)^)) 
= -(-A)-^p, A](5,) + {-A)-\n{curl{A),)). 

Thus, using the boundedness of 5(— A)^/^ and result of proposition 2.22, we obtain: 
\\dD{B,)\\^.^M) < ||(-A)-^[n,A](A)||L^(^) + MV, (2.73) 
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It remains to estimate ||(— A)~5[n, A](i?j)||j;,2(_yv(). We rely on the commutator formula 
(2.69) to write, 

[D,A]B = -4r^VaV5(5o5)+4n-V„nVb(ao(5o5))-2Vor''V„V„5 (2.74) 

F(2) = dd{Ao) + AdA + A\ 
with B any component (Bi), I = 1,2,3. It is easy to check that. 



We write. 



II (-A)-^ [□, A]{Bi) \\l.^m) <Ni + N2 + Ns + N, (2.75) 

N, = Wi-Ay^k^'WaVbidoiBimimM) 

7V2 = \\i-A)-'^n-'V,nV,{doido{Bi)))]\\L2^M) 

Ns = \\{-A)-^Vok'''WaVb{Bi)]\\L2^M) 

N, = \\F^'^\\LrL^^M\d\Bi)U.iM) + \\F^'^\\^^^i^^Jd{Bi)h^L.^^,) 
Using the estimates (2.67) we easily infer that 

< Me^. (2.76) 
To estimate A^i we proceed as follows, using the constraint equations (1.7) for k, 

k'^'VaVbidoiBi)) = Vjfc"^V6(9o(S0)] 

Together with the Bochner inequahty on (2.65) and the bilinear assumption (2.34) , 
we obtain: 

A^i < \\k^%{do{Bi))]U^M) < MV. (2.77) 

To estimate A'2 we write, 

n-^VbnVb{do{do{Bi))) = V'[n-Wbndo{do{Bi))] - {n-'An - n-^\Vn\^)do{do{Bi)). 

Together with the estimates (2.63) for the lapse n and the Sobolev embedding on 
(2.66), this yields: 

N2 < \\n-'Vbndo{do{Bi))\\L.^M) + \\{n-'An - n-^\Vn\^)do{do{Bi))\\^,j^^^^ 

< (llVnIUoo + ||n-^An - n-^\Vn\''\\LrLHj:t))\\do{do{Bi))\\L2iM) 

< Me\\do{do{Bi))\y^My (2-78) 

Remark 2.24. Note that there is no room in the estimate (2.78). Indeed the sharp 
estimate ||Vn||ioo(^) < Me given by (2.63) is crucial as emphasized in remark 2.17. 
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Finally, we consider N3. Recall from (1.10) that the second fundamental form satisfies 
the following equation: 

Vokab = - n-^VaVhn - kach ' = + l-O-t. (2.79) 

where E is the 2-tensor on defined as Eab = R-a06o- In view of the estimates (2.31) for 
k and (2.63) for n, 

II l.o.t.||L-L3(so < ||V^n||L-L3(2^) + ||/i:||iooi6(E,) < Me. 

Thus instead of estimating VoA;"''VaV;,i? in the definition of A^3 it suffices to estimate the 
term 

E-^VaVbB = Va(E-'VbB)-VaE'''VbB 

Using the maximal foliation assumption, the Bianchi identities and the symmetries of R, 
we can write, schematically, V°i?a6 = AR and therefore, together with the bootstrap 
assumptions (2.29) for A and (2.30) for Aq, and the bootstrap assumption (2.27) for R 
yields: 

ll^^^-^llLrLicEO ^ l|A|UooLa(E.)||R|UrL^{EO < MV. (2.80) 

We thus have, 

^okab = VaiE'^'VbB) + l.o.t. 
Now, in view of the bilinear assumption (2.35), 

\\{-Ay/'Va{E-'Vt,B)\\L2^M) < ||R0a06Vbi?|U2(^) < 

Hence, putting all the above together we infer that, 

N, < MV. 

Together with (2.75), (2.77), (2.78) and (2.76), wc derive, 

\\dnB\\L2^M) < M'^e^ + Me\\do{doB)\\L2(^M)- (2.81) 

To close the estimate for ||c} □-B||l2(x) it remains to control the right-hand side of (2.81). 
This is achieved relying in particular on the following formula 

do{doB) = -DB + AB + n" • VB. 

□ 

2.5. Energy estimate for the wave equation on a curved background 

Recall that cq — T, the future unit normal to the foliation. Let tt be the deformation 
tensor of eo, that is the symmetric 2-tensor on M. defined as: 

In view of the definition of the second fundamental form k and the lapse n, we have: 

TTab = -'2kab: -^aO = T^Oa = TToo = 0. (2.82) 
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In what follows "H denotes an arbitrary weakly regular null hypersurface^ with future 
normal L verifying g(L, T) = —1. We denote by y the induced connection on the 2- 
surfaces n E^. 



Lemma 2.25. Let F a scalar function on M., and let 0o o-nd 0i two scalar functions 
on Eq. Let (p the solution of the following wave equation on M.: 

U<t> = F, 

0|so = 00, ^o(0)|so = 01- 

Let Sq,Si denote the energy quantities, 

£o[0] ||90||L,-L2(St) + sup(||y0||i2(^) + \\L{(f))\\L2^n)) 



(2.83) 



H 



||9(a0)|UrL^(sO + \\do{docf>)\\LHM) + sup (11^(90) + ||L(90)|U2(^)) 



where the supremum is taken over all weakly regular null hypersurfaces % (satisfying 
assumptions (2.25) and (2.26) ). The following estimates hold true, provided that sM"^ is 
sufficiently small, 

^0 < \\VM\lHj^o) + \\MlHj^o) + \\F\\lHm), (2.84) 
^1 < I|V'0o||l=^(So) + I|V0i|U2(s„) + ||VF|U2(^). (2.85) 

Proof. We introduce the energy momentum tensor Qa,3 on M. given by: 

In view of the equation (2.83) satisfied by 0, we have: 

Now, consider the divergence of the 1-tensor Pa — Qa^eg = Qao, 



where vr is the deformation tensor of cq. Integrating over well-chosen regions oi M., we 
easily obtain: 



^0 < l|V0o||i2(So) + II0i|Il2(So) + 



FdoHM 



M 



M 



(2.86) 



< 



l|Wo||i2(So) + ||0l||i2(So) + ||i^|U2(.M)||9o0|U2(^) + 



M 



Next, we deal with the last term in the right-hand side of (2.86). In view of (2.82) 
and our maximal foliation assumption, we have: 



Qa^TT'^^dM 



-2 



M 



da(l)db(j)k^''dM + I n-'W''nda(pdo(t>dM 



M 



\.e. it satisfies assumptions (2.25) and (2.26) 
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Together with the bihnear bootstrap assumptions (2.36) and the estimates (2.63) for the 
lapse n, this yields: 



M 



< 



< 



\\ka-d''(l)\\L^M)\\d(l)\\L^M) + l|Vn||L^(^)||a0||^2(^) 



which together with (2.86) concludes the proof of the (2.84). Though more technical the 
proof of (2.85) follows the same ideas, and we refer to [27] for the details. □ 

2.6. Improvement of the bootstrap assumptions (psirt 1) 

In this section, we discuss the proof of Proposition 2.13. More precisely, we derive 
estimates for R, Aq and A which allow us to improve the basic bootstrap assumptions 
(2.27), (2.28), (2.29) and (2.30). 

2.6.1. Curvature estimates. We derive the curvature estimates using the Bell- 
Robinson tensor. 



Let 



(2.87) 



Then, we have: 

where tt is the deformation tensor of cq. We introduce the Riemannian metric, 

hap = 5'a/3 + 2(eo)a(eo)/3 (2.88) 
and use it to define the following space-time norm for tensors U : 

Given two space-time tensors U, V we denote hjU-Va given contraction between the 
two tensors and by \U -Vl the norm of the contraction according to the above definition. 

Let be a weakly regular null hypersurface with future normal L, g{L,T) = —1. 
Integrating (2.87) on a well-chosen, causal, space-time region, we have: 



/ |R|^+/ 



Jm Jm 



We need to estimate the term in the right-hand side of the previous inequality. Note that 
since ttoo = 0, ttoj = n^^V^n, and tt^- = kij, the bootstrap assumption (2.27) for R, and 
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the estimates (2.63) for n yield: 



\'R'-L\^ ^ + II V^-IU"" 11^-111°° z,2{St) + 



M 



< 



+ {Mef 



M 



The term in the right-hand side of the previous inequahty is dangerous. Schematically 
it has the form A;R^| . Typically this term is estimated by: 



Jm 



< 



L2L°°(St)||R'llLf'L2(St)' 



which requires a Strichartz estimate for k which is false even in fiat space. It is for this 
reason that we need the trilinear bootstrap assumption (2.43). Using it we derive, 



ist Jn 

which, for small £, improves the bootstrap assumptions (2.27) and (2.28). 



(2.89) 



2.6.2. Improvement of the bootstrap assumption for Aq. Recall (2.23) 

AAo = AdA + AdAo + A\ (2.90) 

Then, using (2.90), elliptic estimates on E^, and commuting (2.90) with dt in order to 
control dtAo, we are able to obtain the improved estimate (2.47) (see [27] for the details). 

2.6.3. Improvement of the bootstrap assumption for A. Using the estimates 

for HBi derived in Lemma 2.23, the estimates for B on the initial slice Eq, and the energy 
estimate (2.85) derived in Lemma 2.25, we have: 



\\d'B\\L^L-HE^)<s + Mh''. 

Using then (2.91) with Lemma 2.20, we obtain: 

||9A|Uoo^.(s^) < Wd'BW^^LHi:,) + ||9£;|UrL2(E.) < e + MV. 

which proves corresponding estimate in (2.46). 

To estimate OqA. we recall that, do{Aj) — dj{AQ) + Roj... Thus, we have: 



(2.91) 
(2.92) 



||<9o^||Lj°=L2(St) ^ II f^^O II Li°°L2 (Si 



|R| 



which together with the improved estimates for R and Aq yields: 



||9oA|Uoc^2(s,) <£ + (M£)2. 



(2.93) 
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2.7. Parametrix for the wave equation 

Let u± two families, indexed by a; e of scalar functions on the space-time M. 
satisfying the Eikonal equation for each u E S^. We also denote ^u±{t,x) = u±{t,x,u). 
We have the freedom of choosing '^u± on the initial slice Sq, and in order for the results 
in [43], [45] to apply, wc need to initialize '^u^ on Eg as in [42] (sec also Chapter 6). 

Let "H ^u± denote the corresponding null level hypersurfaces. Let '^L± its normal. 
'^L-t is null, and we fix it by imposing g( '^L±,T) = —1. Let the vectorfield tangent to 

'^A^-i- be defined such as to satisfy: 

We pick ( ^e±)A, A = 1,2 vectorfields in Ef such that together with '^N± we obtain 
an orthonormal basis of Ef. Finally, we denote by y/± derivatives in the directions 

Remark 2.26. Note thatU i^u± satisfy assumptions (2.25) and (2.26) from the results 
in [44] (see Theorem 2.15 and section 3.4 in that paper). 

For any pair of functions f± on M.^, we define the following scalar function on Ai: 

p poo p poo 

j/j[f+,f.]{t,x)= / e'^ ""+(*'^)/+(Au;)A2dAda; + / / e'^ ""-(*'^)/-(Au;)A'dAda;. 

J§2 Jo J§2 Jo 

We appeal to the following result from [43] [45] (see also Chapters 3 and 5): 

Theorem 2.27. Let 0o o,nd 0i two scalar functions on Eq. Then, there is a unique 
pair of functions (/+,/_) such that: 

V'[/+,/-]|so =00 anc? 9o(V'[/+,/-])|so = 0i- 
Furthermore, f± satisfy the following estimates: 

||A/+||l2(R3) + ||A/_||l2(]k3) < ||V0oi|L2(Eo) + ||0li|L2(Eo), 

and: 

||A^/+||l2(]R3) + ||A^/_||i2(]K3) < ||V^0o|U2(So) + ||V0i||l2(So). 

Finally, □■0[/+, /_] satisfies the following estimates: 

|PV'[/+,/-]|U2(^) < M£(||V(/)o|U2(s„) + ||0l|U2(Eo)), 

and: 

||9n^[/+,/_]||i2(^) < Me(||V20o||L2(Eo) + l|V0i||z.2(s„)). 

Remark 2.28. The existence of f± and the first two estimates of Theorem 2.27 are 
proved in [43] (see also Chapter 5), while the last two estimates in Theorem 2.21 are 
proved in [45] (see also Chapter 3). 
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We associate to any pair of functions (po, 0i on Sq the function 'ifom[<Po, 4>i] defined for 
{t,x) e as: 

where (/+, /_) is defined in view of Theorem 2.27 as the unique pair of functions associated 
to (00) 0i)- In particular, we obtain: 

||A/+||l2(r3) + ||A/_||l2(ib>3) < ||V0o||l2(So) + ||0i||l2{Eo), 

||AV+IU2(]R3) + ||A^/_||i2(]K3) < ||V^0o|U2(So) + ||V0i||l2(So), 

\\D^om[<i>o,<Pi]\\LHM) < M£(||V0o||l2(So) + WMl^L^o)), (2-94) 

and: 

< M£(||V20o|U2(s„) + ||V0i|U2(Eo))- (2.95) 

Next, let ^'^u± two families, indexed by a; e and s e R, of scalar functions on 
the space-time Ai satisfying the Eikonal equation for each G and s G M. We have 
the freedom of choosing ^'^u± on the slice Eg, and in order for the results in [43] [45] to 
apply, we need to initialize '^'*'m± on as in [42]. Note that the families '^u± correspond 
to '^'^u with the choice s — 0. For any pair of functions f± on M^, and for any s G R, we 
define the following scalar function on M.: 

f' f'OO f' /*CXD 

V^,[/+,/_](i,x,s)= / / e*^"''"+(*'"V+(Ac^)A2(iA(icu+ / / e'^^'"''-^''''^ f-{\u)\^d\du. 
Js2 Jo J §2 Jo 

We have the following straightforward corollary of Theorem 2.27: 

Corollary 2.29. Let s G R. Let (pQ and 4>i two scalar functions on E^. Then, there 
is a unique pair of functions (/+, /_) such that: 

■ips[f+J-]ks = 00 and do{^s[f+,f-])\j:s = 4>i- 
Furthermore, f± satisfy the following estimates: 

||A/+||l2(]r3) + ||A/_||l2(]r3) < ||V0o|U2(S,) + ||0l|U2(s^), 

and: 

||A^/+||l2(]r3) + ||A^/_||i2(IR3) < ||V^0o|U2(S3) + ||V0i||l2(s,). 

Finally, nipsif+yf-] satisfies the following estimates: 

||nV^,[/+,/_]|U2(^) < M£(||V0o||l2(s.) + ||0i||l2(s.)), 

and: 

\\dDi;s[f+, f-]\\mM) < M£(||V^0o||l2(e.) + ||V0i|U2(s,)). 
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Next, for any s G M, we associate to any function F on tiie function s)F 
defined for {t,x) & M. as: 

where (/+, /_) is defined in view of Corollary 2.29 as the unique pair of functions associated 
to the choice {(j)Q, (pi) — (0, F). In particular, we obtain: 

||A/+||l2(]k3) + ||A/_||i2(]K3) < ||F||i,2(s^), 

\\yu\\ 

L2(]R3) + II A /-||l2(]k3) < ||VF||£2(S^), 

||OT(t, s)F\\l2^m) < M£||F||i2(s,), (2.96) 

and: 

||an^(t, s)F\\l2^m) < Mell VF|U2(s,). (2.97) 
Now, we are in position to construct an exact parametrix for the wave equation (2.83): 

Theorem 2.30 (Representation formula). Let F a scalar function on M., and let 0o 
and 01 two scalar functions on Sq. Let (j) the solution of the wave equation (2.83) on 
Ai. Then, there is a sequence (p^^\ j > 0, of scalar functions approximations of (p and a 
sequence F^^\ j > 0, of scalar functions on M., with of the form: 

= *on^[0o, 0i] + f *(i, s)F^'\s, .)ds, F(o) = F 
Jo 

and for all j > 1: 

00) = / ^i){t,s)F^^\s,.)ds, 
Jo 

such that, 

+0O 

= ^0(^-), 

j=0 

and such that and satisfy the following estimates: 

W^^^Il^L^I:,) + lli^^^'^llL^(^) < (M£)^(||V0o|U^(Eo) + I|0i||l^(Eo) + lli^llL^(^)), 

and: 

im^'^h^LHT,) + ||9F(^-)|U2(^) < (M£)^-(||VVo||l2(Eo) + l|V0i|U2(So) + WdFh^t^M)), 
Proof. Let us define: 

F(°) = F and 0(°) = ^om[4>o, 0i] + / *(i, s)F^''\s, .)ds. 

Jo 

Then, we define iteratively for j > I: 

F(^) = -□0(^-^) + F^^-^) and 0^^) = / ^{t, s)F^^\s, .)ds. 

Jo 
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The proof follows from the estimates (2.94), (2.95), (2.96) and (2.97), together with the 
energy estimates for the wave equation of Lemma 2.25. We refer to [27] for the details. □ 

2.8. Improvement of the bootstrap assumptions (part 2) 

The goal of this section and next section is to prove Proposition 2.14. This requires 
in particular to write B using the representation formula of Theorem 2.30. In this sec- 
tion we derive the improved bilinear estimate (2.48), (2.49), (2.50), (2.51) and (2.52) of 
Proposition 2.14. We also derive the improved trilinear estimate (2.57). 

2.8.1. Improvement of the bilinear bootstrap assumptions I. In this section, 
we give the main ideas on the how we derive the improved bilinear estimate (2.48), (2.49), 
(2.50), (2.51) and (2.52) of Proposition 2.14. These bilinear estimates all involve the norm 
in L^{A4) of quantities of the type: 

C{U,dcf>), 

where C{U, dcj)) denotes a contraction with respect to one index between a tensor U and 
90, with being a scalar function which is solution to the wave equation (2.83) with F, 0o 
and 01 satisfying the estimate: 

||V'0o||l2(So) + ||V0i|U2(s„) + \\dF\\L2(M) < Me. 
In particular, we may use the parametrix constructed in Lemma 2.30 for 0: 

+0O 

= J]0(^-), 

3=0 

with: ^ 

0^°^ = *om[0O,0l]+ / ^it,s)F{s,.)ds, 

Jo 

and for all j > 1: 

00) = / ^(t,s)FO)(s,.)ds. 
Jo 

Thus, we need to estimate the norm in L'^{M) of contractions of quantities of the type: 

c(t/,a(^„^[0o,0i])) + V / c(t/,a(*(t,s)FW(s,.)))(is. 

i=o 

After using the definition of "^om and ^'(t, s), and the estimates for provided by 
Lemma 2.30, this reduces to estimating: 

/ / C(t/,a(e*^""+(*'^)))/+(Acu)A2dA(icu+ / / C(t/, a(e*^ ""-(*'")))/_(Acu)A'(iAdcu, 

J§2 Jo Jo 

where f± in view of Theorem 2.27 and the estimates for F, 0o and 0i satisfies: 

||AV±IU2(R3) <M£. 
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Since both half waves parametrices are estimated in the same way, the bihnear estimates 
(2.33), (2.34), (2.35), (2.36) and (2.37) all estimate the norm in L'^{M) of contractions of 
quantities of the type: 




C{U, d{e'^ ^"(*'")))/(Acc;)A2(iA(icc;, 



where / satisfies: 

||AVIU^(m3) < Me. (2.98) 

Furthermore we observe that 9j(e*'^"") = iXe^^'^'^dj{ ^u), and that the gradient of 
on Et is given by: V( '^it) = '^b~^ '^N, with '^6 = |V( '^u)\-^ is the null lapse, and 
'^N — 1^ is the unit normal to H '^uf^^t along S^. Thus, the bilinear estimates (2.33), 
(2.34), (2.35), (2.36) and (2.37) all reduce to L^(A4)-estimates of expressions of the form: 



/ e'^ ""b-^CiU, ''N)f{Xuj)X^dXduj, (2.99) 

_2 Jo 

where / satisfies (2.98). To estimate (t[U, f] we follow the strategy of [21]. 

mUJ]\\LHM) (2-100) 



< 



< 



< 



f{Xu)X^dX 



^/i^l II 



+0O 



(M)\\C{U, '^iV)|U-^L2(^„^) 

supll ''b-'h^^M)] fsup||C(^7, -iV)|Uc» 



ei^^<t^^)f{Xuj)X^dX 



dui 



m^u)] / \\^'f\^^)\\Lidu 



< fsupll ^b-^\\L^(M^ fsup||C(C/, '^iV)||L^„L2(«.j) IIAVII 



1,2 



Now, since has been initialized on Eq as in [42], and satisfies the Eikonal equation 
on the results in [44] (see Theorem 2.15 in that paper, and also (4.42)) under the 
assumption of Theorem 1.14 imply: 



supll "^6 ||lcx=(_a4) < 1. 
Together with the fact that / satisfies (2.98), and with (2.100), we finally obtain: 



§2 ^0 



^i\-u{t,x) -^jj-^dU, ''N)f{Xiu)X^dXdu; 



(2.101) 



< 



M£ sup||C(t/, ^iV)|U^^L2(H. 

Vu;eS2 
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It remains to estimate the right-hand side of (2.101) for the contractions appearing in the 
bihnear estimates (2.33), (2.34), (2.35), (2.36) and (2.37). Since all the estimates in the 
proof will be uniform in a;, we drop the index cu to ease the notations. 

Remark 2.31. In the proof of bilinear estimates (2.48), (2.49), (2.50), (2.51) and 
(2.52), the tensor U appearing in the expression C{U,N) is either R or derivatives of 
solutions (j) of a a scalar wave equation. In view of the bootstrap assumption (2.28) for the 
curvature fiux, as well as the first energy estimate for the wave equation in Lemma 2.25, 
we can control ||C(t/, A^)||ic»i2(^^^ as long as we can show that C{U, N) can be expressed' 
in terms of, R • L, ^0 and L(0). 

2.8.1.1. Proof of (2.48). Since A = curl {B)+E in view of Lemma 2.20 and bootstrap 
assumption (2.29), we have: 

\\A^d,{A)\\L^^M) < \\{curliB)yd,{A)\\LHM) + ||i?|U|L^(E,)||5^||LrL^(E,) (2.102) 
< ||(c«r/(5))^a,(yl)|U.(^) + MV, 

To estimate ||(c-ur/ {B)y dj{A)\\i^2(_j^) we write, {curl {B)ydj{A) =Ejmn dm{Bn)dj{A). We 
are now ready to apply the representation theorem 2.30 to B. Indeed, according to Lemma 
2.23, and proposition 2.19, we have 

DB = F, \\dF\\L2t^M) < (2.103) 

||9i^(0)|U^(Eo) + ||5'i^(0)|U^(Eo) + ||5(5oi?(0))|U^(Eo) < Me. 

Wc are thus in a position to apply the reduction discussed in the subsection above 
and reduce our desired bilinear estimate to an estimate for, 

CiU,N) = ejm.NmdjiA) 

Now, we decompose dj on the orthonormal frame N,fA,A= 1, 2 of E^, where we recall 
that /a. A — 1,2 denotes an orthonormal basis of Hu n E^. We have schematically: 

dj = NjN + f, (2.104) 

where y denotes derivatives which are tangent to n ^t- Thus, we have: 

e,m. N^dj{A) ^Ejm. N^NjdN{A) + f{A) = f{A), 

where we used the antisymmetry of e^^. in the last equality. Therefore, we obtain in this 
case: 

l|C(c/,7V)IUs=L^(«„)<||y(^)IUs.i2(«„). 



In other words, our goal is to check that the term C{U, N) does not involve the dangerous terms of 
the type a and Lcp, where L is the vectorfield defined as L = 2T — L, and a is the two tensor on D Hu 
defined by a^^ = Klalb- 
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It remains to estimate ||y(A)||iooj;^2(-^^). Since A = curl {B) + E we have, using Lemma 
2.20 again, followed by Proposition 2.23 and Lemma 2.25: 



iy(^)IU?>L2(w„) < \\f{dB)\\L^L2^n^) + \\f{E)\\L^LHn^) 



< \mdB)U^i^.^n.) + Me, 

< Ms. 

Therefore, 

\\A^dM)\\LHM) < \\{curl{B)ydM)\\L^iM) + Mh'<Mh', 

as desired. 

2.8.1.2. Proof of (2.49). The proof of (2.49) is similar to the one of (2.48) in view of 
Lemma 2.20. 

2.8.1.3. Proof of (2.50). Since B satisfies a wave equation in view of Lemma 2.23, the 
quantity C{U, N) is in this case^, 

CiU,N)^NjRoj..^RoN.. = ^Rll.. 

which together with the bootstrap assumption for the curvature flux (2.28) improves the 
bilinear estimate (2.35). 

2.8.1.4. Proof of (2.51). We have kj. = and A = curl (B) + E m view of Lemma 
2.20. Arguing as in (2.102), we reduce the proof to the estimate of: 

\\{curlBydj(f)\\L^M)- 

Then, the proof proceeds as the one of (2.48). 

2.8.1.5. Proof of (2.52). The proof of (2.52) proceeds as in (2.51). 

2.8.2. Improvement of the trilinear estimate. In this section, wc shall derive the 
improved trilinear estimate (2.57). To estimate the trilinear quantity \Jj^ Qij-ysk'^-' (^o^o] ■ 
we first write, according to Lemma 2.20, A = curl (B) + E by. Arguing as in (2.102), we 
reduce the proof of (2.57) to an estimate for: 

Q.j^s{curl {B))jeleo 

M 

Making use of the wave equation (2.103) for B we argue as in the beginning of section 
2.8.1 to reduce the proof to an estimate of the following: 

/ / / e'^^''^'^''^ '^b-^^jm- ''NmQj...)f{Xuj)X^dXdcudM 
Jm Js'^ Jo 

where / satisfies: 

||AV||l^(m3) < Me. 
^Use also L = T + N, L = T - N and the symmetries of R 
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Arguing exactly as in (2.100) (2.101), we can estimate the latter integral by the quantity 
suPa,es2||ejm. NmQj...\\L'^^^L^(ji^^)M£. In Other words, 



M 



< sup II e,-^. NrnQj...\y^ ii(^„jM£ + MV. 



(2.105) 



Next, we estimate the right-hand side of (2.105). Since all the estimates in the proof 
will be uniform in uj, we drop the index uu to ease the notations. The formula (2.42) for 
the Bell-Robinson tensor Q yields: 

Qj... — Hj^ ■ -H.x.. + dual 

1 1 
— — -RjL..R. L.. — -Rj L..R.L.. -|- Rj^..R.A-- + dual, 

where we used the frame L, L, fA,A — 1,2 in the last equality. Thus, we have schemati- 
cally: 

^jm- -^mQj--- R(R ■ L-\- Gjr'm- -^mRjr A") 

Decomposing ej with respect to the orthonormal frame A^, J'b, B = 1,2, we note that: 

On the other hand, decomposing Hba - further and using the symmetries of R, one easily 
checks that Hba-- must contain at least one L so that it is of the type R • L. Thus, we 
have schematically: 



NmQj... = 'RiR-L). 



(2.106) 



Thus, in view of (2.105), making use of the bootstrap assumptions (2.27) on R and (2.28) 
on the curvature flux, we deduce. 



M 



< {Mef + MsWRRlWlilh-h^) 

< {Mef + Me\\-R\\L2^M)\\^L\\L^L^{H^) 



In other words. 



Jm 



{Mef. 



(2.107) 



which yields the desired improvement of the trilinear estimate (2.43). 



2.9. Improvement of the bootstrap assumptions (part 3) 

In this section, we conclude the proof of Proposition 2.14. More precisely, we give the 
main ideas in the improvement of the bilinear bootstrap assumptions II. We start with a 
discussion of the sharp L'^{M.)- Strichartz estimate. 
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2.9.1. The sharp Strichartz L'^(A^) estimate. To a function / on and a family 
indexed by a; e of scalar functions on the space-time M. satisfying the Eikonal 

equation for each a; e S^, we associate a half wave parametrix: 

is2 Jo 

Let an integer p and a smooth cut-off function ip on (0, +oo) supported in a shell. We 
call a half wave parametrix localized at frequencies of size A ~ 2^ the following Fourier 
integral operator: 

/ / e*^""(*'"V(2"*'A)/(Aa;)A2rfArfa;. 

We have the following L'^{Ai) Strichartz estimates localized in frequency for a half wave 
parametrix which are proved in [46] (see also Chapter 7): 

PROPOSITION 2.32 (Corollary 2.8 in [46]). Let f a function on R^, let p e and let 
ip a smooth function on (0, +00) compactly supported in the shell 1/2 < A < 2. Let '^u a 
family indexed by cu G E)^ of scalar functions on the space-time Ai satisfying the Eikonal 
equation for each G and initialized on the initial slice Sq as in [42] . Let (f)p the scalar 
function on M. defined by the following oscillatory integral: 

p poo 

(f)p{t,x)= / e^^ ""(*'^V(2~^A)/(Aa;)A2dAda;. 

i§2 Jq 

Then, we have the following L'^[Ai) Strichartz estimates for (pp: 

UpWlHm) < 2i||^(2-^'A)/|U2(j,3), (2.108) 

WOMlHM) < 2^||V'(2-^'A)/|U2(M3), (2.109) 

\\d%\\LHM) < 2^||V'(2-^'A)/|U2(M3). (2.110) 
Note that this Strichartz estimate is sharp. 

2.9.2. Improvement of the non shcirp Strichcirtz estimates. Here, we derive 

the improved non sharp Strichartz estimates (2.55) and (2.56). In view of Lemma 2.20, 
(2.55) easily follows from (2.56), so we focus on the later improved estimate. 

Corollary 2.33. B satisfies the following Strichartz estimate: 

l|5i?l|LfLr(E0 < Me. 

Proof. Recall (2.103) which allows us to apply the representation formula of Theorem 
2.30 to 5. By a straightforward reduction the proof then reduces to the following non- 
sharp Strichartz estimate for a half wave parametrix: 



d ^ e'^ ""(*'^V(Aa;)A2ciAda;^ 



< ||AVIU^(R3). (2.111) 

L2L7(St) 

Then, the proof of Corollary 2.33 follows in particular from the sharp Strichartz estimate 
of Proposition 2.32. We refer to [27] for the details. □ 
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2.9.3. Improvement of the bilinear bootstrap assumptions II. In this section, 
we sketch the proofs of the improved bilinear estimates (2.53) and (2.54) of Proposition 
2.14. Based on the decomposition A = curl {B) + E oi Lemma 2.20 it is easy to show 
that that the proof of the bilinear estimates (2.38) and (2.39) reduces to: 

\\{-A)-HQ,j{dB,dB))\\r..^M) < MV. (2.112) 
Decomposing B according to Theorem 2.30, 

\\{-A)-HQ^J{^B,^B))U.^M) < E ll(-A)-^(g.,(0("^)>(")))||L^(A,). (2.113) 

m,n=0 

Thus it suffices to prove for all m, n > 0: 

||(-A)-i(g,,(0(-),</.(")))|U.(^) < (Me)-+i(Me)-+^ (2.114) 

The estimates in (2.114) are analogous for all m,n, so it suffices to prove (2.114) in the 
case {m,n) — (0,0). In view of the definition of 0^*^^ the estimates for B on the initial 
sfice Eg, estimate (2.71) for dUB, and the definition of and '^{t,s), (2.114) reduces 
to the following bilinear estimate for half wave parametrices: 

|(-A)-^(g,,-(0/,,0/j 

with, 

/»oe 

/ e'^""(*'"V(Acc;)A2rfA(icc;. 

We then decompose /i,/2 with respect to frequency and reduce the desired estimate to 
L'^{M) Strichartz estimate localized in frequency of Proposition 2.32, see details in [27]. 
This concludes the proof of Proposition 2.14. 



< 



1 1 ''^/l 1 1 L2 (RS ) 1 1 A/2 1 1 L2 (K3 ) 



(2.115) 




CHAPTER 3 



Control of the error term 

In this chapter, we consider the Fourier integral operator E given by (1.27) in which 
corresponds to the error term of a plane wave type parametrix. Recall that E is given by: 

Ef{t,x)= / e'^''^^'''''^^ngu{t,x,uj)f{Xuj)X^dXduj, {t,x) eM, 

where u{.,.,ijj) is a solution to the eikonal equation g^^daudi^u = on such that 
u{0,x,uj) X -uj when \x\ +oo on Eg (see section 3.1.1). The goal of this chapter is to 
outline the main ideas allowing us to obtain the control for the error term E in [45] . 

3.1. Geometric set-up and main results 

3.1.1. Geometry of the foliation of M. by u. Recall that is a solution to the 

eikonal equation g'^^daudjju = on depending on a extra parameter G S^. The level 

hypersufaces u{t,x,uj) = u of the optical function u are denoted by Let L' denote 
the space-time gradient of u, i.e.: 

L' = -gT^d^uda. (3.1) 
Using the fact that u satisfies the eikonal equation, we obtain: 

Di,L' = 0, (3.2) 

which implies that V is the geodesic null generator of l-Lu- 

We foliate the space-time M. by space-like hypersurfaccs St defined as level hypcrsur- 
faces of a time function t and we denote by T the unit, future oriented, normal to Ej. We 
have: 

T{u) = ±\Vu\ 

where |Vm|^ = Yl,i=i relative to an orthonormal frame on E^. Since the sign of 

T{u) is irrelevant, we choose by convention: 

T{u) = |Vm|. (3.3) 

We denote by Pt,„ the surfaces of intersection between E^ and 

Definition 3.1 {Canonical null pair). 

L^bL' + N, L^2T - L^T - N (3.4) 
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where L' is the space-time gradient ofu (3.1), b is the lapse of the null foliation (or shortly 
null lapse) 

b-^^-<L',T>^T{u), (3.5) 

and N is a unit normal, along E^, to the surfaces Pf^u- Since u satisfies the eikonal 
equation g'^^daudpu = on Ai, this yields L'{u) — and thus L{u) — 0. In view of the 
definition of L and (3.3), we obtain: 

(3-6) 



Definition 3.2. A null fram,e 61,62,63,64 at a point p G Pt,„ consists, in addition to 
the null pair 63 = L, 64 = L, of arbitrary orthonormal vectors 61, 62 tangent to Pt^u- 

Definition 3.3 {Second fundamental form). Let 61,62,63,64 be a null frame on Pt^u 
as above. The second fundamental form on Pt^u associated to our canonical null pair is 
given by 

Xab =< 0^64, cb > ■ 
We decompose x into its trace and traceless component. 

trx^Z Xab, Xab = Xab - ^trxZAB- 



Recall that trx satisfies a transport equation called the Raychaudhuri equation: 

1 
2 



L(trx) + ^(trx)' = -|x|' + --- (3.7) 



(see precise equation in (4.19)). 

We conclude this section with the identification of the symbol DgW of the error term. 
We have (see for example [44] for a proof): 

n^u = b~hrx. (3.8) 

Thus, we may rewrite the error term E as: 

Ef{t,x)^ / e'^''^^'''''^^b-\t,x,uj)trx{t,x,uj)f{Xuj)X^dXduj. (3.9) 

is2 Jo 

3.1.2. Some norms. We define some norms on H. For any 1 < p < +00 and for 

any tensor F on we have: 



I^IUp(h„) = ( / dt [ \F\Pdfit,u 

\J0 JPt,u 



where d/itu denotes the area element of Ptu- 
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Let x' a coordinate system on Po,«- By transporting this coordinate system along the 
null geodesies generated by L, we obtain a coordinate system (i, x') of "H. We define the 
following norms: 



x'ePo.u \Jo J 



x' 



sup 

0<t<l 



i2(Po,„) 



3.1.3. Estimates for the space-time foliation. In this section, we collect the 

estimates that are needed to follow the discussion of the control of the error term contained 
in this chapter. An outline of the proof of these estimates will be given in Chapter 4 (see 
[44] for the complete proof). 

We start with the regularity in (i, x) of the lapse h and the second fundamental for x- 
We need: 

lltrxlUoc(^) + II Vtrxllioc^a + 1|6- 1||l<-(^) + 1| V6||L^i2(-f^^) + 11x11^00^2 + \\fx\\L^L\Hu) ^ ^■ 

(3.10) 

Remark 3.4. In this section, all estimates hold for any a; e with the constant 

in the right-hand side being independent of u. Thus, one may take the supremum in u 
everywhere. To ease the notations, we do not explicitly write down this supremum. 

We also need an estimate for two derivatives of trx with respect to Vat 

\\VNPA'^Ni^x)\\mHu) ^ ^2^' + (3.11) 
where // in a function satisfying: 

||/^||l2(m) ^ 

Next, we consider the regularity with respect to ui. We have: 

Wh'-iM) < e, (3.12) 

\N{t,x,uj) - N{t,x,uj')\ ~ \uj-uj'\, y{t,x) e M,uj,uj' e §^ (3.13) 

and 

\\d^N\\Lo.iM) < 1. (3.14) 
Furthermore, we have the following decomposition for x- 

X = Xi + X2, (3.15) 

where Xi and X2 are two symmetric traceless Pt^^-tangent 2-tensors satisfying in particular, 
for any 2 < p < +oo: 

IIxiIIl^l- + II^^X2|lL«-m„) ^ ^- (3-16) 



Remark 3.5. The point of decomposition (3.15) is that Xi has a better regularity with 
respect to {t,x) than x? while X2 has a better regularity with respect to cu than x (see 
explanation in section 
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Finally, we need to compare quantities evaluated at two angles u and u in satisfying 
|a; — z/| < 2~2. We have the following decomposition for N{t,x,uj) — N{t,x, v): 

2i{N{t,x,u) -N{t,x,p)) = Fi{t,x,i^) + Fi{t,x,uj,p) (3.17) 
where the tensor Fl does not depend on cu and satisfies: 

l|i^lU^<i, 

and the tensor F| satisfies: 

Here L'^L'^{Tiu) is defined with respect to u — u{t,x,u). We also have the following 
decomposition for tr^: 

irx{t,x,uj) = f({t,x,u) + f^{t,x,uj,u) (3.18) 
where the scalar ff does not depend on ou and satisfies: 

ll/i'lk-<^, 

and where the scalar /2 satisfies: 

11/2 < £2-2. 

3.1.4. Main result. The main result of this chapter is the following. 

Theorem 3.6. Let u be a function on x satisfying suitable assumptions (we 
refer to [45] for the complete set of assumptions, and to section 3.1.3 for some typical 
assumptions). Let E the Fourier integral operator with phase u{t,x,uj) and symbol Hgu: 

Ef{t,x)= / / e'^'''-''''''^^b-\t,x,uj)trx{t,x,uj)f{Xuj)X^dXdw. (3.19) 

Jo 

Then, E satisfies the estimate: 

\\Ef\\mM)<e\\Xf\\L^^^sy (3.20) 

3.1.5. Geometric Littlewood-Paley projections on the 2-surfaces Pt,u- Through- 
out the paper, we will use the geometric Littlewood-Paley projections on 2-surfaces {Pt,u 
in our case) constructed in [24]. In that paper, the following properties are proved 

Theorem 3.7. The LP-projections Pj verify the following properties: 
i) U'-boundedness For any 1 < p < 00, and any interval I (Z'Z, 

\\PiF\\LP{Pt,u) ^ \\F\\LP{Pt,u) (3-21) 

a) Bessel inequality 

ll^i-^llL2(Pt,„) ^ ll-^llz,2(pj_^) 

j 
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in) Finite band property For any 1 < p < oo. 

\\^Pj^\\LPiPt,u) ^ '^'^^\\^\\LP{Pt,u) 
\\PjF\\LP{Pt,u) ^ '^~'^^\\^\\LP{Pt,u)- 

In addition, the estimates 

hold together with the dual estimate 

\\PjfF\\LHP,,.)<'^iF\\LHP.,.) 
iv) Weak Bernstein inequality For any 2 < p < oo 

\\P<oF\\LP{Pt,u) ^ l|-^IU2(Pt,„) 

together with the dual estimates 



(3.22) 



(3.23) 



< IIFI 



3.2. Control of the error term 

3.2.1. The basic computation. We start the proof of Theorem 3.6 with the fol- 
lowing instructive computation: 



WEfhHM) < 



S2 



b{t,x,uj)-hrx{t,x,uj) / e*^"/(Acj)A2rfA 



+00 



du) 



< / \\b{t,x,uj) hTx{t,x,u) 



e'^''fiXi^)X'dX 



duj 



< £||aVI|l2(m3), 



(3.24) 

where we have used Plancherel with respect to A, Cauchy-Schwarz with respect to uj, the 
estimates (3.10) for b and tr^. (3.24) misses the conclusion (3.20) of Theorem 3.6 by a 
power of A. Now, assume for a moment that we may replace a power of A by a derivative 
on b{t,x,u)~HTx(t,x,(x!). Then, the same computation yields: 



§2 Jo 



V(6(i, x, uj)~'trx{t, X, uj))e'^''f{Xuj)XdXduj 



< 



\\V{b{t, X, uj) hrx{t, X, uj)) ||l-l2(«„) 



r+oo 

/ e'^''f{Xuj)X^dX 
Jo 



duj 



(3.25) 



L?. 



< £||A/||l2(]r3), 
where we used the fact that 

\\W{b{t,X,Uj)-hTx{t,X,Uj))\\L^L^Hu) ^ £ 



(3.26) 
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in view of (3.10). Note that the estimate provided by (3.25) is consistent with the control 
of the error term (3.20). This suggests a strategy which consists in making integrations 
by parts to trade powers of A against derivatives of the symbol b{t,x,uj)~Hrx{t,x,uj). 

3.2.2. Structure of the proof of Theorem 3.6. The proof of Theorem 3.6 pro- 
ceeds in three steps. We first localize in frequencies of size A ~ 2-^. We then locahze the 
angle u in patches on the sphere of diameter 2"-'/^. Finally, we estimate the diagonal 
terms. 

3.2.2.1. Step 1: decomposition in frequency. For the first step, we introduce ip and if) 
two smooth compactly supported functions on R such that: 

99(A) + ^ V(2"^'A) = 1 for all A G M. (3.27) 

We use (3.27) to decompose Ef as follows: 

Ef{t,x)= J2Ejf{t,x), (3.28) 

i>-i 

where for j > 0: 

p r+oo 

Ejf(t,x)^ / e'^'"b(t,x,oj)-hrx(t,x,uj)'ijj(2-^X)f(Xuj)X'^dXduj, (3.29) 

Jo 

and 

E_if{t,x)^ / e'^''b{t,x,uj)-hTx{t,x,uj)ip{X)f{Xuj)X'^dXdw. (3.30) 

Jo 

This decomposition is classical and is known as the first dyadic decomposition (see [39]). 
The goal of this first step is to prove the following proposition: 

Proposition 3.8. The decomposition (3.28) satisfies an almost orthogonality prop- 
erty, from which it follows that: 

< J2\\^^f\\lHM)+^'\\f\\hiR^y (3-31) 
j>-i 

A discussion of the proof of Proposition 3.8 is postponed to section 3.3. 

3.2.2.2. Step 2: decomposition in angle. Proposition 3.8 enables us to estimate ||£'j/||i,2(^) 
instead of ||£^/||i2(_A4). The analog of computation (3.24) for ||i?j/||i2(_A4) yields: 

WEjfh^M) < £||AV'(2-^A)/|U2(,) < £2^||V'(2-^A)/|U2(n,3), (3.32) 

which misses the wanted estimate by a power of 2-'. We thus need to perform a second 
dyadic decomposition (see [39]). We introduce a smooth partition of unity on the sphere 



J2. 



Tj^icj) = 1 for all CO e (3.33) 
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where F is a lattice on of size 2~ i , where the support of rjj is a patch on of diameter 
~ 2"-'/^. We use (3.33) to decompose Ejf as follows: 

Ejfit,x) = J2E^fit,^), (3.34) 

where: 



We also define 



§2 Jo 



7-1 = lk(A)/|U2(R3), = \\i){2 ■'A)/||l2(r3), j > 0, 
7j = ||V'(2-^A)77j(a;)/|U2(M3), j > 0, e T, 



(3.36) 



which satisfy: 



Inm = E 7? = E E(Tr)^- (3-37) 

i>-i i>-i !^er 
The goal of this second step is to prove the following proposition: 

Proposition 3.9. The decomposition (3.34) satisfies an almost orthogonality prop- 
erty, from which it follows that 

WmlHM) ^ T.\\^jf\\lHM) + ^'ll (3-38) 

A discussion of the proof of Proposition 3.9 is postponed to section 3.5. 

3.2.2.3. Step 3: control of the diagonal term. Proposition 3.9 allows us to estimate 
||-E'J/IU2(A^) instead of \\Ejf\\L2(^j^y The analog of computation (3.24) for ||-Ej'/||L2(A^) 
yields: 

mf\WiM) (3.39) 



< / \\h{t,X,Uj) hYx{t,X,Uj)\\L^L-^(Hu) 



+00 



e^^X2-^A)77,"(a;)/(Aa;)A"dA 



dui 



< 2%^vol(supp(r;J))||AV^(2-^A)r;J(^)/|U2(K3) 



where the term Y'vol(supp(r7j)) comes from the fact that we apply Cauchy-Schwarz in oj. 

Note that we have used in (3.39) the fact that the support of r/J is 2 dimensional and has 
diameter 2~^/^ so that: 

^vol(supp(77j)) < 2-^/1 (3.40) 

Now, (3.39) still misses the wanted estimate by a power of 2^^^. Nevertheless, using more 
refined techniques, we are able to estimate the diagonal term: 
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Proposition 3.10. The diagonal term E^f satisfies the following estimate: 

WE^fh^M) < ei^. (3.41) 
A discussion of the proof of Proposition 3.10 is postponed to section 3.4. 

Remark 3.11. Note that Proposition 3.9 together with Proposition 3.10 yields the 
estimate: 

\\Ejf\\L-iM) < e^j. (3.42) 

Now, since the proof of Proposition 3.9 and the proof of Proposition 3.10 do not depend on 
the proof of Proposition 3. 8, we are allowed to use the conclusion of Proposition 3. 9 and 
Proposition 3.10 in the proof of Proposition 3.8. In particular, the estimate (3.42) will he 
used for the proof of Proposition 3.8. In the same spirit, since the proof of Proposition 
3.10 does not depend on the proof of Proposition 3.9, we are allowed to use the conclusion 
of Proposition 3.10 in the proof of Proposition 3.9. 

Convention. In the rest of this chapter, we will use several integration by parts. In 
turn, these integration by parts will each generate a large number of terms. For the sake 
of simplicity, we will only discuss few typical terms. We will constantly use the notation 

" H " in various identities and estimates in order to refer to the additional terms. That 

is not to say that these additional terms are lower order or estimated in the same way, but 
simply that the typical terms that we exhibit allow for a simple exposition of the main 
ideas of the proof. We refer the reader to [45] for a complete proof which contains the 
control of the typical terms discussed here as well as the numerous additional terms. 

3.2.2.4. Proof of Theorem 3.6. Proposition 3.8, 3.9 and 3.10 immediately yield the 
proof of Theorem 3.6. Indeed, (3.31), (3.37), (3.38) and (3.41) imply: 

WEfWh^M) ^ Y.\\^^f\\l^iM) + ^'\\f\\hm 

j>-l y&T j>-l (3.43) 

^-^EE(7jr+-'E7|+-^ii/iiiw 

j>-i uer j>-i 

^ ^ II/IIl2(R3), 

which is the conclusion of Theorem 3.6. 

The rest of this chapter is dedicated to a discussion of the proof of Propositions 3.8, 
3.9 and 3.10. The details of the proofs being very involved, we only give a very sketchy 
summary of the main ideas. We refer the reader to [45] for the details. 



3.3. ALMOST ORTHOGONALITY IN FREQUENCY 

3.3. Almost orthogonality in frequency 

We have to prove (3.31): 

This will result from the following inequality using Shur's Lemma: 



Ejf{t,x)Ekf{t,x)dM 



M 



<e^2 ^'^"'^j^k for \j-k\>2. 



In turn, (3.45) will follow from integrations by parts in u. 



55 



(3.44) 



(3.45) 



3.3.1. A first integration by parts. From now on, we focus on proving (3.45). We 
may assume j >k -\-2>. We have: 



Ejf{t,x)Ekf{t,x)dM 



(3.46) 



M 



+ 00 



+ CXD 



e""" "'^ b{t,x,uj) ^trx{t,x,uj)b{x,uj')~Hrx{t:X,uj')dM 
xi^{2-^\)f{\uj)\^il:{2-^\')f{\'uj'){>!fd\dujd>!duj'. 

We consider the coordinate system {t,u,x') on J\4, and we would like to integrate by 
parts with respect to du in (3.46). Since Vu = b~^N and Vu' — b'~^N', we have: 



AXu—iX'u' 



f) f iXu—iX'u'\ 



X-X'l,g{N,N- 



(3.47) 



where we use the notation u for u{t, x,u), b for b{t, x,u), N for N{t. x. uj). u' for u{t, x, cj'), 
b' for 6(t, ,T, uj') and N' for N{x, oj'). We will also use the notation trx for trx(t, x, uj) and 
trx' for trx(t, x, a;'). Using (3.47), we obtain: 



JM 



iXu—iX'u' 7 7/ 



bb'dM 



M 



X-X'l,g(N,N'y 



M X-X'^g{N,N') 



(3.48) 



where the additional terms in (3.48) arise when 9„ falls on the volume element of M. or 
on the denominator in the right-hand side of (3.47). Note that: 



Y 6 
Jb' 



< 



J 



<- + 0{e)<l, 
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where we used the estimate (3.10) satisfied by b and b' and the fact that j > k + 3 so that 
A' < A/2. Thus, we may expand the fraction in (3.48): 

p 



1 



1 X'^g{N,N') 



A 



A 



(3.49) 



Remark 3.12. The expansion (3.49) generates quantities of the type 

+00 



e^^X2-^A)/(Aa;)(2-^A)^'A"dA. 
where p e Z. For simplicity, we omit the index p and denote them by 



+00 



e^^"z^(2-^A)/(Aa;)A'dA. 



since they are essentially equivalent. Note that Plancherel yields: 

< ||^(2-^A)/(Aa;)A|U.(M3) < 2^7,. 
Also, using Cauchy-Schwarz in X, we have 

WFjWlIl^ < 2i||V'(2-^A)/(Aa;)A|U.(M3) < 2^1r 
(3.46), (3.48) and (3.49) imply: 



E,f{t,x)Ekf{t,x)dM 



(3.50) 

(3.51) 
(3.52) 

(3.53) 



IM 

+2- 



trxNFj{u)duj 



V{b'~\Tx')Fk{u')duj']dM 



IM \J&' 

where we only kept the first term in the expansion (3.49) in order to simplify the exposi- 
tion^. 

Remark 3.13. The second term in the right-hand side of (3.53) is easier because the 
derivative falls on the low frequency term. This is why we estimate this term directly while 
the other term requires a more elaborate treatment which is explained in section 3.3.2. 

We estimate the second term in the right-hand side of (3.53). We have: 

iYXNFj{u)du 



< 2- 



trxNFj{u)duj 

S2 



S2 



V{b'-hix')Fk{u')duj' )dM 
V{b'~hTx')Fk{u')du}' 

"L-^iM) 



(3.54) 



-"^note that in the last term in the right-hand side of (3.53), we wrote Vjv(^' ^trx') as N-'V{b' ^trx') 
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We have the foUowing analog of (3.42): 



tYxNFj{u)duj 



(3.55) 



Indeed, one can show that the symbol tr^; satisfies regularity assumptions which are at 
least as good as b^^tix (see [44], and also section 6.1.5), so that the proof of (3.42) may 
adapted in a straightforward manner to obtain (3.55). 

Next, we consider the second term in the right-hand side of (3.54). Then proceeding 
as in the basic computation (3.24), and using the estimate (3.26), we obtain 



V{b'-\rx')Fk{u')duj' 



S2 



< 



£||AV^(2-^A)/|U2(K3) < £2S- (3.56) 



Together with (3.54) and (3.55), we finally obtain: 



2-^1 (I tTxNFj{u)du] ■ ( I V {b'- 'trx') Fk{u')du' ]dM 



< 2-^-+S7., (3.57) 



which is consistent with (3.45). 



Remark 3.14. Estimating the first term in the right-hand side of (3.53) in the same 
way would only yield: 



{ / VNtrxFj{u)dw 

M 



which is not sufficient to obtain (3.45). 



b'-hrx'Fk{u')duj']dM 



(3.58) 



3.3.2. A more precise estimate. In this section, we estimate the first term the 
right-hand side of (3.53). Using the geometric Littlewood-Paley projections on the 2- 
surfaces Pt^u, we decompose Vjvtrx as: 

Vjvtrx = -P<i±fc(Vivtrx) + P^i±k{V N^rx)- 

— 2 ^2 

In turn, this yields a decomposition for the first term in the right-hand side of (3.53): 



2-M / VNtTxF,{u)du) / b'-\Tx'Fk{u')duj' ]dM^Ai + A2, 

IM \JS2 / \J§2 



(3.59) 



where: 



Ai = 2" 



P^i±k{VNtTx)Fj{u)du Ekf(t,x)dM, 



A2^2-^ [ ([ P^i+k{WNtrx)Fj{u)duj] Ekf{t,x)dM. 



(3.60) 
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We first estimate the easier term Ai. The definition of Pi imphes Pi = 2~^'^Pi, and 
thus 

P>i±i(Vivtrx) = Yl ^KVivtrx) 

= 2-2'APKVivtrx), 

2 

which yields the following decomposition for Ai. 

Ai^ Y ^1.' (3.61) 



i> 



j+k 



where Ai^i is given by: 

A,,i = 2-^-^' [ ([ ^i{VNtrx)Fj{u)dw] Ekf{t,x)dM. 
Integrating by parts ^ on Pt^u and using the fact that yFj_i{u) — 0, we obtain: 

Ai,i = -2-^-2' I I \ I fPii^Ntrx)fiEkfit, x)b)dfit,u I Fj{u)dudtdiu + ■■■ , 

J§2 .Jt,U \JPt,u J 

where the additional term corresponds to the case where the derivative falls on the volume 
element of Ai. Next, we apply Cauchy-Schwartz to the integral on A4 and obtain: 

\A^,l\ < 2-^-2' / \\fPi{Vr,trx)F, 

< 2-^-2' / ||yP,(Vivtrx)||Ls=L2(H„)||i^,(«)||Lsliy^fe||L2(M)C^^ 

< 2-^-' / \\VNtrx\\L^L^(Hj\Fj{u)\\Ll\\fEk\\L2^M-)dw 

< £2-^'-' / \\Fj{u)\\Li\\VEkh2^M)duj, 

where we used the finite band property for Pi and the estimates (3.10) for trx. In view 
of (3.62), we also need to estimate || V£^a;||l2(a^). We have: 

p r+oo 

VEkf{t,x)^ I I e'^''V(b-hrx)ij{2-''X)f{Xoo)X'^dXduj 

„+oo (3.63) 



'§2 Jo 



i> r+oo 

+i2^ / / e'^''h-hYxVuil){2-^X){2-''X)f{Xu)X^dXdu. 

is2 Jq 

Using the basic computation (3.24) for the first term together with the estimate (3.26), and 
(3.42) for the second term together with the fact that trx-^ satisfies the same regularity 
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(3.64) 



(3.62), (3.51), and (3.64) yield: 



Together with (3.61), this yields: 



(3.65) 



j + k 



which is consistent with (3.45). 



3.3.3. A second integration by parts in u. To estimate A2, we perform a second 
integration by parts relying again on (3.47). This leads to: 



12 - ^ - / / VNP<i±k(^NtTx)Fj(u)EJ(t,x)dM + 

Im Js2 - 2 



(3.66) 



where we only keep the worst term, which is the one containing two derivatives of tr^. It 
is at this stage that we need the estimate (3.11) for VnPi(^n^^X) which we recall now. 
We have: 

||V;vPKVivtrx)IL2(^„) < 62' + 2^/x(ii), (3.67) 
where in a function satisfying: 

||/^||l2(m) < e. 
In view of the estimate (3.67), wc have: 

VivP<i±^(VArtrx) ^^^^ ^ < J2 \\^NPii^Ntrx)\\mH^) 

L \Hu) 



2 



< 5^(2'£ + 2i/.(x.)) 



K 



j+k 



^ J + fc j + k 

< £2^+2^M«). 



In view of (3.66), this yields after applying Cauchy- Schwartz: 

IA2I < 2-'^\Ekf\\LHM) [ 

< 2-2^' 



duj + 



duj + 



< 2-2^' 



VAr-P<i±fc(VArtrx) 
ejk I {82"^ +e2'^fi{u))Fj{u) 
sjk(^e2'^ Jj\Fj{u)\\L2duj + 2^^ Jj\ii\\L2^^)\\Fj{u)\\L^du?j +■■■ 



(3.68) 



j—k 



< 2 4 £^7^7^. + 
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where we used (3.42) for E^f, Cauchy-Schwarz in cu, and the estimates (3.51) and (3.52) 
for Fjiu). 

3.3.4. End of the proof of Proposition 3.8. In view of (3.53), (3.57), (3.59), 
(3.65) and (3.68), we obtain: 

< £22-^7-7, for \j -k\>2. (3.69) 



Ejf{t,x)Ekf{t,x)dM 

M 

Finally, (3.69) together with Shur's Lemma yields: 

WEfWh^M) ^ E WEjfWlHM) + ^l/lli^(R3). (3.70) 
j>-i 

This concludes the proof of Proposition 3.8. 

3.4. Control of the diagonal term 

Since the orthogonality argument in angle is the core of this chapter, we choose to 
deal first with the control of the diagonal term in this section. We will then proceed with 
the orthogonality argument in angle in the rest of the chapter. 

In order to control the diagonal term, we have to prove (3.41): 

\\E^f\\LHM)<e^^. (3.71) 

Recall that Ej is given by: 

E^f{t,x)^ I b-\t,x,uj)tTx{t,x,uj)Fj{uX{u)dw, (3.72) 
where Fj{u) is defined by: 

r+oo 

Fj{u) = / e*^X2-^A)/(Aa;)A^dA. (3.73) 
Jo 

The proof of the estimate (3.71) will proceed in four steps: 
Step 1. We first consider a decomposition roughly of the type: 



E^fit, x) = b-\t, X, iy)trxit, x, v) [J^^ F, 



so that we have to prove estimate (3.71) with b ^trx replaced by 1. 
Step 2. That estimate is obtained by considering the transport equation along L^,: 



-'1/ 



Fj{u)rij{uj)dw 

Step 3. A certain term in the transport equation of Step 2 needs to be estimated using 

an energy estimate for the wave equation. 
Step 4. We conclude the proof using the estimates obtained in Step 2 and Step 3. 
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3.4.1. Step 1: freezing the u dependance in b ^trx- In view of the estimate 
(3.12) for dojb, the estimate (3.10) for 6, and the decomposition (3.18) for trx, we have: 

6"^ {t,x,u)trxit,x,u) = fl {t, x,iy) + f^{t,x, uj) , (3.74) 

where fl only depends on (t, a;, v) and satisfies: 

\\fi\\L^<s. (3.75) 



and where fl satisfies: 



ll/2lU-L2(^^) <2 26, 



(3.76) 



with u — u{.,uj). (3.74) yields the following decomposition for the diagonal term: 
E;;f{t,x)^fi{t,x,iy)[ Fj{u)rj^{u;)dw+ [ Fj{u)fi{t,x,u;,u)ri^{u;)duj, 

^§2 ^S2 

which implies: 

\l^{M) 



(3.77) 



< 



< 



\\fi\\L^{M) 



Fj{u)r]j{u)du 



+ 



§2 



Fj{u)r]':^{uj)duj 



L2(A1) 



\Fj{u)\\Li\\fi\\L^LHnu)Vj{<^)d<^ 



where we used in the last inequality the estimates (3.75) and (3.76), Cauchy-Schwarz in 
u!, the size of the patch, and the estimate (3.51) for Fj{u). 

Remark 3.15. The point of the decomposition (3.74) is to allow us to replace in the 
diagonal term (3.72) the symbol b~^trx with 1. An obvious way to achieve this is to write 
the following decomposition: 

b~^trx{t, X, to) = b~^trx{t, x, v) + (b~^trx{t, x, cu) — b~^trx{t, x, u)). (3.78) 

The first term clearly satisfies (3.75) in view of the estimate (3.10) for b and trx- On 
the other hand, we obtain in [44] (see also (4.49)^ the estimate d^jtrx G L^LP'il-Lu) which 
together with the estimate (3.12) for d^b yields: 



\mb-Hrx)\\ 



(3.79) 



Now, we have: 



b ^trx{t,x,oj) -b ^trx{t,x,u) ^ {oj -u) / doj{b ^trx){t,x,oj^)da 

Jo 

which together with (3.79) is not enough to conclude since L^LF'il-Lu) CLnd L'^^L?' ((Hua) 
are not comparable. We refer the reader to [44] where the decomposition (3.74) as well 
as several others are proved (see also the discussion in section 4-5) ■ 

The following proposition allows us to estimate the right-hand side of (3.77). 
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Proposition 3.16. We have the following bound: 



Fj{u)rfAijj)dw 



~ h 



(3.80) 



Remark 3.17. In order to control the diagonal term, it suffices to have a hound of 
the L'^{M) norm for the left-hand side of (3.80). The improvement to a hound for the 



Hi; tXi^ t 



Lf norm will he crucial when proving the almost orthogonality in angle. 



Assuming the result of the proposition, estimates (3.80) and (3.77) yield: 
which together with (3.72) and (3.77) implies: 

which is the wanted estimate (3.71). This concludes the proof of Proposition 3.10. 

3.4.2. Step 2: A transport equation in the Ly direction. We still need to prove 
Proposition 3.16. Note that it suffices to show: 



Fj{u)rfAuj)du: 



< 7^. 



(3.81) 



Now, since the space-time gradient of u is given by 6 ^L, we have: 

F^{u)ri'<{uj)du^ b-^g(L{t,x,Lj),L{t,x,u)^Fj{u)ri'^{Lj)dcj, 
where Fj has been defined in (3.50). In view of (3.82), we have: 

Fj{u)r)^{uj)du?j (3.82) 
= 2^b-\t,x,iy) / g{L{t,x,uj),L{t,x,iy))Fj{u)riAuj)duj 

+2^ / {b~\t,x,uj) -b-\t,x,u))g{L{t,x,u),L{t,x,u))Fj{u)r]Au)duj. 

Next, we estimate the second term in the right-hand side of (3.82). We have: 

g(^L{t,x,cu),L{t,x,iy)^ = g{N{t,x,u) - N{t,x,u),N{t,x,uj) - N{t,x,p)). (3.83) 
Thus, the estimate (3.14) for dt^N and the size of the patch yields: 

\\g{L{t,x,u),L{t,x,v))\\L^^Hu) ^ (3.84) 
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which imphes: 

{b~^{t, X, co) — b~^{t, X, i'))g{L{t, X, co), L{t, x, v))Fj{u)rfj {uj)duj 
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(3.85) 



< / \\h \t,x,uj)-b \t,x,u)\\LooL2^^^)\\g{L{t,x,u),L{t,x,u))\\L^\\Fj{u)\\LiVj{i^)^ 

where we used in the last inequahty (3.84), the estimate (3.12) for d^^b, Cauchy-Schwarz 
in u, the size of the patch, and (3.51) for Fj(u). (3.82) together with (3.85) and the 
estimate (3.12) for d^b yields: 



(3.86) 



< 2^ 



g{L{t, X, uj), L{t, X, u))Fj{u)r]j{u)du 



+ ^7 



L-^{M) 



Next, we estimate the right-hand side of (3.86). Using (3.83), the decomposition (3.17) 
for N — N', and arguing as in (3.77), we obtain: 



g(L(t, X, u), L{t, X, iy))Fj{u)r]j{u)du 
Together with (3.86), this implies: 



< 



{u - uyFj{u)r]'^{u)du 



Q^J,{u)7j^{uj)du 



< 



i2i{co-u)fF,iuXico)dw 



+ 7- • (3.87) 



Finally, we need to estimate the first term in the right-hand side of (3.87). We will 
rely on the energy estimate for the wave equation. 

3.4.3. Step3: The energy estimate for the wave equation. Recall from (3.8) 
that: 

DgU — b~^trx- 

Thus, we have: 

□g (^J^ {u - vfFj{u)ri'<{uj)duj^ = b-'\t,x,uj)trx{t,x,uj){22{uj - iy)fFj{u)r]^uj)duj. 



Arguing as in (3.77), we may replace b Hrx by 1: 

b~^{t,x,uj)trx{t,x,uj){2i{uj - p)fFj{u)rij{uj)duj 



(3.88) 



S2 



< 



{2^uj-u)yFj{uX{uj)duj 
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which together with (3.88) imphes: 



S2 



{co - vfFj{u)rf<{(jj)duj 



(3.89) 



L'\M) 



< 



{2i {CO -u)fF,{uX{uj)du 



Let (j) be the scalar function in the left-hand side of (3.88), i.e.: 

= / {u;-ufFj{uX{uj)duj. 

Then, using the energy estimate for the wave equation (2.86) we obtain: 
< l|V0(O,.)|| 

L2(So) + 11^0(0' •)lli2(Eo) + IPg0||L2(>l)||D0||i:,2(_^) 



+ 



M 



(3.90) 



where tt is the deformation tensor of T 

and Qai3 on A4 is the energy momentum tensor associated to 

The control of the parametrix at initial time in [43] (see also Chapter 5) yields 

||V0(O, .)|U2(s„) + ||T0(O, .)|U.(s„) < 7j. (3.91) 

Next, we consider the last term in the right-hand side of (3.90). Prom the maximal 
foliation assumption, tt is traceless, so that 



Using (3.14), one obtains 
and thus 

Q,^7r"'3~7r^^^^(D0)2 + ... . 
It turns out that we have a trace estimate for tcn^n^ (see details in [45]): 

which together with (3.92) implies 



(3.92) 



(3.93) 



< 
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Finally, (3.89)-(3.93) implies 



\B<l>U.^M)<Sm\\Ll^,.^Lr+7j: 



which is equivalent to 



{2Huj-u))'Fj{u)r,Uuj)dco 



< 



L2{M) 

In view of (3.94) and (3.87), we obtain: 



/ {2Huj-u)fFj{uX{uj)duj 



7" 2 r cx 



S2 



Fj{u)riUuj)duj 



< e 



+ 7j- 
(3.94) 



J 2 r oc 



r 2 r oc 



and thus: 

/ F,iuX{uj)(Lj 

which is the desired estimate (3.80). 

3.5. Almost orthogonality in angle 

We have to prove (3.38): 

\\Ejf\\hiM)<J2\\^jf\\lHM)+^"^]- 



This will result from an estimate for: 



E'^f{t,x)E-'f{t,x)dM 



M 



(3.95) 



(3.96) 



Let us introduce integration by parts first with respect to tangential directions, and 
then with respect to L. 

3.5.1. Integration by parts. 

3.5.1.1. Integration by parts in tangential directions. By definition of y, we have y/i = 
V/i — {'VNh)N for any function h on Ej. In particular, we have y/{u) ~ and ^(m') — 
V~^N' - b'~^g(N', N)N. Now, since g(A^' - g(A^, N')N, N') ^ 1 - g{N', and Vu' = 
b''^N', we deduce: 



ji\u—iX'u' 



X'g{N' - g{N, N')N, Vu') 
ib' 



N'-g{N,N')N[G 



iXu~iX u 



(3.97) 



VAr'-g(Ar,A?"')iv(e 



iXu—iX'u' 



A'(l-g(iV',iV)2) 

where we have used the fact that N' — g{N, N')N is a tangent vector with respect of the 
level surfaces of u. We consider an oscillatory integral of the following form: 



/ / 



h{t, x)Fj {u)Fj {u')r]'<{uj)rjf {u')dujduj' dM, 



2 vS2 



66 3. CONTROL OF THE ERROR TERM 

where h is a scalar function on JH. Integrating by parts once using (3.97) yields: 

h{t,x)Fj{u) Fj {u')r]^ {uj)rf<' {uj' ) dwdw'dM 

Jm i§2x§2 1 - g{N,N'y 

X (u) F.j {u')r]'^ {u^' {tj')dwdw'dM, 

where we only kept the term where the derivative falls on h, and neglected for the sim- 
plicity of the exposition the terms when the derivative falls on the denominator of the 
right-hand side of (3.97) or on the volume element oi M.. In view of (3.13), we have: 

N' - g{N, N')N ^ N' -N ^\uj-uj'\r.\u- (3.98) 

and: 

1 - g(7V, TV') = ^ ^ - |cj - cjf r^\u- u'f, (3.99) 

and we thus obtain: 

h{t, x)Fj{u)Fj{u')r)'^{uj)r)^' {uj')dujdu'dM (3.100) 



1 

= i 



b'f{h) Fj (u) Fj {u'X i^K' {oj')dujduj'dM + 



2i z/ - z/' 



Remark 3.18. In the formula (3.100), we neglect two types of terms for the simplicity 
of the exposition. First, we neglect the term when the derivative falls on the denominator 
of the right-hand side of (3.97) or on the volume element of M.. Next, make the following 
approximation: 

N' - g{N, N')N 1 
1 -g(iV, N'y 2j\u-u'\' 
In the actual proof, we use (3.99) to derive the following expansion: 

^ f N-N, y f N' -N,, Y\ 
l-g(7V,7V')2 - |7V,- Ar,,p \2^^''' \\N^-N^^\) \\N^ - N^>\) ) ' ^^'^^^^ 

for some explicit real coefficients Cpg such that the series 

p,g>o 

has radius of convergence 1. Then, (3.100) corresponds to the first term in the expansion 
(3.101) with the additional simplification which consists in replacing \Ni, — Nyi\ with — 
again in view of (3.99). While these approximations greatly simplify the exposition, they 
still allow us to exhibit typical terms in the proof of the almost orthogonality in angle. 
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3.5.1.2. Integration by parts in L. Next, we also introduce integrations by parts with 
respect to L. Since L{u) — and L{u') — 6'~^g(L, L'), we have: 

We consider an oscillatory integral of the following form: 

h{t,x)Fj{u)Fj{u)7]^{uj)r]'^'{u)dujduj'dM, 
where /i is a scalar function on M.. Integrating by parts once using (3.97) yields: 

h{t, x)Fj {u)Fj {u')ri^{ujX' {u')dujduj'dM 

b' 



= -i2- 



— — (L(/i) + . • . )Fj{u)Fj{u')v^{u)vf{u;')dwdu;'dM, 

where we only kept the term where the derivative falls on h, and neglected for the sim- 
plicity of the exposition the term when the derivative falls on the denominator of the 
right-hand side of (3.102) or on the volume element. Using the fact that: 

g(L,L')^-l + g(N,N') (3.103) 

together with (3.99), and keeping only the first term in the expansion (3.101), with the 
additional simplification which consists in replacing \Ni, — N^i\ with |i/ — we obtain: 

h{t, x)Fj{u)Fj{u)7]'^{uj)7]'^'{uj')dujduj'dM (3.104) 

^ [ [ b'L{h)Fj{u)Fj{u')ri'^{uj)ri'^'{uj')du;duj'dM + --- . 

3.5.2. Presence of a log-loss. Let us explain why proceeding directly by integration 
by parts in (3.96) results in a log-loss. Let us define as: 



2^|z/ - u 



Sj,,y^ [ E^f{t,x)E-'f{t,x)dM. 

JM 



We have: 



^,>,.' = / / / / e^'"-^''"'rHrx6'-Hrx'ciA1 

J§2xS2 Jo Jo \Jm J 

xr/J(w)77/(cj7?,/^(2-^A)?,y(2--^A')/(Aw)/(AV)A2A'^dAdA'da;da;'. 

We integrate by parts tangentially using (3.100). Consider the term where the tangential 
derivative falls on trx, which is of the form: 



-1 p p+oo p+oo / p \ 

ji / / / / e'^''-'^''''b-ytrxb'~hrx'dM] 

xr]'^{u)r]f{u')ilj{2-^X)ilj{2-^X')f{Xuj)f{X'uj')X^X'^dXdX'dwduj'. 
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Since Ly/trx is the only derivative of y/trx for which we have an estimate, our next 
integration by parts must be with respect to L, that is we use (3.104). Consider the term 
where the L derivative falls on trx', which is of the form: 



1 



/ / / / e^^"~'^'"'6-^ytrx6'"^L(trxO(^A^ (3.105) 

is2xS2 ,/o Jo \J M J 



Now, note in view of (3.103), (3.99) and the estimate (3.14) for d^N , that: 

g(L, L') ~ , g(L, e'j,) = g(L-L', e'^) ~ and g(L, L') = -2+g(L, L') ~ -2. 

Thus, decomposing L on the frame L', A^', e^, we obtain: 

L L' +\y -u'\f' + \y -v'\^N'. (3.106) 
Together with (3.105), we finally obtain the sum of three terms: 

Sj,uy (3.107) 

= £j^i,y [1] + Sj^i,y [2] + Ej^i,y [3] 

-I r f+oD p+oo / n 

e'^''^'^''''b-ytrxb'~^L'{tTx')dM 



22j|i, _ i,/|3 

xr]^{u)r]f{u')ij{2-^X)ilj{2-^X')f{\u)f{\'u')\'^\'''d\dX'dwduj' 



■^xS^JO Jo \JM 

I, .f\\2\l2. 



+ ^ / / / / e^'''-^''"'?^~'ytrx6'"Y'(trx')rfA^ 

xr]^{u)r]f{uj')ij{2-^X)ilj{2-^X')f{Xuj)f{\'uj')\^\'^dXdX'dwduj' 

-1 /■ /'+00 /•+00 / P 

+ 7i / / / / e'^«-^'"'6-iytrx6'"'iV'(trx')ciM 

X77j(c^)r7/(a;')V'(2"^A)V'(2-^A')/(Ac^)/(AV)A2A'^dAdA'dc^da;'. 

We consider the second term in the right-hand side of (3.107) which is of the form: 
W[2] = ^^-^;^^(^^riytrxF,(x.),7j(a;)da;^ 



X b'~ytrx'Fj{u'X'{u;')dw'^ dM. 
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We claim that such a term leads to a log-loss. Indeed, we have: 
l^w'[2]| 



< 



< 



1 



22j|j,_ j,/|2 
1 



b-ytrxFj(u)ri^(uj)duj 
\b-ytrxFj{u)\\mM)Vj{^)duJ 



I, J\ 



< 



22J|Z/- Z/'|2 



(3.108) 



(2i\u-u'\y 

where we used in the last inequality Cauchy-Schwartz in cu and uj' which gains the square 
root of the volume of the patch, the estimates (3.10) for b and trx, and the estimate (3.51) 
for Fj{u) and Fj{u'). This leads to a log-loss since we have: 

1 



T?(2./>-.'|)^ 



(3.109) 



Indeed, note that v' runs on a lattice on of basic size 2 ■'/^ so that (3.109) corresponds 
to the sum 

/ez2,i<|i|<2J/2 



|/| 



3.5.3. Strategy of the proof of Proposition 3.9. Let us explain informally the 
strategy of the proof. As we noticed in the previous section, the second term in (3.107) 
contains a log-loss. Let us start by showing that the first and the third term in the 
right-hand side of (3.107) do not contain a log-loss. 

3.5.3.1. Control of the first term in the right-hand side of (3.107). We have 



uM{u)du 



X 



b''^L\tix')Fj{u'X\Lj')duj'^ dM. 



In view of the Raychaudhuri equation (3.7), we have: 

^w) = -ix'r+-- 
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where we keep only the worst term. Thus, we obtain 
1 



22j|i/-I/'|3 



M \J§ 



b-ytTxFj{uX{uj)du?j (^j^ b'-^\)('fFj{u'X'{uj')cLu'^dM. 



Let us decompose: 

b'-'\x'\' = b-,'x.'X + {b'-'\x'\' - h-}xu'Xl (3.110) 
and let us assume for the moment that we can control the second term in (3.110). Then, 
we are led to control: 



22i|z/ 
We have: 



< 



< 



2^3\u - i/'|3 
1 



/ b-'xftrxF,{u)7^^{uj)duj b-}x.' f Fj{u'X' {uj')du;' 



\\b-'xftrxFj{u)\\L2^M)Vj{^)duJ 



T.' I t 



Fj{u)r]j' {u')du' 



L2 



u I ,x' I 



< 



22j|i/-i/'|3 L2 



\\b ^||L-(M)||xlU'^L2||ytrx||Loc^2||Fj(u)||L277}(a;)cia; 



^2 iy r/' 

(22 |z/ — v'\y^ 

where we used the estimate (3.10) for trx, £ and 6, the estimate (3.80), Cauchy-Schwartz 
in the size of the patch, and the estimate (3.51) for Fj{u). Note that the right-hand 
side of (3.111) does not contain a log-loss since: 

^^pE ^o./.u.^ (3-112) 



-Z_.(2.72|^_^.|) 

Remark 3.19. While the estimate obtained in (3.111) is correct, one has to modify 
slightly the method leading to it. Indeed, x does not have enough regularity with respect to 
CO to be able to handle the second term in the decomposition (3.110). The way to overcome 
this is to make use of the decomposition (3.15) for x' 

X = Xi + X2- 

Then, we exploit the fact that, in view of the estimate (3.16), xi has better regularity than 
X with respect to {t,x), while xi has better regularity than x with respect to uj. We refer 
to [45] for more details. 
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^j,uy [3] 



1 



b-ytTxFjiu)riUoo)duj 



X (^j^ b'-^N'{tTx')Fj{u'X'{uj')duj'^ dM. 



Now, wc have: 
1 



< 



< 



22j|i/-i/' 
1 

22i|i/-i/'| 
1 



b-ytTxFjiu)r]^{uj)du^ (^j^ b'~^N'tTx'Fj{u')r]f'{uj')duj'^ dM 

b'-^N'trx'Fj{u')r)^'{uj')dw' 



M \^S2 

b-ytrxFj{u)r]Uco)du; 



L-^{M) 



b-ytrxFj{u)\\mM)Vj{^)dio 



< 



22i|z/ - u' 

X (^jjb'-'N'trx'Fj{u')\\L^M)Vj{^y^' 



X 



7§2 



< 



22(22|l/-I/'|) 



(3.113) 



where we used in the last inequahty Cauchy- Schwartz in lo and lo' which gains the square 
root of the volume of the patch, the estimates (3.10) for b and trx, and the estimate (3.51) 
for Fj{(u) and Fj{u'). Note that the right-hand side of (3.113) does not contain a log-loss 
since we have: 

^, 22 (22 |z/ - z/'l) 

3.5.3.3. A decomposition for E J f . To remove the log- loss exhibited in (3.108) (3.109), 
we rely on a decomposition of trx using the geometric Littlewood-Paley projections Pj. 
We have: 

trx = P<j/2(trx) + 

l>j/2 

which in turn yields the following decomposition for Ej f: 



E^fit,x)= Y.Eff{t,x), 

l>j/2 



(3.115) 
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where: 

Eff{t,x)= [ b{t,x,uj)-'Pitrx{t,x,uj)Fj{u)7i'^{uj)dw, V/ > ^ (3.116) 

and: 

E'^''/^f{t,x)= [ b{t,x,uj)-^P<j/2trx{t,x,uj)Fj{uX{uj)dw. (3.117) 
In order to prove almost orthogonality in angle, i.e. (3.95), we will estimate: 

J2 I Eff(t,x)E;'''^f(t,x)dM . 

3.5.3.4. The mechanism to remove the log-loss. In order to explain the mechanism 
which allows us to remove the log-loss, let us assume for convenience that m < / in 
(3.118). Then, notice first from (3.107), (3.109), (3.112) and (3.114) that the only term 
in the right-hand side of (3.107) which contains a log-loss is the second one, i.e the term 
which contains only tangential derivatives. In order to remove the log-loss, our goal will be 
to always put more tangential derivatives on the lowest frequency, i.e. P^^trx' (as opposed 
to the higher frequency P^trx). This is achieved as follows (see [45] for the details): 

(1) Integrate by parts with respect to L using (3.104). 

(2) One term corresponds to the case where the L derivative falls on the largest 
frequency P/tr^, while the other term corresponds to the case where L falls on 
the lowest frequency P^trx'- For the second term, decompose the L derivative 
on the frame L', A^', e'^ as in (3.106). 

(3) Notice that the terms involving L, V or N' are estimated in the spirit of (3.111) 
and (3.113), and should in principle contain no log-loss in view of (3.112) and 
(3.114). 

(4) Finally, the last term is the one containing the y derivative. This term is the 
only one which contains the log- loss exhibited in (3.109). Now, we have achieved 
our goal since after integration by parts, the tangential derivative fell on Pmtrx' 
which is the lowest frequency. 

Remark 3.20. Due to the decomposition (3.115), we now not only need to obtain 
summability in (u,!/'), but also in {l,m). This creates additional difficulties, in particular 
when estimating the terms £j^i,y[l] and £j,vy^] in (3.107). We refer to [45] for more 
details. 
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CHAPTER 4 



Control of the space-time foliation 

The goal of this chapter is to prove the estimates on the control of the space-time 
foliation by the optical function u which are needed for the proof of Theorem 3.6 (see 
section 3.1.3), i.e. for the control of the error term. Here, we outline the main ideas and 
we refer to [44] for the details. 

4.1. Geometric set-up and main results 

4.1.1. Geometry of the foliation of Ai by u. Recall from section 1.2.1 that the 
space-time M is foliated by space-like hypersurfaces defined as level hypersurfaces of a 
time function t, where T denotes the unit, future oriented, normal to Ej and k its second 
fundamental form. Recall also that it is a solution to the eikonal equation g^'^daudpu — 
on Ai depending on a extra parameter u G S^. The level hypersufaces u{t,x,u) = u of 
the optical function u are denoted by 1-Lu- Let L' denote the space-time gradient of u, i.e.: 

L' = -g'^^dpud^. (4.1) 

Using the fact that u satisfies the eikonal equation, we obtain: 

T>L'L' = 0, (4.2) 

which implies that L' is the geodesic null generator of Hu- 
We have: 

T{u) = ±\Vu\ 

where \Wu\'^ — Ylt=i relative to an orthonormal frame Cj on E^. Since the sign of 

T{u) is irrelevant, we choose by convention: 

T(u) = \Vu\. (4.3) 

We denote by Pt_„ the surfaces of intersection between E^ and Hu- They play a funda- 
mental role in our discussion. 

Definition 4.1 (Canonical null pair). 

L = bL' = T + N, L = 2T - L = T - N (4.4) 

where L' is the space-time gradient ofu (4.1), b is the lapse of the null foliation (or shortly 
null lapse) 

b-^^-<L',T>^T(u), (4.5) 
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and N is a unit normal, along T,t, to the surfaces Pt^u- Since u satisfies the eikonal 
equation g'^^daudjsu — on A4, this yields L'{u) — and thus L{u) — 0. In view of the 
definition of L and (4.3); we obtain: 

N=-^. (4.6) 

Remark 4.2. u is prescribed on Sq as in [42]. For any {0,x) on Sq, L is defined as 
L = T + N where T is the unit normal to Sq at {0,x) and N = — Vm/|Vm| at {0,x), and 
b is defined as b"^ = |Vm|. Let nxif) denote the null geodesic parametrized by t and such 
that Kj;(0) = (0, x) and k^(0) = b~^L. Then, we claim that 

K'.it) = b{K,{t))-^L,^(^t) for all t. (4.7) 
Indeed, L' — b~^L is the geodesic null generator ofHu (see (4.2) 

Definition 4.3. A null frame 61,62,63,64 at a point p G Pt^u consists, in addition to 
the null pair 63 = L,e4 — L, of arbitrary orthonormal vectors ei, 62 tangent to Pt^u- 

Definition 4.4 {Ricci coefficients). Let 61,62,63,64 be a null frame on Pt,u above. 
The following tensors on St^u 

Xab =< 0^64, 6s >, =< DAe3, cb >, (4.8) 

Ca = ^ < 0364, eA>, = ^ < 0463, ca >, 

= ^ < D363,6A > . 

are called the Ricci coefficients associated to our canonical null pair. 
We decompose x and x into their trace and traceless components. 

trx = S''''xAB, t^ = S^''xAB^ (4-9) 

XAB = XAB - ^trxgAB, Xab ^ Kab ~ \^^Sab, (4.10) 

Definition 4.5. The null components of the curvature tensor R of the space-time 
metric g are given by: 

aab = R{L,ea,L,eb) , = ^R(6„, L, L, L), (4.11) 

p = ^R{L,L,L,L), a^^*R{L,L,L,L) (4.12) 

= ^R(e„, L,4^), a„;. = R(L,e„, L,e6) (4.13) 
where *R denotes the Hodge dual of R. 

Observe that all tensors defined above are Pj^-tangent. 
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Remark 4.6. Note that a is the only null component which does not contain a con- 
traction o/R with L. With the notation of Chapter 2 (see for instance (2.28) ), we have: 

R-L = (a,/3,p, a,p). 

Definition 4.7. We decompose the symmetric traceless 2 tensor k into the scalar 5, 
the Pt,u-tangent 1-form e, and the Pt,u-tangent symmetric 2-tensor rj as follows: 

kNN = S 

kAN = (4.14) 

kAB = TlAB- 

The following Ricci equations can be easily derived from the definition of T, the fact 
that L' is geodesic (4.2), and the definition (4.8) of the Ricci coefficients (see [9] p. 171): 



(4.15) 



0^64 = 


XAB^B — ^A^i, 


0^63 






0464 = 


-Se4, 


0463 ^ 


= 2C^eA + 5e3, 




0,364 = 


"^Ca^a + (5 + n""^VArn)e4, 


0363 : 


= 2^^eA -{S + n 


VArn)e3 


046^ = 


f4eA + C^e4, 


036^ 


= fs^A + (Aes + 






^ 1 1 

-- fBSA + ^XAB 63 + ^Xab ^4 









where, ^3, ^4 denote the projection on Pt^u of D3 and D4, ^ denotes the induced covariant 
derivative on Pt^„ and S, e are defined by: 

S^S-n-^N{n),eA^eA-n-^VAn. (4.16) 

Also, 

Kab ^ ~ 

= -^A, (4.17) 

= eA + n'^fAn- Ca- 



4.1.2. Null structure equations. Below we write down our main structure equa- 
tions (see [9] chapter 7 or [44] for a proof). 
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Proposition 4.8. The components trx, X, C ^'^^ the lapse b verify the following equa- 



tions^: 



L{b) = -bS, (4.18) 

L{trx) + ^itrxT = -\x\^ - Strx, (4.19) 

fiX + trxx^-Sx-c^, (4.20) 

%Ca + \{trx)CA = -{eB + Cb)xab - \trxtA - /9a, • (4.21) 



Remark 4.9. Equation (4.19) is known as the Raychaudhuri equation in the relativity 
literature. 

To obtain estimates for x- we may use the transport equations (4.19) (4.20). However, 
this does not allow us to get enough regularity. Instead, we follow [9] [20] [22] and 
consider (4.19) for tr^ together with an elliptic system of Hodge type for x- 

Proposition 4.10. The expression {dfvxjA — )^^Xab verifies the following equation: 
{dfvx)A + XABeB = \{yAtrx + eAtrx) - Pa- (4.22) 
See [9] chapter 7 or [44] for a proof. 

Finally, we consider the control of ( and Ltrx- To this end, we follow again [20] [22]: 
we derive an elliptic system of Hodge type for ( and a transport equation for Ltr^. 

Proposition 4.11. We have: 

L{trx) + \trxtrx = 2c//< + (5 + n'^V Nn)trx - x • X + 2C • C + 2p. (4.23) 
Also, the expressions dfvC = y^Cs (^''^d cuflC, =E^^ J^aCb verify the following equations: 

dH=l(^l^ + ltrxtrx + X-X-m''^ - P, (4-24) 

cu/-/C = -^XAx + <7, (4.25) 

where for F, G symmetric traceless Pt^u-tangent 2-tensors, we denote by F AG the tensor 
F A Gab —&ab FacGbc- Finally, setting, 

A* = Utrx) - (5 + u'^Vnu) trx (4-26) 



^which can be interpreted as transport equations along the null geodesies generated by L. Indeed 
observe that if a P(^„-tangent tensor H satisfies the homogeneous equation ^411 = then 11 is parallel 
transported along null geodesies. 
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we find 

- trx (2dfvC + 2C ■ C + 4(e - C) ■ n-yn - 2S{S + n" V^n) + 4p (4 27) 

- \trxtrx + 2|e|2 + 3|x|' + 4% • ?7 - 2\n-'N{n)\'^^ . 
See [20] or [44] for a proof. 

4.1.3. Commutation formulas. We have the following useful commutation formu- 
las (see [9] p. 159): 

Lemma 4.12. Let Ua be an m-covariant tensor tangent to the surfaces Pt_„. Then, 
fsfAUA-f^fBUA = XBcfcUA-n-yBnfiUA (4.28) 

i 

fsfzUA - fsfsUA = XecfcUA - i^f^UA - b-^VebfsUA + J^^-XaM^ 

i 

+ XBci^^-XA,B^C + XscCA.+ ^A,C*^j,)UA,..C..Ar^: (4.29) 

%f4UA-f4f3UA = -SfsUA + i6 + n-W^n)%UA + 2{CB-CJfBUA{^.m 

+ 2 E(^Cc - CcCa+ &a,c *<t)Ua,..c..a^. 

i 

4.1.4. Bianchi identities. In view of the formulas on p. 161 of [9], the Bianchi 



equations for a, ^, p, a, p are: 

fL{/3) = d^a-5/3 + (2e-e) -a (4.31) 

fM = fp+{far + 2x-P + {S + n-'VNn)(3 + ^-a + 3{Cp + XcT) (4.32) 

L{p) = d^/3-^x-«+(e-2e)-/3 (4.33) 

Up) = -dyV^-^x-« + 2e-/5 + (e-2C)-^ (4.34) 

L{a) = -cu/]^ + ^x*Q; + (-e + 2e)*/3 (4.35) 

L{a) - -cu/]^- 1^*^-2^/3 +(6 -20*^ (4.36) 

= -yp+(y<j)* + 2x-^ + 5^-3(Cp- V) (4.37) 
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4.1.5. Main results. We introduce tlie curvature flux TZ relative to the time 
foliation: 

^ = {Minn.) + mihin.) + Mhin.) + Minn.) + mlnn.)) ' ■ (4-38) 

In view of Remark 4.6, we have TZ = ||R-L||x,2(^^). Thus, the we may rewrite the bootstrap 
assumptions of Chapter 2 on R as: 

||R'||l-l2(so < Me, sup7^ < Me. 

u 

To ease the notations, we drop the bootstrap constant M: 

llR-IUfLaCEt) < SUp 7^ < £. (4.39) 

u 

The goal of this part is to control the geometry of the null hypersurfaces Tiu of u up to 
time t — 1 when only assuming the smallness assumption (4.39). 

Remark 4.13. In the rest of the chapter, all inequalities, except the ones of Theorem 
4-16 below, hold for any u with the constant in the right-hand side being independent of 
u. Thus, one may take the supremum in u in these inequalities. To ease the notations, 
we do not explicitly write down the supremum in u in these estimates. 

ti is a solution to the eikonal equation g^^^daudpu = on depending on a extra 
parameter a; G Now, for u to be uniquely defined, wc need to prescribe it on Sq 
(i.e. at t = 0). This issue has been settled in [42] (see also Chapter 6). From now on, 
we assume that u is the solution to the eikonal equation g'^^daudjsu = on Al which is 
prescribed on Eq as in [42] . 

Remark 4.14. In the rest of the chapter, all inequalities hold for any a; e with the 

constant in the right-hand side being independent ofcu. Thus, one may take the supremum 
in u everywhere. To ease the notations, we do not explicitly write down this supremum. 

We define some norms on For any 1 < p < +oo and for any tensor F on 1-Lu, we 
have: 



Wo JPt,u I 



where dijbt,u denotes the area clement of Pt^u- We also introduce the following norms: 
MF) = \\F\\ 

U2{F) = M(F) + \\Tf\\l2(^^^ + wyyLFWmnuy 

Let x' a coordinate system on Po,«- By transporting this coordinate system along the 
null geodesies generated by L, we obtain a coordinate system (t,a;') of %. We define the 
following norms: 



I^IIl-l?= sup (! \F{t,x')\Ht] 
x'ePo,u \Jo J 
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sup \F{t,x'))\ 

0<i<l 



The following theorem investigates the regularity of u with respect to {t,x): 

Theorem 4.15. Assume that u is the solution to the eikonal equation g'^^daudjsu = 
on M. such that u is prescribed on Eq as in [42]. Assume also that the estimate (4.39) 
is satisfied. Then, null geodesies generating Hu do not have conjugate points and distinct 
null geodesies do not intersect. Furthermore, the following estimates are satisfied: 

\\n - + ||Vn||iooi2, + ||V'n|| rooi2, + || VDrnH roc^2, < e, (4.40) 

^fl{k) + WfLeh^n^) + II ^(5)IU2(«„) + llelL^L? + PIIl-l? < (4-41) 

||6 - 1|U^ + AA2(6) + II m\\Li,lr < ^, ' (4.42) 

Wtrxh^ + WftrxhiL^ + II ^^^IIl2,l- < e, (4.43) 

I|x||l2,l- + ^(x) + WfLXhHHu) ^ (4.44) 

\L-^,Lr+WO<s. (4.45) 



Wc introduce the family of intrinsic Littlcwood-Paley projections Pj which have been 
constructed in [24] using the heat flow on the surfaces Pt^u (see also section 3.1.5). This 
allows us to state our second theorem which investigates the regularity of LLtrx and 

fLC- 

Theorem 4.16. Assume that u is the solution to the eikonal equation g°'^daud^u — 
on Ai such that u is prescribed on Sq as in [42]. Assume also that the assumption (4.39) 
is satisfied. Then, there exists a function A in L^(]R) such that for all j > 0, we have: 

WPjLLtrxhHn.) < 2^e + 2ix{u), (4.46) 

and 

\\PjfL{C)\\LHH.)<e + 2-iX{u). (4.47) 

The following theorem investigates the regularity with respect to the parameter a; e 

Theorem 4.17. Assume that u is the solution to the eikonal equation g,°'^do.ud^u — 
on M. such that u is prescribed on Sq as in [42]. Assume also that the estimate (4.39) is 
satisfied. Then, we have the following estimates: 

||-9.iV|U- < 1, (4.48) 

||Da,7V||^2,^oc + \\d^b\\L^ + Wfd^LlLr + II^-x||l^,lc« + \m\\Ll,Lr < e. (4.49) 

Furthermore, we have the following decomposition for x- 

X = Xi + X2, (4.50) 

where xi o-nd xi o.re two symmetric traceless Pt^u-'to.ngent 2-tensors satisfying: 

M(Xl) + liyLXl||L2(«„) + l|9c.Xl|U?°L^+M(X2) + ||yLX2|U2(«„) + ||9^X2||L-L2, < S (4.51) 
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and for any 2 < p < +oo, we have: 

IIxiIIl^l- + \\du;X2\\^p^^- + ||te||L6-(H„) + liytelU2(«„) < £, (4.52) 

where for any real number a, a_ = a — 5 for any 8 > 

Remark 4.18. Notice from (4.51) that Xi o-nd X2 have at least the same regularity as 
X- Now, the point of the decomposition (4.50) is that both xi md X2 have better regularity 
properties than £. Indeed, in view of (4.52), xi has better regularity with respect to {t, x) 
while X2 has better regularity with respect to uj. 

Next, the following theorem contains estimates for second order derivatives with re- 
spect to UJ. 

Theorem 4.19. Assume that u is the solution to the eikonal equation g,'^^da.udpu — 
on M. such that u is prescribed on So as in [42]. Assume also that the estimate (4.39) is 
satisfied. Then, we have the following estirnMes: 

WdlNW^lLr ^ 1, (4-53) 
||y^n(9»|U.(«„)<£. (4.54) 

||p,-yz,n(9^7V)||^.^. + ||p,.^(9^x)llLo^L^ + ||P,n(9^c)IU?L^ < 2^£, (4.55) 

where p is any real number such that 2 < p < +00, and where U denotes the projection 
on Pt,u-to,ngent tensors. 

Finally, we need to compare quantities evaluated at two angles uj and i/. The following 
decompositions are used in sections 3.4 and 3.5 

Theorem 4.20. Let u and v in such that — v I < 2-2. Let u = u{., u), N = 
N{.,uj) and Ni, — A^(.,i/). For any j > 0, we have the following decomposition for N — N^. ■ 

2-2{N-N^)^ Fi + Fi (4.56) 
where the tensor Fl does not depend on uj and satisfies: 

\\Fl\\L^<h 

and where the tensor F^ satisfies: 

We also have following decomposition for trx-' 

trx = fi + fi (4.57) 
where the scalar fl does not depend on u and satisfies: 

\\fl\\L^<e, 

and where the scalar /| satisfies: 

11/2 IU-L2(^^) < £2-2. 
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Let us conclude this section by mentioning several ingredients of [44] that have been 
omitted here for the sake of simplicity: 

• estimates for the transport equations along L, and the elliptic systems of Hodge 
type on u involved in the null structure equations 

• embeddings on "H^, and „ 

• geometric Littlewood-Paley projections and Besov spaces on 

• control of the Gauss curvature of Pt^u 

• Bochncr inequalities on Ej and Pf.u 

• estimates for various commutator terms of the type: [Di,y], [D^, y], [DL,Pj], 
[Dl,P,], ... 

4.2. Regularity of the foliation with respect to {t, x) 

In this section, we outline the main ideas of the proof of Theorem 4.15. We assume 



the following bootstrap assumptions: 

11^ - 1|U-CK0 + \\b - l|U-(«„) < ^, (4.58) 

\\Vn\\L^Ll, + W^^MLrLl, + WVBttiWl^lI, +M2{h) + \\L{h)\\L2^^Lr < De, (4.59) 

J^i{k) + WfLeWm-Hu) + I|Dl5||l2(^„) + llelUoo^i + \\S\\l^l^^ < De, (4.60) 

l|trx|U-(H„) + \\ftTx\\Ll,L^ + II Ltrx|U2,^^^ < De, (4.61) 

llxl|L^,Lr +-^i(x) + WfLXhHH^) < De, (4.62) 

\\C\\L^L^+^fl{0<De, (4.63) 
where D > is a large enough constant. We will improve on these estimates. 



4.2.1. Non intersection of null geodesies on Hu- The control we obtain on the 
geometric quantities associated to our foliation is only valid as long as there are no con- 
jugate points and null geodesies do not intersect. The goal of this section is to prove that 
this holds at least until t = 1. In addition to the bound (4.39) on the curvature tensor R 
of g, we make the following regularity assumption on g. There exists a coordinate chart 
such that 

||g||c2(^) < M, (4.64) 

where M is a very large constant. 

Remark 4.21. The assumption (4.64) is only used to prove the absence of caustic and 

that null geodesies do not intersect at least until t = 1, which is a qualitative property. On 
the other hand, we only rely on the bound (4.39) on R to prove the various quantitative 
bounds of Theorems 4-15, 4-^6, /^.i7 and 4-19. 

For {0,x) in Eq, recall the definition in Remark 4.2 of the null geodesic K^(t). For all 
< i < 1, let $t : E ^ E^ defined by $t(0,a;) = K^{t). We have ^o{0,x) = {0,x) on Eq. 
We define to >0 as the supremum of < i < 1 such that is bijective from Eq to E^. 
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Remark 4.22. ^45 long as < t < to, there are no conjugate points and no distinct 
null geodesic intersections. Thus, we may assume that the u-foliation exists and satisfies 
the bounds (4.58) -(4.63) given by the bootstrap assumptions. Furthermore, we may assume 
the identity (4.7) for the null geodesies Kx{t). 

Our goal is to show that we have in fact to — 1. We proceed in three steps (see [44] 
for the details): 

Step 1. As noticed in Remark 4.22, the bound for trx given by (4.61) holds for 
< t < to- Furthermore, using the Raychaudhuri equation (4.19) and the bound 
(4.64), we obtain the existence of a constant 5 > depending on M such that 
the L°° bound for trx given by (4.61) holds lor Q < t < to -\- 5. This control for 
trx allows us to prove that there are no conjugate points on < t < to + ^• 

Step 2. Next, we prove that = UuPt,u for < t < to + ^ where 5 > is a constant 
depending on M. This requires the bound (4.64), and the control it induces on 
forward and backward light cones for small time intervals with a size depending 
on M. 

Step 3. Assume now that < to < 1- In view of Step 1 and Step 2, the only thing that 
can go wrong at t = to is that two distinct null geodesies intersect in . Assume 
by contradiction that this is indeed the case so that there exists (0, a;i) (0, 2:2) 
two points in Eq such that K^^ito) — K^^ito) — {to,Xo)- Since 

= bC^x^it))'^!^^^,^^^, j = 1,2, 

in view of Remark 4.2, the regularity of b and L yields k,'^_^ (to) = k'^^ (to). Prom the 
classical uniqueness result for ODEs, we deduce that 

(t) = Kx^it) for all t. In 

particular, taking t = 0, we obtain {0,xi) = (0,^2) which yields a contradiction. 
Finally, Steps 1, 2 and 3, yield to > 1. In particular, we have: 

On < t < 1, there are no conjugate points and no intersection of distinct 

null geodesies. In particular, u exists on < t < 1 and the bootstrap (4.65) 

assumptions (4.58)-(4.63) hold. Furthermore, St = UuPt,u for all < t < 1. 

4.2.2. Lower bound on the volume radius of S^. In this section, we prove the 
lower bound on the volume radius of E^ given by the estimate (2.32). We use the global 
coordinate system x' — {x^,x'^) on Po,m which has been constructed in [42] (see also 
Proposition 6.13). Transporting this coordinate system along the null geodesies generated 
by L, we obtain a coordinate system x' of Pt^u, which in particular satisfies 

(1 - 0{em' < iAB{p)i^e < (1 + 0{e)m\ uniformly for all p e P,,„, (4.66) 

where 7 is the metric induces by g on Pt^„. We denote by x' this global coordinate system 

on Pt^u- 

Next, we obtain a global coordinate system on E^ as follows. First, recall from (4.65) 
that Ej = yjPt,u so that u is defined on E^. To any p G E^, we associate the coordinates 
{u{p),x'{p)) where u{p) is the value of the optical function u at p, and x'{p) are the 
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coordinate of p in the coordinate system of Pt^u- In this coordinate system, the metric gt 
on Et (i.e. the restriction of g on E^) takes the following form: 

where 7 is the induced metric on Pj^„. Together with the estimate (4.58) for b and (4.66) 
for 7, we obtain the following lower bound on the volume radius of E^ at scales < 1: 

rU^t,l)>\, (4.68) 

which is the estimate (2.32). 

4.2.3. Estimates for the second fundamental form k and the lapse n. We 

first estimate A; on E^. k satisfies the following symmetric Hodge system on E^: 

curl/cjj H,^ii,jT^T , 

V^kij = 0, (4.69) 
trk = 0, 

where cmlkij = |(G-™ Vikmj+ 6^™" Vikmi) and irk — g^^hj. Using an eUiptic estimate 
for the Hodge system (4.69), we easily obtain: 

||VA;|Uoc^.(s^) < e. (4.70) 

Recall from (1.9) that the lapse n satisfies the following elliptic equation on E^: 

An = \k\^n. (4.71) 

Using (4.71) and (4.70), together with elliptic estimates on E^, we improve the estimate 
for n in the bootstrap assTimptions (4.58) (4.59). We also prove the following estimate 
which is needed for the estimate (2.63) 

llVnlUoc(^) < e. (4.72) 

Using (4.71) and (4.70) together with the Sobolev embedding on the three dimensional 
riemannian manifold yields An G L'^L^illt)- Together with elliptic estimates, this 
implies V^n e L^L^(Ej), and thus Vn misses to be in L°°(A^) by a log divergence. 
However, one can overcome this loss by exploiting the Besov improvement with respect to 
the Sobolev embedding on E^. This requires to introduce a geometric Littlewood-Paley 
theory on Ef^. We refer the reader to section 4.4 in [44] for the details. 

Finally, we estimate k on To this end, we use the decomposition of k (4.14) in 
5, e and ?], and obtain a Hodge system for 5, e and rj on "H^. This allows us to derive the 
following estimate 

Mk)<mri) +M(e) +M(5) < e. (4.73) 



^Note that we use a geometric construction based on the heat flow on since we don't have enough 
regularity for the metric in order to use a coordinate dependent Littlewood-Paley decomposition 
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Deriving an estimate for T{6) and yre, together with (4.73), then yields 



\T>LS\\L^L^(nu) + \\f Le\\L^L2^nu) ^ ^- (4-74) 



4.2.4. Time foliation versus geodesic foliation. While we work with a time foli- 
ation, we recall that the estimates corresponding to the bootstrap assumptions on % and 
C have already been proved in the context of a geodesic fohation in [22] [24] [23]. One 

may reprove these estimates by adapting the proofs to the context of a time foliation. 
However, this would be rather lengthy and we suggest a more elegant solution which con- 
sists in translating certain estimates from the geodesic foliation to the time foliation, and 
in obtaining directly the rest of the estimates. More precisely, we wish to obtain the L°° 
bound from trx, and the trace bounds for £ and ( by exploiting the corresponding esti- 
mates in the geodesic foliation. We will obtain the trace bounds for 6 and e by reducing 
to estimates in the geodesic foliation in section 4.2.5. Finally, these trace bounds and the 
null structure equations will allow us to get all the remaining estimates in section 4.2.6. 
We start by recalling some of the results obtained in the context of the geodesic foliation 
in [22] [24] [23]. 

4.2.4.1. The case of the geodesic foliation. Recall that L' = —^^^dpuda is the geodesic 
null generator of "H^. Let s denote its affine parameter, i.e. L\s) = 1. We denote by 
the level surfaces of s in Hu- 

Definition 4.23. A null frame e[,e2,e'^,e'^ at a point p G P^'^ consists, in addition to 
64 — L' , of arbitrary orthonormal vectors 6^,63 tangent to P^^ and the unique vectorfield 
e'g — Jj satisfying the relations: 

g{e',, e',) = -2, g(e'3, e'3) = 0, g{e'^, e[) = 0, g(e^, e',) = 0. 

Definition 4.24 {Ricci coefficients in the geodesic foliation). Let 6^,63,63,64 be a 
null frame on Pj as above. The following tensors on Pj^ 



(4.75) 



x'ab =< De^e^, c'b >, =< De^6'3, c'b >, 

Ca = i < De'^6^, 6^ > 

are called the Ricci coefficients associated to the geodesic foliation. 
We decompose x' and x' into their trace and traceless components. 

trx! = S^'^'ab, trx; = g^'^x!^^. (4.76) 

Xab = Xab - ^^r\:'gAB, = X^s - ^^'^'gAB- (4.77) 
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Definition 4.25. The null components of the curvature tensor R of the space-time 
metric g in the geodesic foliation are given by: 

= R{L\e',,L',ei), = ^R(e;, L', L', L'), (4.78) 

p' = \r(I/,L',IJ,L'), a'^\*R(lJ,L',lJ,L') (4.79) 

= 2^(el,L',L',L'), a:„, = R(L',e;,L',e;) (4.80) 
where *R denotes the Hodge dual of R. 

We now recall the main estimates obtained in [22] [24] [23] . We have: 

l|trx'||Loc(«„) + Wx'hl^Lr + IIC'I|l^,l- < e (4.81) 

and 

\W\\Ll,Lr+Kix')+KiC')<e, (4.82) 
where the norm J\f[ is given by 

KiF) = \\F\\L^iH.) + Wf'FhHH.) + WfuFU^H.)- 

Remark 4.26. Note that the norm L°°{'Hu) does not depend on the particular foliation. 
Now, this is also the case for the trace norm L^iLf'. Indeed, recall the definition of the 
null geodesic Remark 4-2. Then, we have: 



\F\\l^^2= sup f \F{K^{t))\^dt^ sup f \F{K^{s))\^n~^h-^ds ^ \\F\ 

^' * (0,i;)eSoJo (0,i;)eSo7o 
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where we used the fact that ^ — n ^ and the fact that nb 1 by the bootstrap assump- 
tion (4.58). 

In the next section, we will obtain the estimates corresponding to (4.81) in the time 
foliation. For now, we conclude this section by recalling the definition and some properties 
of the Besov spaces constructed in [22] [24] [23]. For P^,^-tangent tensors F on Tiu, 
< a < 1, we introduce the Besov norms: 

= J]2^'^ sup \\p;f\\l2^p,j + sup WP'^.F^^^p^j, (4.83) 
'rrz o<s<i ' o<s<i 

j>0 — — 

= J2^^^\\PjPhHnu) + \\P<oF\\LHn.) (4.84) 
j>o 

where Pj arc the geometric Littlewood-Paley projections on the 2-surfaces Pg^- Using the 
definition of these Besov spaces, we have (see [22] [24] [23]) 

llxlls'o <£■ (4-85) 
Furthermore, we have for scalar functions on Hu (see [22] section 5): 

< < WfhrL^ + Wf'fWs'O- (4.86) 
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Finally, we have the following version of the sharp classical trace theorem (see Corollary 
4.21 in [44] for a proof). 

Proposition 4.27. Assume F is an P'g ^-tangent tensor which admits a decomposition 
of the form, f'F = Af^'P + E. Then, 

\\F\\L-Li<K{F) + M[{P){\\A\\r.^ + WAW^iLr + liyL'^IU^,Lsc) + \\E\\^,o. (4.87) 

4.2.4.2. Estimates in the time foliation. In this section, we obtain the L°° bound for 
trx, and the trace bounds for x ^-iid C by relying on the corresponding estimates in 
the geodesic foliation (4.81). We start by establishing the relation between the Ricci 
coefficients in the time and in the geodesic foliation. Recall from (4.4) that L = bL'. 
Since (ci, 62) and {e[, e'2) are both orthonormal vectors in the tangent space of T-Lu which 
are both orthogonal to L, we may chose these vectors such that there is a tensor F' on 
satisfying: 

eA = e'^ + F'^L',A = l,2. 
We then easily express L in the frame (L', Lf, e'^). Finally, we have the following relations: 

L = bV, 

eA = e'j, + F'^L', A = 1,2, (4.88) 
L = b-^lJ + ^b-^F^A + b-^\F'\^L'. 

Next, using the definition (4.8) and (4.75) of the Ricci coefficients respectively in the 
time and geodesic foliation, and the identities (4.88), we easily obtain 

X = fcx', trx = 6trx', X = U = C'a + ^acF'c- (4-89) 
(4.89) together with the bootstrap assumption (4.58) and the estimate (4.81) yields: 

l|trx||L-(«,) < ||6|Uoo(^„)||trx'||L-(«,) <e, 

llxlU-L? < ||6||loc(-^„)||x'||l2,l- < e, (4.9O) 

IICIIl^l? < IIC1U^L? + llx1U°?L?ll^1U- <£ + £||^1U-, 

x' f x' ° x' ° 

where we have used the fact that the trace norms Ll,L^ and L^iLf^ are equivalent by 
Remark 4.26. 

In view of the trace estimate for ( given by (4.90), we need to estimate ||F'||^oo. To 
this end, we estimate J/'F'. Using the definition (4.8) of x and (4.75) of x\ and the 
identities (4.88), we obtain: 

g(De^F',e's) = -^x'^^ + --- , 

where we only kept the main term. Together with the estimate for x' (4-85), this yields 

\\f'F'\\s,o<De 

which together with (4.86) implies: 

llF'lUoc < De. (4.91) 
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In particular, (4.90) and (4.91) imply: 

IICIIl^l? < e. (4.92) 

Note that (4.90) and (4.92) are improvements of the corresponding estimates in the boot- 
strap assumptions (4.61)-(4.63). 

4.2.5. Trace norm bounds for 6 and e. The goal of this section is to improve the 
estimate for ||(5||2,o<^l2 and ||e|| ^00^,2 given by the bootstrap assumption (4.60), where S and 

e are defined in (4.16). Let us first define and kiA- 

kLL = -g(DiT, L), kLA = -g(DiT, e^), A = 1, 2. (4.93) 
Then, using in particular the definition (4.16), we have: 

I = kLL and Za = kLA- (4.94) 

We also define kyv and kL'A'- 

kuL' = -g(Di.r, L'), kL'A = -g(Di.r, e'j,),A^ 1, 2. (4.95) 

Then, the relations (4.88) between L, ei,e2 and L', 6^,62 together with the definitions 
(4.93) and (4.95) yield: 

kLL = h'^kL'L' and kLA = bkL'A + hF'j^kL'v- (4.96) 
Thus, (4.94) and (4.96) imply: 

II^IIl-L? ^ \\bkL'L'\\L^,Ll<\\kL'L'\\L^,Ll (4.97) 

II^||l°7L2 ^ \\bkL'A\\L°°,L'i + ll^-^yl^L'L'IU°?ii ~ || ^L'L' || L°?L2 + \\k U a\\ L°^, L'i 

where we used the bootstrap assumption (4.58), the L°° bound for F' (4.91) and Remark 
4.26. 

In view of (4.97), it is enough to bound the trace norms ||A;l'l'IU°?l2 and ||A;i:,'^||£,oo^2. 
To this end, wc would like to apply the trace estimate (4.87), which requires to show that 
y/'kL'L' and yki'A admit a decomposition of the form, A^^l'P + E. We only discuss the 
estimate for kL'u, and we refer the reader to [44] for /cl'A- We have: 

f'e'^kL'L' = -B,'^g{BL'T,L') = -g{B,'BL'T,L')-g{BL'T,D,>^L'). 

Introducing the commutator term [De/^,Di/], and decomposing the corresponding com- 
ponent of R, we obtain 

%,^kL'L'^-b-'Fi,a'^^ + -- - , (4.98) 

where we only kept a typical term for simplicity. Relying on the Bianchi identities, the 
following decomposition for a' was obtained in [22] : 

a'^fL'{P) + E, (4.99) 

where P = V'^^f]', and 

Af[{P) + \\E\\-p,o < e. (4.100) 
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Together with (4.98), we obtain a decomposition of the following form 

= AiyL'Pi + ^i, (4.101) 

where 

ii^iIl- + iiy%iiL^,L- + iiyL'AiiL^,Ljc +j^[{Pi) + \\e,\\^,o < e. (4.102) 

Using (4.101), (4.102) and the trace estimate (4.87), we deduce 

||A;l'l'IU°?l2 < £, 



which together with (4.97) allows us to improve the estimate for ||5||^oo^2 and ||e||£^oo^2 
given by the bootstrap assumption (4.60). 



4.2.6. Remaining estimates for tr^, and b. In order to improve the re- 

maining estimates in the bootstrap assumptions (4.58)-(4.63), we use the null structure 
equation of section 4.1.2, which consists of transport equations along L and Hodge sys- 
tems on We refer the reader to section 4.8 in [44], where using the L°° bound of trx, 
the trace estimates for £, S and e, and the estimates for the lapse n, we easily obtain the 
remaining estimates. Thus, there exists a universal constant D > such that (4.58)-(4.63) 
hold. This yields (4.40)-(4.45) which concludes the proof of Theorem 4.15. 

4.3. An estimate for LLtry 

In this section, we outhne the main ideas of the proof of Theorem 4.16. Let ^i = h L{ijl) . 
Then, we first derive a transport equation for /ii, and a Hodge system for J/lC- For 
simplicity, we only discuss the transport equation for fii. We differentiate the transport 
equation (4.27) satisfied by ii with respect to L and multiply it by nb. We also use 
commutator formulas of section 4.1.3, the Bianchi identity (4.34) for p, the curvature 
bound (4.39) and the estimates (4.40)-(4.45) obtained in Theorem 4.15. We obtain 



nL(/^i) + ntrx/^i = -2bnf l{C) ■ ftrx - 2bnx ■ [f®fL{C) + b'^fbf l{C) + 2^ lC®C 

+2ntrxbn-yn ■ f l(0 + d^(Fi) + /s, 
where the tensor Fi and the scalar /2 satisfy 

||-^i||l2(h„) + II/2||li(?^„) ^ £■ 

This yields: 

; ( [ {bnx ■ {f^fL{0)dr] +■■■ . (4.103) 

Here, the term 2% comes from the estimate for Fi and /2 together with Bernstein and 
the finite band property for Pj, and the term 22 A(m)£ comes from the initial data term 
for the transport equation - i.e. /ii at t — which is estimated in [42] - together with 



\PM)\\mH.)<2iX{u)e + 2^e + 
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Bernstein for Pj. We have only kept one typical term in the right-hand side of (4.103) for 
the sake of simplicity. 

In view of the desired estimate (4.47), we are bootstrapping an estimate of the type 

\\PjfL(0\\LHH.)<De + D2-h(u) 

for some large enough bootstrap constant D, and where A is a function in L^(]R). Thus, 
estimating directly the term in the right-hand side of (4.103) would yield an upper bound 
of the type 

J2 2'2-^7('Vf \ where 7^^^ e £'(N) and 7^ ^ e £~(N) (4.104) 

l,q 

which is not summable. Instead, we rely on the following decomposition for bnx- 

f{bnx) = fniP + E (4.105) 
where P, E are Pj^^-tangent tensors, and P, E satisfy: 

M(P) + ll^llpo < e. 

Remark 4.28. A similar decomposition has been proved in the geodesic foliation in 
[22], and adapted to the time foliation in the spirit of section 4-2. 4- In order to obtain 
(4.105), we use the fact that the proof in the geodesic foliation relies on a specific structure 
of certain commutators and of the Bianchi identities, which can be recovered in the time 
foliation. We refer to [44] for the details. 

Using the decomposition(4.105), we decompose the term in the right-hand side of 
(4.103) in a sum of two terms which are estimated as follows (see [44] for the details): 

• For the term involving J^ulP, we integrate by parts in y^L, and consider the 
term where the L derivative falls on y lC- Differentiating the transport equation 
(4.21) satisfied by ( with respect to l, commutators formula, and the Bianchi 
identity (4.32), we obtain 

for some elliptic operator of order P° on Pj^^. Wc then integrate by parts the 
L derivatives, and obtain for this term an upper bound of the type 

2Wi(P)PU.(«„) + ---. 

Then, using the estimate for P and the the curvature bound (4.39) for /3, this is 
enough to bound the term involving y„L-P in the right-hand side of (4.103). 

• For the term involving E, we we rely on the Besov improvement for E, and we 
derive an upper bound of the form 

J]2^2-^7W7f\ where 7^^) G e\N) and 7^ ^ G £°°(N), 

l,q 

which is summable unlike (4.104). This is enough to bound the term involving 
E in the right-hand side of (4.103). 
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Remark 4.29. The reader may wonder why the estimate for LLtrx is not better than 
LP' with respect to the variable t - instead of L°° as one should expect since we rely on a 
transport equation for the corresponding quantity Hi. The reason is that boundary terms 
arise from several integration by parts in in the course of the proof. The point is that 
we do not have better estimates than LP with respect to the variable t for these terms. 

4.4. Reguleirity of the foliation with respect to u 

4.4.1. First order derivatives with respect to a;. In this section, we outline the 
main ideas of the proof of Theorem 4.17. Let us first explain how to control d^t^x- 
Differentiating the Raychaudhuri equation (4.19) with respect to a;, we obtain 

L{d^i^x) = [L, dJ\tYx H ■ 

Now, we have 

[L, dJ[tYX = -du:N{tYx) = -^a^jvtrx 
where we used the fact that g(A'", A^) = 1, which differentiated with respect to u implies 
that d^N is a vectorfield tangent to Pt,u- Thus, we obtain 

L{d^trx) = - fd^N^'^X H 

which together with the estimate (4.43) for ^trx immediately yields 

l|5cutrx||L2,L- % £■ 

x' 

Remark 4.30. Ln view of the commutator 

[L, 9J = -fd^N, 

a derivative with respect to ou has essentially the same regularity as a -derivative. 

The estimates in (4.48) and (4.49) are obtained in the same way, i.e. by differentiating 
the ricci equations (4.15) and the transport equations (4.18) (4.19) (4.20) and (4.21) 
with respect to u, computing the commutators [9(^,0^] and [9(^,^1,], and estimating the 
corresponding transport equations (see [44] for the details). 

Next, let us explain how to derive the decomposition (4.50) for £: 

X = Xi + X2, 

where Xi and X2 are two symmetric traceless .^-tangent 2-tensors satisfying the estimates 
(4.51) and (4.52). Recall the Codazzi type equation (4.22) satisfied by 

dj/^ = ^ytrx -/? + •••. 

This is an elliptic system on P^^^, and we may write formally 

x = ^v-ytTx-v-'p+-- - , 
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where P is a pseudodifferential operator of order -1 on Pt^u- This allows us to define xi 
and X2 as 

Xi = \v-^ytix + • • • and X2 = -V'^^. 

The estimate (4.51) corresponds to the estimate (4.44) for x and is prove similarly, so we 
focus on the estimate (4.52). For the sake of clarity, we only explain why, compared to Xi 
Xi has better regularity with respect to {t, x) while X2 has better regularity with respect 
to u, which is the point of the decomposition (4.50) (see Remark 4.18). Indeed, since 
the estimate (4.43) for trx is better than the estimate (4.44) for Xi and since is a 

pseudodifferential operator of order on „, we are able to obtain better regularity in 
{t, x) for xi compared to x- Next, we focus on X2- Now, note that the curvature tensor R 
does not depend on u. Thus, when differentiating (3 with respect to ou, the ou derivative 
falls on the frame (L, L, ca), and we obtain schematically 

d^P^{a + p + a)d^N. 

In particular, we have 

Wdoj^h^ L\H^) < \\d^N\\Lc.{\\a\\Lo.L2^nu) + WpWl^l^Hu) + hh^L^in^)) ^ 

where we used the curvature bound (4.39) for a, p and a, and the estimate (4.48) for di_jN. 
Thus, dijp has the same regularity with respect to (t, x) than In view of the definition 

of X2i we obtain that di^X2 has essentially the same regularity as X2) while the estimate 
(4.49) for di^x looses one ^-derivative with respect to the estimate (4.44) for x- Thus, 
the regularity of X2 with respect to u is better than the corresponding regularity for x- 

4.4.2. Second order derivatives with respect to u. In this section, we outline 
the main ideas of the proof of Theorem 4.19. We focus on the estimate (4.55) for d^C 
which is typical. Differentiating twice with respect to uu the transport equation (4.21) for 
C, and computing the commutator [^/l, d^], we obtain 

fLmolO) = -X ■ n(a^c) + f{F,) + F2 + • • • , (4.106) 

where the Pt,„-tangent tensors Fi and F2 satisfy 

11-^1 II + ||-F2||l1,l2 < e. 

We first get rid of the first term in the right-hand side of (4.106) which is troublesome. 
To this end, we use the following lemma. 

Lemma 4.31. Let 7 denotes the metric induced by g on Pt^u- Let M the Pt^u-tcmgent 
2-tensor defined as the solution of the following transport equation: 

^lMab = MacXcb, Mab = lAB on Po,u, (4.107) 

Then, Mab satisfies the following estimate: 

||M-7||ioo + ||yM||so <£. (4.108) 
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Using the transport equation (4.106) for iK^d'^Q and the transport equation (4.107), 
for M allows us to get rid of the troublesome term x ■ n(9^C)- 

fL{M ■ U{dlO) = f{M ■ Fi) - f{M) . Fi + M • F2 + • • • . 

Together with the finite band property and the Bernstein inequality for Pj, the estimates 
for Fi and F2, and the estimate (4.108) for M, we obtain for M ■ Il{d^() the estimate 
corresponding to (4.55). Then, we obtain the wanted estimate (4.55) for Il{d^() by 
proving that the estimate (4.108) for M is enough to ensure that the multiplication by 
preserves the estimate (4.55). 

4.5. Additional decompositions 

In this section, we outline the main ideas of the proof of Theorem 4.20. We need to 
compare N and trx at two different angles cu and u. The basic tool is the following lemma. 

Lemma 4.32. Let cu and uj' in Let u — u{t,x,uj) and u' — u{t,x,uj'). Then, for 
any tensor F, we have: 

J 1 ( ( /-u+lw-w'l 

\\F\\l^,l\h^,) < \\F\\L^LHHu) + \^-^'\M\F\\l^L^iHu) ( I \\^F\\lHnr)dr 

In order to compare the norms L'^ L"^ {Hu') and L^L^(?^„), we need coordinate sys- 
tems. We define ^t,cj : ^ defined by: 

%,u}{t,x) :— u{t,x,u)u -\- d^u{t,x,u). (4.109) 

Then we claim that ^t,uj is a global diffeomorphism from to and therefore 
provides a global coordinate system on (see Proposition 6.6 for a related result on Eq). 
Next, we prove that 

WnhiH.)^ t [ \F{<l>T,UyWdy'dt. (4.110) 
Jo Jr^ 

This formula allows us to compare the norms L'^L'^(T-Lui) and L'^L'^iT-L^)- In turn, one 
needs to evaluate 

In particular, we need to estimate 

\\d^[KAy')\\\L^- 

We refer the reader to [44] for details on the proof of Lemma 4.32. 

Using Lemma 4.32 as well as commutator estimates for [D, P;] among others, we may 
prove the following corollary. 

Corollary 4.33. Let f a scalar function and uj,uj' in S^. Then, for any I > 0, we 
have: 

\\Pif\\Lo^,LHH^,) < (2-' + \uj- a;'|^2-^)(||/|Ug=i.(^„) + \\BfU^i^.^nJ, 
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and 

\\P<if\\L^,L2(n^,) 
< (1 + la; + la; - ^'|^ 

X (sup 5^ + 

\ " q<l 

We also need the following non sharp commutator lemma. 
Lemma 4.34. Let f a scalar function and uj,uj' in Then, for any I >0, we have: 

\\[du^:P<l]f\\L^,L^H^,) < ||D/||i^i2(^„). 

Using Corollary 4.33 and Lemma 4.34 together with the estimates (4.43) and (4.49) 
for trx and the fact that |a; — i/| < 2~2, we are able to prove the decomposition (4.57) for 
trx- 

Next, we consider N — N^. We have 

2i(7V-7V^)= / d^N{.,J')dw"{2^uj-u)), 

where [a;, i/] denotes the arc of §^ joining uj and v. Since |a; — z/| < 2~2^ we want to 
proceed as for the decomposition of tr^. More precisely, we want to use Corollary 4.33 
and Lemma 4.34 together with the estimates (4.48), (4.49), (4.53), (4.54) and (4.55) for 
dujN , in order to prove the decomposition (4.56) for N — Ny. Now, unlike trx which is a 
scalar, d^jN is a tensor. Since Corollary 4.33 and Lemma 4.34 only apply to scalars, we 
need one last ingredient to prove the decomposition (4.56) for N — and conclude the 
proof of Theorem 4.20. Namely we need to scalarize d^N using a basis of the tangent 
space of which does not depend on u. We refer to [44] for the details. 




CHAPTER 5 



Construction and control of the parametrix at initial time 



In this chapter, and the next one, we will only consider the leave Eq of the foliation 
of A^, and we denote it by S for simplicity. Recall the plane wave type parametrix 
given by (1.26)^ 

"OO 



where u{.,.,uj) is a solution to the eikonal equation ^"^daudj^u = on such that 
u{Q,x,uj) X ■ uj when \x\ +oo on S. The goal of this chapter is to outline the main 
ideas allowing us to obtain the control for that parametrix restricted to E in [45] . 

5.1. Geometric set-up and main results 

5.1.1. Presentation of the parametrix. In this section, we construct a parametrix 
for the following homogeneous wave equation: 

□g0 = on M, /_ 

0|e = 00, T{(l))\^ =01, 

where 0o and 0i are two given functions on E and T is the future oriented unit normal 
to E in the space-time M.. 

We recall the plane wave representation of the solution of the flat wave equation. This 
corresponds to the case where g is the Minkowski metric. (5.1) becomes: 

□0 = on M^+^ 
0(0,.) = 00, 9i0(O,.) = 01 on Ml 

The plane wave representation of the solution of (5.2) is given by: 



(5.2) 



^ (^0o(A.) +^^) .Ad. 



(5.3) 



where F denotes the Fourier transform on R^. 

We would like to construct a parametrix in the curved case similar to (5.3). We 
introduce two solutions u± of the eikonal equation 

Z'^^do^u^dpu^ = on M, (5.4) 

-"^This is actually a half wave parametrix. See (5.6) below for the full parametrix 
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such that: 

T{u±) = t|Vm±| = Tag^ on E, (5.5) 

where V is the covariant derivative on E associated to the metric g induced by g on E, 
I • I is the length associated to g for vectorfields on E, and a± is the lapse of u± onT,. We 
look for a parametrix for (5.1) of the form: 

S{t,x) = SMt,x) + S_f4t,x), (5.6) 

where 

r r+oo 

S±f±{t,x)= / / e'^"±(*'"'")/±(Acu)A2dAda;. (5.7) 
Js^ Jo 

In the next two sections, we specify the parametrix (5.6) by prescribing u± on E and by 
making our choice for f± explicit. 

5.1.1.1. Prescription of u+ and m_ on E. (5.4) and (5.5) are not enough to define u± 
in a unique manner. Indeed, we still need to prescribe ii± on E. To motivate our choice, 
we need to introduce some geometric objects connected to u±. Let N± the vectorfield on 
E defined by: 

iV± = = a±Vii±, (5.8) 

\\/u±\ 

and L± the vectorfield on Ai which is given on E by: 

L± = a^g'^^dau^d^ = a±{-T{u±)T + Vu±) = iT + N±. (5.9) 

Let = {x G E/ u±{x) = u±} denote the level surfaces of u± in E. Since N± is the 
unit normal to P^^, the second fundamental form of P^j. in E is given by: 

0±ie% e±) = g{D^±N±, e±), A,B^1, 2, (5.10) 
where {e^,ef) is an arbitrary orthonormal frame of TP^^. Let 

Uu^ = {{t,x) e M/u±{t,x) = u±} 

denote the null level hypersurfaces of u± in A^. Since L± is null and orthogonal to Pu± 
in T-Lu±, the null second fundamental form x± is given on P^^ by: 

X±{e% e±) = ^(D,±L±, e±), A,B = 1, 2. (5.11) 

Taking the trace in (5.10) and (5.11), and using (5.9) and the fact that k is the second 
fundamental form of E, we obtain: 

trx± = itrA; + tr^±. (5.12) 

Note that tr gk — trk + k^N, where tr^ denotes the trace for 2-tensors on E. In view of the 
maximal fohation assumption (1.6), we have tTgk — 0. Together with (5.12), this yields: 

trx± = l^kN±N± + tr6'±. (5.13) 
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Now, in [44] (see also Theorem 4.15), we prove that trx± belongs to L°°(A^) using a 
transport equation (the Raychaudhuri equation, see (4.19)) provided that it belongs to 
L°°(E) a,t t — 0. Thus, one needs the following estimate 

trx± e L-(E), (5.14) 

which in view of (5.13) is equivalent to: 

^kN^Ni + tre± e L°°(E). (5.15) 

We construct in [42] (see also Chapter 6) a function u{x,u!) on E x such that 

-/c^TV + tr^ e L~(S). (5.16) 

Note that —u{x, —u) satisfies: 

fcATjv + tr^G L°°(S). (5.17) 
Thus, in view of (5.15), (5.16) and (5.17), wc initialize u± on S by: 

M+(0,a;,cj) = u{x,uj) and 1x^(0, x, cj) = — m(x, —oj) for {x,uj) G E x S^. (5.18) 

Remark 5.1. Note that in the particular case where k = - the so-called time sym- 
metric case-, we may take 

M_|_(0, X, cu) = u_(0,x,u) = u{x,u!) for {x,u) G E X S^. 

In particular, we have m+(0, x, uj) — m_(0, x,uj) — x • uj in the flat case. 

5.1.1.2. The choice of f+ and /_. Having defined u±, wc still need to define f± in the 
parametrix (5.6). According to (5.1), the half wave parametrix S+ and S- should satisfy 
on E: 

' s+f+{o,x) + s.f4o,x) = Mx), 

T{S+f+){0,x) + T{S^f^){0,x) ^ M^). 
Let us introduce the following operators acting on functions of M^: 

M±f{x) = e^'^''^'''^^^f{\u)\^dXduj (5.20) 
Js^ Jo 

and 

r r+oo 

Q±f{x)^ / / e^'^''^'''^''^a{x,±uj)-^f{Xuj)X^dXdw, (5.21) 

where a{x,uj) — \'Vu{x,uj)\~^ is the lapse of u. Using (5.5), the definition of S± in (5.6), 
(5.18), the definition (5.20) of M± and the definition (5.21) of Q±, we may rewrite (5.19) 
as: 

' M+/+ + M_/_ = 00, 

g+(A/+)-g_(A/_) = i0i. 

The goal of this chapter will be to show that there exist a unique (/+, /_) satisfying (5.22), 
and that (/+,/_) satisfies the following estimate: 

||A/+||l2(]r3) + ||A/_||l2(]r3) < ||V0o||l2(s) + ||0i||l2(s)- (5.23) 



(5.19) 
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Remark 5.2. In the case of the flat wave equation (5.2), we have (S, g, k) = (M^, 6, 0), 
u±{t,x,u!) — + X ■ uj, u{x,uj) — X • uj and a{x,u!) — 1. In particular, the operators 
M± and Q± defined respectively by (5.20) and (5.21) all coincide with the inverse Fourier 
transform. Then, the system (5.22) admits the following solutions: 



which clearly satisfy the estimate (5.23). 

5.1.2. Geometric set-up. Wc define the lapse a{x,uj) = \'Vu{x,uj)\~^, and the unit 
vector N such that Vu{x,uj) = a{x,u!)^^N{x,u!). We also define the level surfaces P„ ~ 
{x I u{x,uj) — u} so that N is the normal to P^- 

For 1 <p,q < +00, we define the spaces LpL^(P„) using the norm 



We assume that 1/2 < a{x) < 2 for all x G E (see (5.24) below) so that L'p^U'{Pu) coincides 
with LP(S) for all 1 < p < +oo. We denote by 7 the metric induced by g on P^, and by 
y the induced covariant derivative. 

Before stating precisely the main results of this chapter, we first record the regularity 
obtained for the phase u{x,ijj) constructed in [42] (see also Chapter 6). 

5.1.3. Regulsirity eissumptions on the phase u{x,l>j). In this section, we collect 
the estimates for the phase u{x, uj) of our Fourier integral operators that are needed to 
follow the discussion of the control of the parametrix at initial time contained in this 
chapter. An outline of the proof of these estimates will be given in Chapter 6 (see [42] 
for the complete proof). 

We start with the regularity in x of the lapse a. We need: 





\\^a\\L^L\Pu) + ||« - l|U-(s) + ||yVa||L2(s) < e. (5.24) 

We also need a decomposition for Vato. For all j > 0, there are scalar functions a\ and 
03 such that^: 

Vato = a{ + 4 where ||a{||L2(2) < ||a^2lU-L2(p„) < £ 25) 



Next, we consider the regularity with respect to u. We have: 




^we choose a\ = P>j/2iy no) and = P<:ji2(y no), and then obtain (5.25) using (5.24) and an 
estimate for V^a (see [42] for the details) 
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and 

||9>|U,^^(E) < 1. (5.29) 

We will need the following global change of variable on E. Let a; e S^. Let : E — >■ 
defined by: 

(f>oj{x) := u{x,u)u + d^u{x,uj). (5.30) 

Then (f)^^ is a bijection, and the determinant of its Jacobian satisfies the following estimate: 

III det(Jac0c^)| - 1||loo(s) < e. (5.31) 

Finally, wc can compare u{x,uj) with a phase linear in ou. Let u & EP' and (pi, the map 
defined in (5.30). Then, we have: 

u{x,uj) — (j)i,{x) ■ uj = 0{e\uj — z/p), 

d^uix, bj) - d^[(^y[x) ■ u) = 0{6\lj - i/|), (5.32) 
dlu{x,uj) - dliM^) ■ ^) = 0{e). 

Remark 5.3. In (5.24) -(5.32), all inequalities hold for any w G with the constant 
in the right-hand side being independent of uj. Thus, one may take the supremum in cu 
everywhere. To ease the notations, we do not explicitly write down this supremum. 

Remark 5.4. In the case of the fiat wave equation (5.2), we have {T,,g) = {M.^,5), 
u{x,uj) — x-uj,a — 1, N — UJ and — IdR3. Thus, (5.24)-(5.32) are clearly satisfied 
with £ = 0. 

Remark 5.5. Recall that the lapse a is at the level of one derivative ofu with respect 
to X. Thus, we obtain from (5.24) that som,e components ofV^u are in L2(S). Note that 
this is not true for all components since (5.25) does not allow us to control V%a in -f'^(E). 
In fact, (5.25) is consistent with 3/2 derivatives of a with respect to N in L^. 

5.1.4. Main results. We first state a result of boundedness on for Fourier integral 

operators with phase u{x,uj). 

Theorem 5.6. Let u be a function on S x §^ satisfying suitable assumptions (we 
refer to [43] for the complete set of assumptions, and to section 5.1.3 for some typical 
assumptions). Let U the Fourier integral operator with phase u{x,uj) and symbol h{x,uj): 

Uf{x)^ / / e'^''^'''''\x,u)f{\u)\^d\duj. (5.33) 

Let D > 0. We assume furthermore that b{x, uj) satisfies: 

||6||l-(s) + \\Vb\\L^LHP^) + ||yV6|U2(E) < D, (5.34) 

WdML^i:) + l|Va,6|U2(s) < D, (5.35) 

and 

Viv& ^bi + bi where ||6i||L2(E) < ^'^D, H^'slU- l^Cp^) ^ D, 
and ||Viv&^IU2(s) + WbihlL^iP^) ^ 22D. 
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Then, U is bounded on and satisfies the estimate: 

||t//||L2(E)<^||/||L2(M3). (5.37) 

Remark 5.7. We intend to apply Theorem 5.6 to the Fourier integral operators M± 
and Q± introduced in section 5.1.1.2 whose symbol are respectively 1 and a^^ . Thus, our 
assumptions on the regularity of the symbol b{x, u) are consistent with the assumptions 
on the regularity ofa{x,u) given by (5. 24) -(5. 29). 

Remark 5.8. Under the additional assumption (5.32) on u, and under some restric- 
tions on the constant D appearing in (5.34), (5.35), and (5.36), we may prove the opposite 
of (5.37); 

||/IU2(M3) < ||[//|U2(E) 

(see Proposition 5.19). This will he a major ingredient in the proof of Theorem 5.9 below, 
and in particular of (5.41). 

Recall the definition of the Fourier integral operators M± and Q± introduced in section 
5.1.1.2: 

M±f{x) = / / e±*^"("'±")/(Acc;)A'(iA(ia;, (5.38) 



S2 ^0 

and 

P P+OO 

Q±f{x)^ / e^'^''^'''^'^^a(x,±uj)-^f{Xuj)X^dXduj. (5.39) 

The following theorem is the main result of this chapter. 

Theorem 5.9. Let u be a function on S x satisfying suitable assumptions (we 
refer to [43] for the complete set of assumptions, and to section 5.1.3 for some typical 
assumptions). Then, there exist a unique (/+,/_) satisfying: 

M+f+ + M_f_ = <P,, 

g+(A/+)-g_(A/_) = z0i. ^^-^^^ 

Furthermore, (/+, /_) satisfies the following estimate: 

||A/+||l2(r3) + ||A/_||i2(jj3) < ||V0o||l2{e) + ||0i||l2(e)- (5.41) 

Proving the estimates (5.37) and (5.41) for the Fourier integral operators U, M± 
and Q± will require taking several integrations by parts. The main difficulty in proving 
Theorem 5.6 and Theorem 5.9 will be to perform these integrations by parts within the 
very low level of regularity for the phase u{x, u) given by (5.24)-(5.32) and for the symbol 
b{x,u)) given by (5.34) (5.35) (5.36). The proof will rely both on harmonic analysis 
decompositions and the geometry of the foliation of E by m. Theorem 5.6 will be reviewed 
in section 5.2 and Theorem 5.9 will be reviewed in section 5.3. 



5.2. CONTROL OF FOURIER INTEGRAL OPERATORS 101 

5.2. Control of Fourier integral operators 

5.2.1. Structure of the proof of Theorem 5.6. The proof of Theorem 5.6 pro- 
ceeds in three steps. We first localize in frequencies of size A ~ 2-^. We then locahze the 
angle u in patches on the sphere of diameter 2"-'/^. Finally, we estimate the diagonal 
terms. 

Remark 5.10. Note that the structure of the proof is analogous to the one on the 

control of the error term in Chapter 3 (see section 3.2.2). However, the proof each step 
( almost orthogonality in frequency, almost orthogonality in angle, and control of the diag- 
onal term) is different, more particularly the last two steps. 

5.2.1.1. Step 1: decomposition in frequency. For the first step, we introduce (/? and -0 
two smooth compactly supported functions on M such that: 

ip{\) + V(2~^A) = 1 for aU A e R. (5.42) 

We use (5.42) to decompose Uf as follows: 

Uf{x) = U,f{x), (5.43) 

where for j >0: 

P P + OO 

Ujf(x)= / e'^''b(x,uj)^jj(2-^X)f(Xu)X'^dXdu, (5.44) 

is2 Jo 

and 

C/_i/(x)=/ / e'^''b(x,oj)ip(X)f(Xuj)X^dXduj. (5.45) 

J§2 Jo 

The goal of this first step is to prove the following proposition: 

Proposition 5.11. The decomposition (5.43) satisfies an almost orthogonality prop- 
erty: 

l|t^/lli2(E) < E + D'WfWh^^sy (5.46) 

i>-i 

The proof of Proposition 5.11 is postponed to section 5.2.2. 

5.2.1.2. Step 2: decomposition in angle. Proposition 5.11 allows us to estimate 1^2 
instead oi \\U f\\L'^T,)- We perform a second dyadic decomposition. We introduce a smooth 
partition of unity on the sphere S^: 

E Vj{^) = 1 for aU CO e (5.47) 
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where F is a lattice on of size 2~ i , where the support of rjj is a patch on of diameter 
~ 2"-'/^. We use (5.47) to decompose Ujf as follows: 

C/,/(x) = (5.48) 

where: 



U:;f{x)= / e'^%{x,uj)ilj{2-^X)7]'^{uj)f{Xuj)X^dXduj. (5.49) 



(5.50) 



\7 

We also define: 

7-1 = WvWfh^R^), 7i = I|V^(2~^A)/||l2(k3), j > 0, 
7j = ||^(2-U)r;J(a;)/||^.(M3), J>0,i^er, 

which satisfy: 

ii/iii2(M3) = E ^1 = E (5-51) 

The goal of this second step is to prove the following proposition: 

Proposition 5.12. The decomposition (5.48) satisfies an almost orthogonality prop- 
erty: 

\\Ujf\\lH^)<J2\Mf\\lH^) + D'^l (5.52) 
The proof of Proposition 5.12 is postponed to section 5.2.3. 

5.2.1.3. Step 3: control of the diagonal term. Proposition 5.12 allows us to estimate 
||^J/||l2(s) instead of ||f/j/||L2(s)- The diagonal term is estimated as follows. 

Proposition 5.13. The diagonal term U^f satisfies the following estimate: 

||C/;/||l2(e) < D^^. (5.53) 

The proof of Proposition 5.13 is postponed to section 5.2.4. 

5.2.1.4. Proof of Theorem 5.6. Proposition 5.11, 5.12 and 5.13 immediately yield the 
proof of Theorem 5.6. Indeed, (5.46), (5.51), (5.52) and (5.53) imply: 

l|(^/lll-(E) ;£ EllVlll-(E) + B=ll/K=( 

sE'Eii''//iii'(.) + «^E7| + B=ii/iik., 

j>-i udV j>-i (5.54j 

^ E E(^P' + E ^1 + D'WfWhm 

j>-i !/er j>-i 

<^2„^„2 



lL2{ 

which is the conclusion of Theorem 5.6. 

The remainder of section 5.2 is dedicated to the proof of Proposition 5.11, 5.12 and 
5.13. 
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5.2.2. Proof of Proposition 5.11 (almost orthogonality in frequency). We 

have to prove (5.46): 



J>-1 



2 

L2(S) 



This will result from the following inequality using Shur's Lemma: 



Ujf{x)Ukf{x)dE 



< D22-^7,-7fe for \j -k\>2. 



(5.55) 



(5.56) 



We consider a coordinate system {u, x') on E where x' denotes a coordinate system 
on Pu, and we would like to integrate by parts with respect to 9„. Since Vu — a~^N and 
Vu' — a'~^N', we have: 



AXu—iX'u' 



iXu—iX'u' 



(5.57) 



X- X'§g{N, N' 

where we use the notation u for u{x,uj), a for a{x,uj), N for N{x,uj), u' for u{x,u'), a' for 
a{x,u!') and A^' for N{x,lj'). Then, the proof of (5.56) is analogous to the proof of (3.45), 
so we skip it for the sake of simphcity. Let us just say that, as for the proof of (3.45), 
most terms require one integration by parts using (5.57) and the estimate (5.34) for b (as 
in section 3.3.1), while one term requires a second integration by parts using (5.57) and 
the decomposition (5.36) for Viv6 (as in sections 3.3.2 and 3.3.3). 

5.2.3. Proof of Proposition 5.12 (almost orthogonality in angle). We have to 
prove (5.52): 



\Ujf\\] 



< 



(5.58) 



5.2.3.1. Presence of a log-loss. Integrating by parts twice in J^Uj f{x)Uj' f{x)dT, 
would ultimately imply: 



un{x)un{x)dE 



This yields to a log-loss since we have: 



< 



(2i/2|z/-Z/'|)2' 



J- 



\u-u'\^0. 



(5.59) 



(5.60) 



Remark 5.14. Recall that there is an analogous log-loss in the almost orthogonality 
argument in angle for the error term (see section 3.5.2). In section 3.5, we removed the 
log-loss in particular by using integration by parts with respect to the null vectorfield L. 
On the other hand, we work here on S which is Riemannian, so there is no equivalent 
of the null vectorfield L. Instead, we will use a second decomposition in A (see section 
5.2.3.2). Note that such a strategy can not be used to control the error term (see Remark 
5.16). 
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To avoid the log- loss present in (5.59), we will instead derive the following inequality: 

J^u;f{x)Wf^dj: 

where a > 0. Indeed, since is 2 dimensional and 1 < 2^/^|i/ — u'l < 2-'/^ for 
and u ^ 1/', we have: 

and 

sup >^ ^. — < Ca < +00 Va > 0. (5.63) 

2^^/2(2^/2 1 1/ - I/' I) 2-a - " ^ ^ 

Thus, (5.61), (5.62) and (5.63) together with Shur's Lemma imply (5.58). 

5.2.3.2. A second decomposition in frequency. To avoid the log-loss present in (5.59), 
we do a second decomposition in frequency. A belongs to the interval [2-'~^,2^+^] which 
we decompose in intervals Ik'- 

[2J-1, 2^+^] ^ [j Ik where diam(4) ~ 2^> - (5.64) 

i<fe<|i^-i/'|-« 

Let (j)k a partition of unity of the interval [2-'"^, 2^+^] associated to the /^'s. We decompose 
Ujf as follows: 

U-f{x)^ J2 fix), (5.65) 

l<k<\u-v'\-" 

where: 

U';'''f{x)= / e^^"6(a;,a;)^(2-^AX(a;)0fc(A)/(Aa;)A^dAda;. (5.66) 

Remark 5.15. The point of this additional decomposition is to exploit the volume in 
A. Indeed, after performing Cauchy-Schwarz in X, we obtain 

V ~ 22 |i/ — i/'l 2 

which displays the crucial gain |z/ — z/'|t. 

We also define: 

= ||V(2-U)r^;(a;)0fc(A)/|U2(M3), j > 0, u e T,! < k < \u - i/'l"", (5-67) 
which satisfy: 

W"- E ('f)'- (5-68) 
i<ifc<|i/-i/'|-« 



< 



2ja/2^2i/2|z/-I/'| 



+ 



\u-u'\^0, (5.61) 



5.2. CONTROL OF FOURIER INTEGRAL OPERATORS 



105 



5.2.3.3. The two key estimates. We will prove the following two estimates: 

< 



u-;''fix)u;'^''f{x)dE 



+ 



7--i2 i^-k v'k 



2W2(2i/2|zy _i/'|)2-a {2il^\u-u'\f 

for |z/ - z/'l 7^ 0, 1 < A; < |i/ - 



and 



v',k' 



< 



\k - A;'|2i/2(i-4«)(2i/2|z/ - z/'|)i+4«' 
for ^0,l<k,k'< k ^ k'. 



(5.69) and (5.70) imply: 



u-;f{x)uff{x)dE 



(5.69) 



(5.70) 



(5.71) 



< 



< 



J2 I u;''f{x)u;'''f{x)dE + J2 f u;''f{x)u;''''f{x)dj: 

-\„-„i\-a l<kjtk'<\v-v'\-" 



l<k<\u-u'\ 



E 



7-i2 y,k v'k 



i<fe<|i^-i/'|- 



2W2(2J72|i/_i/'|)2-a 



+ 



E 



i<fe<|i/-i/'|-« 



T-vO v,k v' ,k 
(2^/2|i/-I/'|)3 



7-,2 v,k i/',fc' 

^ 7,- Ij 



< 



I^^7j< 



2W2(2i/2|;y _ ;y'|) 



2-a 



+ 



i^^j7r 



(2j/2|zy _zy'|)3' 

where we have used (5.68) and the fact that we may choose < a < 1/5, together with 
the fact that: 

Since (5.71) yields the wanted estimate (5.61), we are left with proving (5.69) and (5.70). 
The discussion in the following section will be very informal for the sake of simplicity. We 
refer to [43] for the details. 

5.2.3.4. Proof of (5.69). The estimate (5.69) will result of two integrations by parts 
with respect to tangential derivatives (in the spirit of section 3.5.1.1). Let us consider 
for instance the case where the two tangential derivatives fall on the symbol h of Up^f 
defined in (5.66). This yields a term of the form 



^2 ^ 



e^^X2-^A)(/)jt(A)/(Aa;)A^ciA ) rfAu)du. 



(5.73) 



Then, in view of the estimate (5.34) for 6, we have in particular JP'h G i>^(E) which will 
force us to estimate the A integral in (5.73) in L^. To this end, we do Cauchy-Schwarz 
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and obtain in particular the square root of the diameter of Ik- Due to (5.64), we thus 
gain an additional factor of \v — v'\°' with respect to (5.59), which yields (5.69). 

Remark 5.16. Here, the log-loss is removed by exploiting the size of the diameter 
of Ik- This is possible since we estimate the A integral of (5.73) in using Cauchy 

Schwartz. In turn, this is a consequence of our estimate for J/^b in L^(E). This explains 
why this method cannot be used to remove the log-loss of the error term in Chapter 3. 
Indeed, our estimates for the space-time foliation in Chapter 4 are typically of the type 
L'^L^{T-lu)) so that the integral in A is estimated in L\ using Plancherel. In turn, this 
does not allow us to see the size of the localization in \, so that a second decomposition 
in frequency of the type (5.65) would be useless in that case. 

5.2.3.5. Proof of (5.70). Note that we not only need to gain summability in (z/, u') for 
this term, but also in (A;, k'). This is achieved through the presence of the additional gain 
k — k' m. the right-hand side of (5.70). The estimate (5.70) will result of two integrations 
by parts, one with respect to the normal derivative N , and one with respect to tangential 
derivatives. We obtain a term analogous to (5.73) 

"+00 \ 

e'^X2-^A)0fc(A)/(Aa;)A2(iA rfAuj)duj. (5.74) 



/ y^Nblf 



In view of the estimate (5.34) for 6, we have in particular yv^b G i^^(S). Thus, the 
log-loss of the summation in (z/, z/') is removed as in the previous section, in particular 
using the size of the diameter of Ik. Note also that the gain k — k' in the right-hand side 
of (5.70) comes from the integration by parts in N. Indeed, we use 



A\u-i\'u' _ Y7 („iXu-iX'u'\ /c 7c^ 

--X-X'^g{N,N')^-^' ^'-^'^ 

Now, since A e 7^, A e Ik', we have in view of (5.64), the assumption (5.27) for d^a, and 
the assumption (5.28) for d^jN 



a' 



\k-k'\2^\i 

which yields the gain k — k' in the right-hand side of (5.70). 



z/T, 



5.2.4. Proof of Proposition 5.13 (control of the diagonal term). We have to 
prove (5.53): 

< D^^. (5.76) 

Recall that Uj is given by: 

U^fix) = / bFj{u)r)f{uj)dw, (5.77) 

where Fj{u) is defined by: 

r+00 

Fj{u) = / e*^X2-^A)/(Ac^)A^ciA. (5.78) 
Jo 
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We decompose Uj in the sum of two terms: 

mf{x)^b{x,u) I Fj{u)r]Uuj)duj + I {b{x,uj)-b{x,iy))Fj{u)r]Uu;)du;. 

Then, using in particular the assumption (5.34) for b and the assumption (5.35) for di^b, 
we obtain 



L2(E) 



< D 



L2(E) 



Fj{u)vi''j{uj)du 

In order to estimate the right-hand side of (5.79), we use the following proposition. 
Proposition 5.17. We have the following bound: 



(5.79) 



Fj{u)r]'^{uj)duj 



(5.80) 



The proof of Proposition 5.17 is postponed to the next section. Finally, (5.79) and 
(5.80) yield the wanted estimate (5.76) which concludes the proof of Proposition 5.13 
5.2.4.1. Proof of Proposition 5.11. Recall that Jga -^j(ii)?7j(^)c^^ is given by: 

■" + 00 



§2 ^0 



e'^''iP{2~^X)r]'^{uj)f{Xuj)X'dXduj. 



(5.81) 



Relying on the classical TT* argument, (5.80) is equivalent to proving the boundedness 
on -L^(E) of the operator whose kernel K is given by: 



K{x,y) 



|2 ^0 



eiMx,^)-iMy,^)^^2-^XX{uj)X''dXdw, x,ye^. 



(5.82) 



The decay satisfied by this kernel is stated in the following proposition. 



Proposition 5.18. The kernel K defined in (5.82) satisfies the following decay esti- 
mate for all x,y in E; 



\K{x,: 



< 



O-'l 



{l + \2J\u{x,u)-u{y,u)\-2J/^\d^u{x,u)-d^u{y,u)\\)^ 
2^ 



(5.83) 



The proof of Proposition 5.18 is postponed to section 5.2.4.2. In the rest of this section, 
we show how (5.83) implies Proposition 5.17. According to Schur's Lemma, the operator 
whose kernel is K is bounded on I/^(S) provided we can prove the following bound: 

sup / \K{x,y)\dy < +00, sup / \K{x,y)\dx < +00. (5.84) 
ses JE yes Js 

Due to the symmetry of K in x,y, the two bounds in (5.84) are obtained in the same 
way. We focus on estabhshing the first bound. Using in particular (5.83) and the global 
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change of variable on E given by (5.30)^, we are able to obtain: 



L 



where y = y -u + y' and y' -u = 0. Making the change of variable y ^ z where z is defined 
hy z ■ u = 2^y ■ u and z' = 2^1'^y' in the right-hand side of (5.85), and remarking that z ■ v 
is one dimensional, and z' is two dimensional, we obtain: 

f dz 

(5.86) implies the first bound in (5.84). K being symmetric with respect to x,y, the 
second bound in (5.84) is also true. Thus, the operator whose kernel is K is bounded on 
L^(E) which concludes the proof of Proposition 5.17. 

5.2.4.2. Proof of Proposition 5.18. Recall the definition of K: 

K(x,y)= / e^^"(^''^)-^^"(^''^V(2~^'A)r7j(cc;)A2dAda;, X,?/ e E. (5.87) 

is2 Jq 



We need to prove that K satisfies the following decay estimate for all x, y in E: 

2^ 



\K{x,y)\ 



^ {l + \2j\u{x,iy)-u{y,u)\-2j/'^\d^u{x,u)-d^u{y,u)\\)'' 

X- 



(5.88) 



i + 2^/^\dMx, v)-dMy.^W 



For the sake of simplicity, let us just describe the general strategy of the proof of Propo- 
sition 5.18. In view of the regularity for ii(a;, cj) with respect to ui provided by (5.24) and 
(5. 26)- (5. 29), we have 

\u{x,uj)\ + \d^u{x,uj)\ + \dlu{x,uj)\ + \dlu{x,uj)\ < 1 + Vx e E, Vu; e (5.89) 

This regularity allows us to integrate by part three times with respect to a;, while we may 
integrate as much as we want with respect to A. The estimate (5.88) is then obtained after 
performing in (5.87) three integrations by parts with respect to ui and two integrations 
by parts with respect to A. 



5.3. Control of the parametrix at initial time 

In this section, we discuss the proof of Theorem 5.9. To this end, we first show that 
the Fourier integral operator U of Theorem 5.6 almost preserves the norm provided 
we make additional assumptions on its symbol. We then use this observation to prove 
the estimate (5.41). Finally, we conclude the proof of Theorem 5.9 by establishing the 
existence and uniqueness of (/+, /_) solution of the system (5.40). 



'using also the bound on the Jacobian (5.31) 
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5.3.1. A refinement of Theorem 5.6. In Theorem 5.6, wc have proved that the 
Fourier integral operator U with phase u and symbol b is bounded on provided u 

satisfies (5.24), (5.26)-(5.29) and (5.30)-(5.31), and the symbol b satisfies (5.34) (5.35). 
We now would like to prove that U satisfies the following bound from below: 



|l^(M3) < ||t//||L^(E), (5.90) 

provided u also satisfies (5.32) and under additional assumptions on the symbol b. This 
is the aim of the following proposition. 

Proposition 5.19. Let u he a function on E x satisfying suitable assumptions (we 
refer to [43] for the complete set of assumptions, and to section 5.1.3 for some typical 
assumptions). Let U the Fourier integral operator with phase u{x,u!) and symbol b{x,(x!): 

Uf{x)^ I r e'^''^''''^^b{x,uj)f{Xuj)X^dXduj. (5.91) 

We assume furthermore that b{x, u) satisfies: 

||5.&||l2(s) + ||Va^6|U2(s) < 1, (5.92) 
\\b - Ih^iM) + ||V6|U2=i2(p^) + ||yV6|U2(s) < e, (5.93) 

and 

Vjv6 = b{ + bi where \\b{\\L2(S) < 2"2£, \\bi\\L^L^p„) < £ 
and ||V7v6^||l2(s) + ||&2IIl2l->(p„) < 22£. 
Then, U is bounded on L^ and satisfies the estimate: 

U2(M3) < ||C//|U2(E). (5.95) 



(5.94) 



Remark 5.20. Notice that the only difference in the assumptions with respect to The- 
orem 5.6 lies in the fact that u also satisfies (5.32) and in the constant D which has been 
replaced by 1 in (5.92) and by e in (5.93) (5.94). 

The proof of Proposition 5.19 uses the decomposition in frequency and angle of the 
operator U introduced in section 5.2. In order to control the diagonal term in a third step 
(see next section), we have to modify slightly the size of the support of our partition of 
unity rfj on introduced in (5.47). Let 5 > such that: 

< « « 1- (5.96) 

We now require that the support of r]f is a patch on of diameter ~ b2~^l'^ . With this 
modification, the assumptions for b in Proposition 5.19, and by carefully tracking the size 
of the various terms in the almost orthogonality argument in frequency and angle, and 
the control of the diagonal term, we obtain 

= E E / S]f{x)SfJ{x)dT. + O + 5) ||/||i.(«3), (5.97) 
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where the operator is defined on S by: 

S^f{x)^ / / e^^"(^''"V(2"'A)r7j(a;)/(Aa;)A'dAda;. (5.98) 
is Jo 

Recall (5.30)-(5.31) which states that the map : E ^ defined by: 

(puix) :— u{x , u)!/ + di^u{x , ly) , (5.99) 

is a bijection, such that the determinant of its Jacobian satisfies the following estimate: 

III det(Jac0i,)| - 1||l->(x) < £• (5.100) 

Let us note the inverse Fourier transform on R^. We introduce the operator Sj on 
E defined by: 

S^f{x)^J'-\^Pi2-^-Xf)iM^))^ [ e''^-^^'>^,l;i2-^XXiu;)f{Xu)X'dXdw. (5.101) 
The following proposition shows that the term J^Sjf{x)Si'f{x)dT, is close to the term 

J^s^f{x)sr'f{x)dj:. 

PROPOSITION 5.21. We have the following bound: 

\\S^f - S^fh^i:) < 5H- (5.102) 

The proof of Proposition 5.21 relics on the classical TT* argument, and the comparison 
between u{x,uj) and 0i/(x) • uj provided by (5.32) (see [43] for the details). Now, (5.97) 
and (5.102) yield: 

wufwW) = E I s^mWmd^ +0(^+5^) ii/iii.(^3). (5.103) 

\j-l\<2 \i,-u'\<2S2-^/^ ^ 



Making the change of variable y = (f)u{x) in J^Sj f{x)Sf f{x)dT, and using (5.100) and 
(5.101) implies: 



EE/ s^msr'nx)d^ 



\j-l\<2 |i/-jy'|<2(52-j/2 



= E E j /-'m-^-Xf){y)F-'{^{2-^-X'f){y)dy 

\j-l\<2 \v-v'\<2&2-il^ -^^^ C5 ;l04) 

+0(^)||/||i.(M3) ^ ■ ^ 



\ V'(2-^A)r/J(a;)/(Aa;)^(2-'AX('^)/(Aa;)% 



\3-l\<2 |i/-i/'|<2(52-j/2 

+o(£)||/||i.(M3), 
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where we have used the fact that T ^ is an isomorphism on L^(M^) in the last equahty 
of (5.104). Now, we have: 

J2 J2 f V'(2-^A)r;J(^)/(A^)V^(2-'AX(a;)/(Aa;)d|/ = ||/||i.(M3), (5-105) 
which together with (5.103) and (5.104) yields: 

l|t^/lli^(E) = ll/lli^(M3) + o + ll/lli^(M3)- (5.106) 

Choosing 5^ and eS~'^ small enough, we deduce from (5.106): 

||/||l^(m3) < ||C//||l^(e), (5.107) 
which is the wanted estimate. This conclude the proof of Proposition 5.19. 

5.3.2. Proof of Theorem 5.9. The symbol of the Fourier integral operators M± and 
Q± are respectively given by 1 and a{x, ±00)"^. Thus, they clearly satisfy the assumptions 
of Proposition 5.19. Relying on Proposition 5.19, we are then able to prove the estimate 
(5.41). We refer to [43] for the details. 

The uniqueness of (/+, /_) solution of the system (5.40) is an immediate consequence of 
the estimate (5.41), so there only remains to prove the existence of (/+, /_) to conclude the 
proof of Theorem 5.9. Recall that the phase u{x,u) of our Fourier integral operators has 
been constructed in [42] (see also Chapter 6) on S x §^ under the assumption that (S, g, k) 
satisfies the following bounds consistent with the bounded curvature conjecture: 

<e, ||VA;||l2(s) <e, (5.108) 

where the fact that we may take e > small comes from a reduction to the small data 
case. (S, g, k) also satisfies the constraint equations and the maximal foliation assumption 

V^fcij = 0, 

R=\k\^, (5.109) 
TrA; = 0. 

We introduce two sets V and W: 

V = {(E, g, k) such that (5.108) and (5.109) are satisfied}, (5.110) 

and 

1^ = { (E, gi, /c) e y such that (/+, /_) solution of (5.40) exist for all (0o, <t>i) 
such that V(/)o G L^(E) and 0i e -L^(E)}. 

(5.111) 

Not first that W is not empty since (S, g, k) = (M^, 6, 0) belongs to W in view of Remark 
5.2. We then show that V is connected and W is both open and closed in V for a suitable 
topology (see [43] for the details). We infer W = V. This proves the existence of (/+, /_) 
solution of (5.40) and concludes the proof of Theorem 5.9. 



CHAPTER 6 



Control of the foliation at initial time 

In this chapter, we will only consider u{t,x,uj) and at t = 0. Thus, for simplicity, 
we denote in the rest of this chapter u{0,x,uj) by u{x,uj) and Eq by S. The goal of this 
chapter is to prove the estimates on the control of the foliation of E by u{x,u) which 
are needed for the proof of Theorem 5.9 (see section 5.1.3), i.e. for the control of the 
paramctrix at initial time. The estimates obtained for u{x,u) in this chapter must also 
be consistent with the control on J\4 for u{t,x,(x!) obtained in Chapter 4 (see section 
4.1.5). Here, we outhne the main ideas and we refer to [42] for the details. 

6.1. Geometric set-up and main results 

6.1.1. Reduction to small data. Recall from section 1.2.3 that we have reduced 
ourselves to an asymptotically flat initial data set (E, g, k) solution to the constraint 
equations which satisfies the bounds 

||-R|U2{s) < £, ||VA;||l2(s) < £, (6.1) 

and is smooth outside of a small neighborhood U. In order to construct u{x, u) satisfying 
the asymptotic behavior u{x, uj) ~ x - a; when \x\ — >■ +oo on E, we need to modify (E, g, k) 
outside of U. We can glue it in a trivial way to (M^, S, 0) so that the new initial data set is 
still smooth outside of U, satisfies (6.1), and coincides with (M^,5, 0) outside of a slightly 
larger neighborhood. We still denote this initial data set (E,gf,A;). Of course, (Ti,g,k) 
does not satisfies the constraint equations in the annulus where the gluing takes place. 
However, for the construction of u{x, cu), we only require (E, g, k) to satisfy the constraint 
equations in U. Outside of U, (E, g, k) is smooth, so things are much easier. 

6.1.2. Geometry of the foliation of S by a scalar function u. We define the 
lapse a — |Vit|~^, and the unit vector N such that Vu — a~^N. We also define the level 
surfaces P„,j = {x / u — Uq} so that N is the normal to P^- The second fundamental form 
9 of Pu is defined by 

e{X,Y) = g{VxN,Y) (6.2) 

with X, Y arbitrary vectorfields tangent to the -u-foliation P.^ of E and where V denotes 
the covariant differentiation with respect to g. We denote by tr6' the trace of 9, i.e. 
tr^ = S^^9ab where 9ab are the components of 9 relative to an orthonormal frame 
{eA)A=i,2 on P„. 

113 



114 



6. CONTROL OF THE FOLIATION AT INITIAL TIME 



6.1.3. Structure equations of the foliation of S by a scalar function u. We 

recall some of the structure equations of the foliation of E by a scalar function u which 
will be needed for the discussion of the main result of this chapter (see [9] for a proof). 

PROPOSITION 6.1. The orthonormal frame frame N,eA,A — 1,2 of T, satisfies the 
following system: 

VjvAT = -a-ya. ^ ' ' 

Also, the lapse a satisfies the following equation 

a-^^(a) = -VNtr^-\0\^ + Rnn, (6.4) 

where ^ is the Laplace- Beltrami for the metric 7 on P„ induced by g. Finally, the second 
fundamental form 9 satisfies the following Codazzi equation 

f'OAB = fAtr9 + Rna- (6.5) 

6.1.4. Choice of u{x,(x!). We look for u{x,(x!) satisfying the three following condi- 
tions: 

(a) u{x,uj) ~ x.uj when \x\ — )■ +00 on E 

(b) The regularity ofu{x,uj) with respect to x and ou is consistent with the regularity 
of u{t,x,uj) with respect to (t,x) and u obtained in Chapter 4 (see section 4.1.5). 
In particular, we have trd — k^N G L°^(S) (see the discussion in section 5.1.1.1) 

(c) M(a;, uj) has as enough regularity in x and uj to control the parametrix at initial 
time, i.e. to obtain the conclusion of Theorem 5.9 

where the initial data set (E, g, k) satisfies 

v^hj = 0, 

i?=|fc|2, (6.6) 

tTgk — 0, 

in U (see section 6.1.1), and where R and VA; are in -L^(E) and satisfy the smallness 
assumption (6.1). 

In order to motivate our choice of u{x,lj), we investigate the regularity of the lapse a, 
which by (6.4) satisfies the following equation: 

a-^f (a) = -V^tr^ - - Rnn- (6.7) 

Since R is in ^^(S), (6.7) implies that a has at most two derivatives in ^^(S). Thus, 
u{x,uj) has at most three derivatives with respect to x in iv^(E). This is not enough to 
satisfy (c). In fact, the classical T*T argument (see for example [39]) relies on integrations 
by parts in x and would require at least one more derivative since E has dimension 3. 

Alternatively, we could try to use the TT* argument which relies on integrations by 
parts in uu. Indeed, R being independent of uj, one would expect the regularity of u{x,u) 
with respect to ou to be better. Differentiating (6.7) with respect to ou, we obtain: 

a-^f (9^a) = 2fV Na + ■■■ , (6.8) 
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where the term on the right-hand side comes from the commutator [d^^, ^] (see section 6.4). 
Thus, obtaining an estimate for d^^a from (6.8) requires to control VatO- Unfortunately, 
(6.7) seems to give control of tangential derivatives of a only. This is where the specific 
choice of u^x^oj) comes into play. 

Having in mind the equation of minimal surfaces (i.e. tr6' = 0), condition (b) suggest 
the choice tr^ — kj^^ = 0. Unfortunately, this equation together with (6.7) does not 
provide any control of V/va. We might propose as a second natural guess to take instead 
tr^ — /cjvAT = VatC. Plugging in (6.7) yields an elliptic equation for a: V^a + a~^/^{a) ~ 
— I^P — Vn^nn — Rnn- This allows us to control V%a in I/^(E). However, tr6' — k^N — 
Vjvo, and Viva is at most in i^^(S) which does not embed in L°°(E) - since E has 
dimension 3 - so that condition (b) is not satisfied. To sum up, the first guess trO—kNN — 
satisfies (b), but not (c), whereas the second guess tr^ — kNN = Vato might satisfy (c), 
but does not satisfy (b). 

The correct choice is the intermediate one 

tr^ — kNN = I — a. (6.9) 

We will see in section 6.2 that a — 1 belongs to L°°(E) so that (b) is satisfied. Also, 
plugging (6.9) in (6.7) yields the parabolic equation: 

V^a - a"^^(a) = j^j^ + VjvA;jvjv + Rnn- (6.10) 

This will allow us to control normal derivatives of a. In turn, we will control derivatives of 
a with respect to cu using (6.8). Ultimately, we will prove enough regularity with respect 
to both X and cu, such that (c) is satisfied. 

6.1.5. Main results. For I < p,q < +00, we define the spaces L^L'^^Pu) for tensors 
F on E using the norm: 

i/p 



\\F\\lIli{p^) = \\F\\l,^p^)du 

Remark 6.2. In the rest of the paper, all inequalities hold for any a; e with the 
constant in the right-hand side being independent of ui. Thus, one may take the supremum 
in CO everywhere. To ease the notations, we do not explicitly write down this supremum. 

We first state a result of existence and regularity with respect to x ior u. 

Theorem 6.3. Let {'Ljgjk) chosen as in section 6.1.1. There exists a scalar function 
u on'E X Ei^ satisfying assumption (a) and such that: 

\\a - l|Ug°L2(p„) + ||Va||L-L2(p„) + \\a - 1||l-(s) + Wf^aWi^j:) < e, 
\\tr9 - kNN\\Loo(j:) + \\W9\\l2(j:) < e, 

where Pu, a, N and 6 are associated to u as in section 6.1.2. 

Notice that condition (b) is included in (6.11). In order to state our second result, we 
introduce fractional Sobolev spaces H^{Pu) on the surfaces for any 6 e R (see [42] for 
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their precise definition). We have the following estimate for V%a, and improved estimate 
for Vato. 

Theorem 6.4. Let {T,,g,k) chosen as in section 6.1.1. Let u the scalar function on 
E X constructed in theorem 6.3, and let a and N he associated to u as in section 
6.1.2. We have: 

||Viva|U^i4(p„) + l|V^a||^,^_i^^^^ < e. (6.12) 

The third theorem investigates the regularity of u with respect to u: 

Theorem 6.5. Let {T^,g,k) chosen as in section 6.1.1. Let u the scalar function on 
E X constructed in theorem 6.3, and let P^, a, N and 6 he associated to u as in section 
6.1.2. We have: 



||9a,a||L-(E) + \\fdu;a\\L^L^Pu) + \\f d^a\\L^{i:) + l|Vjv9c^a|| 1 



+l|v^a,a||^^^_3^^^^ + ||V9,^|U2(S) < £, ||a,iv|U-(E) < i, 



(6.13) 



II^^AriU.(^) < 1, ^ • ^ 
and 

\\dlu\\Lr,^ii:) < 1. (6.15) 

6.1.6. Additional results. The following proposition establishes the existence of a 
global coordinate system on E. 

Proposition 6.6. Let a; e Let (f)^^ : T, ^ defined hy: 

(l)i^{x) :— u{x,uj)uj + di^u{x,uj). (6.16) 

Then (f)^ is a hijection, and the determinant of its Jacohian satisfies the following estimate: 

III det( Jac0^) I - l|U->(s) < (6.17) 

Below, we state several additional estimates. We start with a first proposition. 

Proposition 6.7. Let (T,,g,k) chosen as in section 6.1.1. Let u the scalar function 
on S X constructed in theorem 6.3. Let £ and (j)^ the map defined in (6.16). Then, 
we have for all x & and G §^ 

u{x, to) — 4>u{x) ■ CO = 0{e\co — z/|^), 

d^u{x, u) - d^{(j)^{x) ■ co) = 0{e\co (6.18) 
dlu{x,uj) - dliMx) ■ uj) = 0{e). 

Using the geometric Littlewood Paley projections Pj on P„ constructed in [24] (see 
section 3.1.5) together with the estimates for Vjvot in (6.11), and the estimate for V%a in 
(6.12), we obtain the following proposition: 
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Proposition 6.8. Let {T,,g,k) chosen as in section 6.1.1. Let u the scalar function 
on E X constructed in theorem 6.3, and let a and N he associated to u as in section 
6.1.2. For all j > 0, there are scalar functions a{ and such that: 

VatO = a{ + a{ where ||a{||L2(s) ^ 2^ie and WV NdiWL'^^T.) ^ S*^. (6.19) 

Remark 6.9. Recall from section 6.1. 4 that we do not have enough regularity in x 
to apply the T*T method. Alternatively, we could try the TT* method which relies on 
integration by parts in cu. But d^u e L°°(E) is also not enough and we would need at least 
one more derivative in u. Nevertheless, using the regularity in x and u obtained for u in 
the present and the previous section, we are able to control the parametrix at initial time 
(see Chapter 5). 

Let us conclude this section by mentioning several ingredients of [42] that have been 
omitted here for the sake of simplicity, and that have to be proved by relying on low 
regularity assumptions for u which are consistent with the results stated in this section 
and the previous one: 

• estimates for the parabolic operator (Vat — 

• estimates for 9 and N 

• product estimates in the Sobolev spaces H^[Py) 

• embeddings on E and Pu 

• a control of the Gauss curvature of P„ 

• Bochner inequalities on P„ 

• estimates for various commutator terms of the type: [Vat, f], [Vat, Pj], ... 

The rest of this chapter is as follows. In section 6.2, we discuss the proof of Theorem 
6.3. In section 6.3, we discuss the proof of Theorem 6.4. In section 6.4, we discuss 
the proof of Theorem 6.5. In section 6.5, we discuss the proof of Proposition 6.13 and 
Proposition 6.6. Finally, Proposition 6.7 and Proposition 6.8 are discussed in section 6.6. 

6.2. Construction of the foliation and regularity with respect to x 

In this section, we discuss the proof of Theorem 6.3. By section 6.1.1, we may assume 
that {T,,g,k) coincides with (R^,5, 0) outside of a compact, say \x\ > 1. Notice that in 

> 1 and for all uj G S^, the scalar function x.u satisfies the equation (6.9) and the 
estimate (6.11), since a = 1,6 = and N = u in this region. Thus, we would like to 
construct a function u solution of (6.9) satisfying (6.11) in a region containing \x\ < 2 
and to glue it to x.u in 1 < \x\ < 2. Now, (6.9) is of parabohc type - see (6.10) - where 
u plays the role of time. For each fixed oj e S^, we start with u = —2 on x ■ oj — —2. 
Then, we propagate with the parabolic equation (6.10), coupled with the equation for 6 
(6.5), to the strip S = {x E T, such that — 2 < u{x,u)) < 2}. This strip covers the entire 
region < 1, and we then glue u to x ■ cu outside of \x\ < 1 (see section 6.2.2). In the 
next section, we prove a priori estimates consistent with the estimate (6.11) and valid on 
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—2 < u < 2 for the solution u of: 

tr^ - kNN = 1 - a, on - 2<u<2, (6.20) 

where u is initiahzed on x.uj = —2 by: 

u{x,uj) = —2 on x.u = —2. (6.21) 

Note that the first equation of (6.25), (6.21) and the fact that {g,k,Yl) coincides with 
(5,0,M2) for \x\ > 2 yields: 

W{a - 1) = 0, V^*^ = 0, V^(A^ - a;) = for all p e N on M = -2. (6.22) 

6.2.1. A priori estimates for lower order derivatives. Let {Ti,g,k) chosen as in 
section 6.1.1. In particular, we assume: 

||VA;|U2(s) + \\R\\mi:) < e. (6.23) 

Let u a scalar function on E x and let P„, a, and 9 be associated to u as in section 
6.1.2. Assume that u satisfies the additional equation (6.9), which we recall below together 
with (6.3) and (6.4): 

VArA = -yiog(a), ^ 
and 

tr6' — /cattv = 1 — a, 



Vjva - a-^^(a) = |6'p + VNkNN + Rnn- ^^'^^^ 
In this section, wc establish a priori estimates for a, N and 6 corresponding to (6.11) in 
the region 5* of S between P_2 and P2 (i.e. S = {x / — 2 < u{x, cu) < 2}) where u is 
initialized on x.cu — —2 by (6.21). In particular, we have (6.22), so that the subsequent 
integrations by parts will not create boundary terms at u = —2. 

For the sake of simplicity, let us just discuss the estimate (6.11) for the lapse a. We 
rewrite the second equation of (6.25) as: 

(Vjv-a-^f)(a-l) = /i, (6.26) 

where h is given by: 

h = VNkNN + RNN + -- - ■ (6.27) 
Using in particular (6.23), we obtain: 



Together with (6.26) and an parabolic estimate for the operator (Vjv — a ^^), we 
obtain 

||a - l|U-L2(p„) + \\fa\\L^L2(Pu) + ||Viva||L2(s) + ||y^a||i2(2) < e. (6.28) 
In order to obtain estimates for ^VAra and V%a, we differentiate the second equation 
of (6.25) by Vjv: 

(Viv - a-^f ) Viva = Vj^kNN + VnRnn + ■■■ , (6.29) 
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where we only kept two typical terms. Note that V%kNN and VnRnn are dangerous 
terms which cannot be estimated directly. We first need to trade a Vat derivative with 
a y derivative using Bianchi identities and the constraints equations. We use the twice- 
contracted Bianchi identity on E 

V^Rij = ^ViR. (6.30) 

In particular, using also the constraint equations (6.6), we have 

VnRnn = —^aRan + Viv|A;|^. 
Also, the constraint equations (6.6) yield 

V^fciVAf = -VN'^AkAN -\ -VAVivA:AAr H ■ 

Together with (6.29), we obtain 

(Vat - a-^^)VNa = d]^{H) + hi, (6.31) 

where 

H^-VNk.N-R.N, (6.32) 

and hi satisfies 

llft.ll,.,|,P.,<^. (6.33) 

Using the smallness assumption (6.23) and the definition of H (6.32), we have 

II^^IU^(E) < e (6.34) 

which together with (6.33), (6.31), and an parabolic estimate for the operator (Vat — 
a~^^), yields 

||VAfa||LS°L2(P„) + ||yVAra||L2(s) < £. (6.35) 

Finally, (6.28) and (6.35) yield the wanted estimate (6.11). 

6.2.2. End of the proof of Theorem 6.3. We briefly sketch the rest of the proof 
of Theorem 6.3, and we refer to [42] for the details. In (6.28) and (6.35), we have obtained 
a priori estimates consistent with the estimate (6.11) and valid on — 2 < < 2 for the 
solution u of (6.20). Then, we also prove on — 2 < < 2 a priori estimates for higher 
derivatives of the solution u oi (6.20). We then use the existence of u solution to^: 

tr^ — k^N — I — a, ona<u<a + T, (a Qfi^ 

u — a on u — a, \ ■ ) 

where —2 < a < 2, m is smooth, and T > is small enough. Together with the a priori 
estimates, this allows us to control the solution of (6.36) on —2 + kT < u < —2 + (A;+ 1)T 
uniformly with respect to /c = 0, . . . , [4/T] in order to obtain a solution u of (6.20) on 
—2<u<2. Finally, we conclude the proof of Theorem 6.3 by showing how to glue the 



-""this local existence result could be proved either using a Nash Moser procedure or a combination of 
Cauchy-Kowalewska and enhanced a priori estimates for all derivatives 
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solution u of (6.20) to x.u in 1 < |a;| < 2 in order to obtain a solution on S satisfying 
(6.11). 

6.3. Estimates for Vatc and V%a 

In this section, we discuss the proof of Theorem 6.4. Recall the decomposition (6.31), 
(6.32), and the estimate (6.33). We introduce the scalar functions on S oi and 02 solutions 
of: 

(Vat - a'^^)ai = hi on S, ai(-2, .) = 0, (6.37) 

and: 

(Vjv - a-^f )a2 = d]^{H) on S, a2(-2, .) = 0, (6.38) 
which yields, in view of (6.31) and (6.22), the decomposition: 

Vato^ 01 + 02- (6.39) 

The idea behind the decomposition (6.39) is to take advantage of the better regularity 
of hi for ai (see (6.33) compared to (6.34)), and to use the structure of d]/^{H) to obtain 
a useful equation for Indeed, in view of the equation (6.38) satisfied by 02, 

satisfies: 

(Vat - a-i^)(V^a2) = VA.(di^(i/)) + • • • , (6.40) 

and using in particular the twice-contracted Bianchi identity on E, the constraint equa- 
tions in the maximal foliation (6.6), and (6.40), we obtain 

(Vat - a-^^)VNa2 = d\^di^{Hi) + d\^{H2) + /i2 + • • • , (6.41) 
where the tensors Hi, H2 and the scalar h2 satisfy 

+ ll^2||^2i|(p„) + II^2||l1(S) < 

These ideas allow us to derive the following two propositions (see [42] for the detailed 
proof of these propositions). 

PROPOSITION 6.10. Let ai be the solution of (6.37), where hi satisfies (6.33). Then, 
we have: 

||ai|U-L4(p^) < £, (6.42) 

and: 

J]2-^■||P,(VA^aO||i.(s)<£^ (6.43) 

j>Q 

Proposition 6.11. Let 02 be the solution of (6.38), where H is defined in (6.32). 
Then, we have: 

||a2||L-L4(p„) < (6.44) 

and: 

sup||P,(V7va2)|U2(s) <£■ (6.45) 

j>0 
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In view of the decomposition (6.39) for Vato, the estimates (6.42) (6.43) for ai, and 
the estimates (6.44) (6.45) for 02, we immediately obtain the estimate (6.12) for Vato and 
VlfO. This concludes the proof of Theorem 6.4. 



6.4. RegulEirity of the foliation with respect to u 

Let u{x,uj) the function constructed in section 6.2. In this section, we discuss the 
proof of Theorem 6.5 which deals with the control of the derivatives with respect to uj 
of the foliation of E provided by u{x^uj). Recall that {Ti,g,k) coincides with (]R^,5, 0) in 
\x\ > 2. Also, u{x,u!) coincides with x.(jj in \x\ > 2, and so a = 1, = a; and 9 = 
in this region. Thus, u clearly satisfies the estimates of Theorem 6.5 in > 2 and it is 
enough to control the derivatives with respect to uj of the function u{x, uj) solution to: 

tr^ — kpfN — 1 — a, on — 2 < u < 2, , . 

u{.,uj) — —2 on x.cu — —2, ^ ' ^ 

in the strip S = {x/ — 2 < u < 2}. 

6.4.1. First order derivatives with respect to u. The goal of this section is to 
prove (6.13). For the sake of simplicity, we only outline of the proof of the estimate for 
d^a. Differentiating the second equation of (6.25) with respect to ou, we obtain: 

(V^ - a-%)d^a = 2fV Na + 2i?jva^^ + • • • (6.47) 

where the first term on the right-hand side comes from the commutator [d^j, ^ (see [42]). 
Since yVAra and it! are in L^(E) respectively by (6.11) and (6.23), we obtain using in 
particular an parabolic estimate for the operator (Vat — o"^^) 

\\^ N^u,a\\L^T.) + ||y<9a;a||Ls°L2(p„) + Wf^ d^a\\L2{T.) < £. (6.48) 

Next, we differentiate (6.47) with respect to Vjv. We obtain: 

(Vat - a-^f ) Vat^^o = 2fWla + 2WNRNd^N + ■■■ ■ (6.49) 

The term ^NRNa^N may be treated using the contracted Bianchi identity for - as we 
did for VnRnn in section 6.2.1 - and turns out to be in L\H~^{Pu). On the other hand, 
in view of the estimate (6.12) for V^a, yV^a belongs to L'^H~2(^P^). We obtain using 
in particular a refined parabolic estimate for the operator (Vat — 

Finally, by interpolation between (6.48) and (6.50), we obtain d^^a in L'^Hi^Py) which 
embeds in L°°(E) since P„ has dimension 2. Together with (6.48) and (6.50), we obtain 
the estimate corresponding to ^^^a in (6.13). 
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6.4.2. Second order derivatives with respect to uj. The goal of this section is 
to prove (6.14). For the sake of simphcity, we only outhne of the proof of the estimate for 
d^a. Differentiating the equation (6.47) for 9^,0 with respect to uj, we obtain: 

(Vat - a"'A)'9'a = 2V%a + fVndu^a + 2Rd^Nd^N + ■■■ (6.51) 

where the first two terms on the right-hand side come respectively from the commutators 
[d^,y] and [5^,^]. Since R is in /.^(S) by (6.23), V%a is in LlH-^Pu) by (6.12), and 
^Nd^a is in LlH^P^) by (6.13), the right-hand side of (6.51) belongs to LlH--2{P^). 
Using in particular estimates for the parabolic operator (Vat — a~^^), we deduce 

+ + l|V^.a>||,,^_.^,^^ < (6.52) 

which is the estimate corresponding to d^a in (6.14). 

Remark 6.12. Note that we may not differentiate the equation (6.51) for d'^a with 
respect to Vat. Indeed, the term Viv-Rs^iva^Ar has no structure: unlike Rnn (^nd RNd^N 
which were involved in the equation for a and d^jO, Ra^Nd^N does not contain any contrac- 
tion with N since d^iN is tangent to P„. Thus, unlike V nRnn and V NRNd^N, we can 
not write V NRd^Nd^^N as a tangential derivative using the contracted Bianchi identities 
for R. Consequently, we can not obtain any estimate for V%d^a. 

6.4.3. Third order derivatives with respect to u. The goal of this section is to 
prove (6.15). Recall that div(iV) = tr6l, N = Vm/|Vm|, a = l/\Vu\ and tr6l = l-a + A;ArAf, 
so that: 

div (^) = 1 - + k^N. (6.53) 
\\\/u\J \Vu\ 

Differentiating (6.53) three times with respect to ou yields: 

(Vat - a-'4^)dlu = fd^ log(a) + • • • . (6.54) 

Using in particular (6.54), the estimate (6.14) for d'^a and refined parabolic estimates for 
the operator (Vat — obtain 

Now, since d^u e L^H'^(Pu) and P„ is 2-dimensional, we obtain that d^u belongs to 
L^^^(E), which is the desired estimate (6.15). 



6.5. A global coordinate system on P„ and E 

The goal of this section is to discuss the proof of Proposition 6.6. We start by con- 
structing a global coordinate system on P„. 
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6.5.1. A global coordinate system on P„. We have the following proposition 
Proposition 6.13. Let a; e S^. Let : P„ -> T^S^ defined by: 

^u{x) --^ diju{x,iju), (6.55) 

where T^^S^ is the tangent space to at ou. Then is a global diffeomorphism from 
Pu to T^S^. 

For the sake of simplicity, we only briefly sketch the proof. We start by showing that 
is a local diffeomorphism. We have 

fJac<|.FJac$ -a-2^^a^'^v.^) 9{d^N,d,N) 
^Jaccp^j jaccp, - a y ^^^^^^ ^^^^ ^^^^^^ ^^^^ 

where (v^,^') denotes the usual spherical coordinates on S^. Using the estimates (6.11) 
and (6.13), we are able to derive the following estimate 

||(Jac$„)^Jac$„ - /||loc(s) < £, (6.56) 

so that is a local diffeomorphism. In turn, this yields: 

|||det(Jac$„)| - l|U-(s) < e. (6.57) 

It remains to show that is onto and one-to-one. The proof relies on the estimate 
(6.56) for the Jacobian of the fact that u coincides with x.ui in the region \x\ > 2 and 
geometric considerations on the level sets of d^pU and d^u. We refer to [42] for the details. 

6.5.2. Proof of Proposition 6.6. Let uj G S^. Recall the definition (6.16) of 0^; : 

where has been defined in Proposition 6.13. The fact that (f)^ is a bijection is an easy 
consequence of the fact that is a bijection for all u. Then, it remains to prove (6.17). 

Wc arc able to obtain 

(Jac0t^)^Jac0tj = 

1 -d^ log(a) -d^ log(a) 

-d^ log(a) {d^ log(a))^ + g{d^N, d^N) log(a)5^ log(a) + g{d^N, d^N) 

-d^ log(a) \og{a)d^ log(a) + g{d^N, d^N) {d^ log(a))^ + g{d^N, d^N) 

Taking the determinant yields: 

det((Jac0(^)'^Jac0(^) = det((Jac$u)"^Jac$„), (6.58) 

which together with (6.56) and the estimate (6.11) for a implies: 

||det((Jac0^)^Jac0^) - 1||loo(s) < s. (6.59) 

(6.59) yields (6.17). This concludes the proof of Proposition 6.6. 



X 



124 



6. CONTROL OF THE FOLIATION AT INITIAL TIME 



6.6. Additional estimates 

The proof of Proposition 6.7 and Proposition 6.8 follow from the estimates of Theorem 
6.3, Theorem 6.4, and Theorem 6.5 using also for some estimates the fact that u coincides 
with x.uu in the region > 2 or the properties of the Littlewood-Paley projections Pj. 
For the sake of simphcity, we skip these proofs and refer the reader to [42] for the details. 



CHAPTER 7 



The Strichartz estimates 



Recall Steps A, B, C and D introduced in section 1.2.4. In this chapter, we perform 
Step D, i.e. we prove Proposition 2.32. More precisely, let j > 0, and let ip a smooth 
function on supported in 

^ < lei < 2. 

Let (pj the parametrix (1.24) with an additional frequency locahzation A ~ 2^ 

/poo 
/ e'^"(*'^''^V(2"^A)/(Aa;)A2rfArfa;, (7.1) 
_2 Jo 

where u{.,.,uj) is a solution to the eikonal equation g°''^daudfsu = which depends on 
an extra parameter cu G E>^. Assume that the space-time A4 is foliated by space-like 
hypersurfaces defined as level hypersurfaces of a time function t. Let {p,q,r) such that 
p,q > 2, q < +00, and 

111 313 

- + -<-, r = . 

p q 2 2 p q 

In this chapter, we outline the proof of the following sharp^ Strichartz estimates 

lk.llL|;_,,L.(Eo < 2nm^-^x)f\\mm- (7-2) 

The proof of Proposition 2.32 will then be a simple consequence of (7.2) with the choice 
p = g = 4. 

Remark 7.1. Even though we only need L^(A4) Strichartz estimates - which corre- 
sponds to p — q — A in (7.2) - to prove Proposition 2.32, it turns out that this particular 
case is not easier to prove than the general case. 

7.1. Assumptions on the phase M(t,a;,a;) and main results 

7.1.1. Time foliation on }A. We foliate the space-time M. by space- like hypersur- 
faces Tit defined as level hypersurfaces of a time function t. We assume < t < 1 so 
that 

[jTt. (7.3) 

0<t<l 

^Note in particular that the corresponding estimates in the flat case are sharp. 
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We denote by T the unit, future oriented, normal to Sj. We also define the lapse n as 

= T{t). (7.4) 

Note that we have the following identity between the volume element of M. and the 
volume element corresponding to the induced metric on 

dM = nd^tdt. (7.5) 

We will assume the following assumption on n 

l<n<2 (7.6) 

which together with (7.5) yields 

dM ~ dEt dt. (7.7) 

Remark 7.2. The assumption (7.6) is very mild. In particular, it is compatible with 
the estimates for n derived in [44] (see also (4.40) 

7.1.2. Geometry of the foliation generated by u on M.. Recall that w is a 

solution to the eikonal equation g'^^daudgu = on depending on an extra parameter 
u G S^. The level hypcrsufaces u{t, x,uj) = u of the optical function u are denoted by 1-Lu- 
Let L' denote the space-time gradient of m, i.e.: 

L' = g'^^dpud^. (7.8) 

Using the fact that u satisfies the eikonal equation, we obtain: 

T>L'L' = 0, (7.9) 

which implies that L' is the geodesic null generator of 1-Lu- 
We have: 

T{u) = ±|Vm| 

where \Vu\^ = Yl^^i \ei{u)\'^ relative to an orthonormal frame on E^. Since the sign of 
T{u) is irrelevant, we choose by convention: 

T(u) = -\Vu\ (7.10) 

so that u corresponds to —t + x • ou in the flat case. 
Let 

L = bL' = T + N, (7.11) 

where L' is the space-time gradient of u (7.8), b is the lapse of the null foliation (or shortly 
null lapse) 

b-^ = - < L', T >= -T{u), (7.12) 
and is a unit vectorfleld given by 

N-^- (7.13) 
Note that we have the following identities. 



7.1. ASSUMPTIONS ON THE PHASE u{t,x,Lo) AND MAIN RESULTS 
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Lemma 7.3. 

L(u) = 0, L{d^u) = (7.14) 

and 

giN,d^N) = 0. (7.15) 
The proof is elementary and can be found in [46] . 

7.1.3. Regularity assumptions for u{t,x,u}). We now state our assumptions for 
the phase u{t,x,u). These assumptions are compatible with the regularity obtained for 
the function u{t,x,uj) constructed in [44] (see also (4.42), (4.48), (4.49)). Let < £ < 1 
a small enough universal constant, b and N satisfy 

\\b-l\\L^ + \\d^b\\L^<e. (7.16) 

\\g{d^N,d^N)-l2\\L'>o<e. (7.17) 

\N{.,uj)-N{.,uj')\ = \uj-u'\{l + 0{s)). (7.18) 

Remark 7.4. In the flat case, we have Ai = (M"^^^,m), where m is the Minkowski 
metric, u(t, x,u) = —t + x ■ u, b = 1, N = u and L = dt + cu ■ dx- Thus, the assumptions 
(7.16) (7.17) (7.18) are clearly satisfied with e = 0. 

Remark 7.5. In terms of the regularity of u{t,x,u}), the assumptions (7.16) (7.17) 
correspond to 

Vu e and Vd^u e L°° 

which is very weak. In particular, the classical proof for obtaining Strichartz estimates for 
the wave equation relies on the stationary phase for an oscillatory integral involving u as 
a phase, and typically requires at the least one more derivative for u (see Remark 7.11). 

7.1.4. A global coordinate system on E^. For all < i < 1, and for all a; e S^, 

{u{t,x,uj),d^u{t,x,uj)) is a global coordinate system on Sj. Furthermore, the volume 
element is under control in the sense that in this coordinate system, we have 

^ < \/detc/ < 2 (7.19) 

where g is the induced metric on Ej, and where det^f denotes the determinant of the 
matrix of the coefficients of g. 

Remark 7.6. In the flat case, we have = {t} x IR'^ and u(t,x,uj) = —t + x ■ u so 
that {u{t,x,uj),di^u{t,x,uj)) is clearly a global coordinate system on Ej and deig = 1 in 
this case. These assumptions are also satisfied by the function u{t, x, ou) constructed in 
[44]. 
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7.1.5. Main results. We next state the main result of tliis cliapter concerning gen- 
eral Strichartz inequalities in mixed space-time norms of the form L^^ -^jL^(Et) defined as 
follows, 

ii^iUf„,,L.(Eo = \\nt,-)\\u^,)dt 

Theorem 7.7. Let (p, q) such that p,q>2, q < +oo, and 

111 

- + -<-. 
P 

Let r defined by 

3 1 3 



2 p q 

Then, the parametrix localized at frequency j defined in (7.1) satisfies under the as- 
sumptions (7.6); (7.16), (7.17); (7.18) and the assumptions in section 7.1.4 the following 
Strichartz inequalities 

lk.l|Lf„ ^L.(E.) < 2^>||V'(2-^A)/|U.(M3). (7.20) 

Proposition 2.32 - which corresponds to Corollary 2.8 in [46] - is then a simple conse- 
quence of Theorem 7.7, see [46] for the details. 

The rest of the chapter is organized as follows. In section 7.2, we use the standard 
TT* argument to reduce the proof of Theorem 7.7 to an upper bound on the kernel K of 
a certain operator. This kernel is an oscillatory integral with a phase 0. In section 7.3, 
we prove the upper bound on the kernel K provided wc have a suitable lower bound on 
(p. Finally, in section 7.4, we prove the lower bound for (p used in section 7.3. 

7.2. Proof of the Stricheirtz estimates 

The goal of this section is to prove Theorem 7.7. We start with the following remark. 

Remark 7.8. Fixing a global system of coordinates x — {x^^x'^.^x^) in T^, such as 
the one described in section 7.1.4, we note in view of (7.19) that (7.20) is equivalent 
with the same inequality where the norm L'^{'Et) on the left-hand side is replaced by the 
corresponding euclidean norm in the given coordinates. More precisely we can assume 
from now on that 




which we will denote by a slight abuse of notation by 

ll^l|Lfo^,L.(M3)- 

Note also that in the {t, x) coordinates M. — [0, 1] x R' 




7.2. PROOF OF THE STRICHARTZ ESTIMATES 129 
For convenience, let us introduce the operator Tj acting on functions / G L^(M^), 

Tjf{t,x)^ / e^^"(*'^''^V(2~^A)/(Aa;)A2(iA(ia;. (7.21) 
Js2 Jo 

Note in particular that 

Tjf = (7.22) 

where Lpj is the paramctrix localized at frequency 2^ defined in (7.1). To prove Theorem 
7.7, we rely on the standard TT* argument for the Fourier integral operator (7.21). Note 
that the operator Tj takes real valued functions h on M. to complex valued functions on 

Jm 

Therefore, the operator Uj := TjT* is given by the formula, 

Ujh{t,x)= / / / e^^"(*'^''")-^^"("*'"V(2~'A)'/i(s,y)A'dAda;dsdy. 
Note, in view of Remark 7.8 and (7.22), that (7.20) is equivalent to the following estimate 

II^.J^I|Lf,,L.(M3)<2^^n|/^IL.;^^..'(M3)> (7-23) 
where p' (resp. q') is the conjugate exponent to p (resp. q). Observe that, 

^^K^' ^) ^ L [ L ^-'<---'^—<---'^^^>^ (^. I) 

with 2^M. — [0,2.^] X relative to the rescaled variables {s,y). Thus, setting, 

J§2 Jo J2iM 

we have 

We easily infer that (7.23) is equivalent to the estimate. 
We introduce the kernel K oi A 

ri(X) 



K{t,x,s,y)= / e'^''''y^'^'''r'^''''y^'^''')^|J{XyX'dXduJ. (7.25) 
Js2 Jo 

Remark 7.9. In the flat case, we have u{t, x,uj) — —t + x ■ uj so that 
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In particular, K is independent of j 

is^ Jo 

We have the following proposition. 

Proposition 7.10. The kernel K of the operator A satisfies the dispersive estimates, 
\K{t,x,s,y)\ < V(i,x) e 2^M, y{s,y) e 2^M. (7.26) 

The proof of Proposition 7.10 is postponed to section 7.3. We now conclude the proof 
of Theorem 7.7. (7.24) follows from (7.26) using interpolation and the Hardy- Littlewood 
inequality according to the standard procedure, see for example [38] and [39]. Finally, in 
view of the discussion above, (7.24) yields (7.23) which in turn implies (7.20) in view of 
(7.22). This concludes the proof of Theorem 7.7. 

7.3. Upper bound on the kernel K 

The goal of this section is to prove Proposition 7.10. Let the scalar function on 
M X M xEP' defined as 

(f){t, X, s, y, uj) — u{t, X, uj) — u{s, y, ou). (7.27) 
In view of (7.25), we may rewrite K as 

K{t,x,s,y)= [ r e'^''<^'^'^'^''')x''dXdw. 

Jo 

After integrating by parts twice in A, and using the size of the support of ip, this yields 

mt,x,s,y)\< [ — —TTtV^^^^^- (7-28) 
1 + 22^0 (I, It, a;) 

The next section is dedicated to the obtention of a lower bound on |0| which will allow 
us to deduce (7.26) from (7.28). 

Remark 7.11. It is at this stage that we depart from the standard strategy for proving 
Strichartz estimates. Indeed, the usual method consists in using the stationary phase 
method to derive (7.26). To this end, one considers the neighborhood in of stationary 
points cuq, i.e. such that d^(l)\^^^^ = 0. One then needs an identity of the type 

(j)= {s- t)A{Lu - uo) ■{u-uo) + o ((s -t){uj~ uof) (7.29) 

foru in the neighborhood of ujq and for some 3x3 invertible matrix A. (7.29) then allows 
to perform a change of variables in uj which ultimately leads to (7.26). In particular, the 
standard method requires at the least^ dt^xd^u £ L°° just to derive (7.29). 



^One also needs to take care of the contribution to K of the angles w S corresponding to the 
exterior of the neighborhood of stationary points which may increase the needed regularity. 
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Our assumptions correspond only to dt^x^uiU G L°°. Thus, in order to obtain (7.26), we 
instead integrate by parts in X to obtain (7.28), and then look for a suitable lower bound 
on 101 . In particular, we obtain lower bounds of the following type (see details in Lemma 
7.19) 

101 > \s-t\\uj-ujo\'' (7.30) 

for uj in the neighborhood of some ujq e S^. The fundamental observation is that, as it 
turns out, the inequality (7.30) requires less regularity than the equality (7.29). 

7.3.1. The key lemma. Let {t,x) and {s,y) in and let lo e S^. In this section, 
we obtain a lower bound on 0(t, x, s, y, uj). We may assume 

< t < s < 1. 

Definition 7.12. For any w G and c g M, let 7a,(cr) denote the null geodesic 
parametrized by the time function and with initial data 

7^(0) = it.x), = b-\t,x,uj)L{t,x,uj). 

Definition 7.13. Let us define the subset S of Eg as 

We also define for all (s, z) G 

m{s, z) — ma.x{u{s, z, uj) — u{t, x,uj)). (7.32) 

We have the following lemma characterizing the zeros of m (see [46] for a proof). 
Lemma 7.14. We have 

S^ipeJ^s, /m(p) = 0}. 
Next, we define the following two subsets of Eg 

Aint = {pe ^s/m{p) < 0}, Aext = {pe /m(p) > 0}. (7.33) 
Note in view of Lemma 7.14 that 

E, = 5 U A,nt U Aexf (7.34) 

Remark 7.15. In the flat case, the picture is the following: 

(1) The null geodesies^ 7^ span the light cone from {t,x). In particular, the null 
geodesies do not intersect except at {t,x). 

(2) S is the intersection^ of the forward light cone from {t,x) with {s} x R^. 

(3) Aint and A^xt correspond respectively to the interior and the eocterior of S. 



■^which are straight lines in this case 
^5 is a sphere in this case 
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Note that we do not need to prove these statements in our case. This is fortunate since 
these statements - while probably true in our general setting - would be delicate to establish 
(see for instance [25] for a proof of (1) on a space-time g) with limited regularity) . 

Next, we introduce some further notations. First, we denote by tuq the value of m at 
(s,y), i.e. 

tuq — m.&x.{u{s,y,uj) — u{t,x,uj)). (7.35) 

We also denote by loq an angle in §^ where the maximum in (7.35) is achieved, i.e. 

mo = u{s,y,UQ) - u{t,x,ujQ). (7.36) 

Remark 7.16. In the flat case, ooq is unique and corresponds to the angle of the 
projection of (s, y) on S. Again, while this may he also true in our general setting, we do 
not need to prove this statement in our case. 

Note that if {s,y) G A^xti the function u{s,y,u!) — u{t,x,u!) may change sign as cu 
varies on S^. We define 

L> = G SV ^> ^) = y, w)}- (7.37) 
The following lemma gives a precise description of D (see [46] for a proof). 

Lemma 7.17. Let (s.y) G A^^t. Let D defined as in (7.37). Let {0,Lp) denote the 
spherical coordinates with axis cuq. Then, there exists a 27t -periodic function 

9i : [0,27r) ^ (0,7r) 

such that in the coordinate system {0,(fi), D is parametrized by 

D^{e^ei(ip), 0<(^<27r}. 

Remark 7.18. In the flat case, recall that u{t,x,u) ~ —t -\- x ■ uj. In this case, one 
easily checks that D is a circle of axis ujq on the sphere S which is generated by the tangents 
to S through y (see figure 1). 

Let u G S^. According to Lemma 7.17, the great half circle on originating at 
and containing ui intersects D at a fixed point oji. Let 9 and 9i respectively denote the 
positive angles between ujq and uj (resp. ujq and uji). 

In order to obtain a lower bound for |0|, we will argue differently according to whether 
(s, y) belongs to the region S, Aint or Aext- 

Lemma 7.19 (Key lemma). |0| satisfies the following lower bounds 
(1) If {s, y) G S, we have 

|0(i, x, s, y,a;)| > -|t — s| |a; — o'ol^. (7.38) 
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Figure 1. Representation of D in the flat case 

(2) // {s,y) e Aint, we have 

\(t){t,x,s,y,uj)\ > l\t - s\\u - uol'^ . (7.39) 

o 

(3) // {s,y) G Aext and 9i < 9 < tt , we have 

\(t>(t,x, s,y,ijj)\ > - \t — s\\uj — Uil"^ . (7.40) 

(4) // (s, y) G Aext and < 9 < 9i, we have 

1 — cos(9 — 9i) 

m,^,s,y,u) > J— ^-jTTY^rno 7.41 

y 1 - cos(6'i) 

Tlie proof of Lemma 7.19 is postponed to section 7.4. 

Remark 7.20. The proof of Lemma 7.19 is inspired by the overlap estimates for wave 
packets derived in [35] and [36] in the context of Strichartz estimates respectively for C^'^ 
and H^'^^ metrics. Note however that the estimates in these papers rely heavily on a direct 
comparison of various quantities with the corresponding ones in the flat case. Such direct 
comparisons do not hold in our framework. Here, the closeness to the flat case manifests 
itself in the small constant e in the right-hand side of (7.16), (7.17) and (7.18), and in 
the existence of the global coordinates systems of section 7.1.4. 

7.3.2. Proof of Proposition 7.10. Recall that we need to show that the kernel K 
defined in (7.25) satisfies the upper bound (7.26). To this end, we will use the estimate 
(7.28) for K together with the estimates provided by Lemma 7.19. We argue differently 
according according to whether {s,y) belongs to S, A^t or Af,xt- 

If {s,y) belongs to S, we have the lower bound (7.38) for |0| 

\4>{t,x,s,y,u)\ > - - Wop, 
where Uq G is an angle satisfying (7.36). Then, we deduce 
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Together with (7.28), this yields 

\K{t,x,s,y)\< [ ^. 

y§2 l + \t- s\^\u-uo\^ 

Using the spherical coordinates {9, tp) with axis cuq, we obtain 

M / r sin(^)rf^ 

Performing the change of variables 

|i-s|(l-cos(^)) 

we obtain 
This implies 

\K{t,x,s,y)\<j^^, yit,x)e2^M, V |) G 5 (7.42) 

which is the desired estimate. 

The estimates corresponding to the cases where {s,y) belongs to Aint or 
similar (see [46] for the details). This concludes the proof of Proposition 7.10. 

7.4. Lower bound for \(f)\ 

The goal of this section is to prove Lemma 7.19. The main ingredients of the proof are 
already present in the flat case. Thus, to simplify the exposition, we will prove Lemma 
7.19 for the phase function u = —t + x ■ co oi the flat case. We will then explain what are 
the modifications in the general case (see Remark 7.22). We refer to [46] for the proof in 
the general case. 

7.4.1. A lower bound for |0| when {s,y) E S (proof of (7.38)). In view of the 
definition of mo in (7.35) and ouo in (7.36), we have in the fiat case 

y — ^ 

u(t, x,uj) — —t + X • uj, ujo — -. r, mo — —{s — t) + \y — x\. (7.43) 

\y - x\ 

If (s, y) e S, we have \y — x\ — s — t. Together with (7.43), this yields 

u{s,y,u!) — u{t,x,uj) = —s + t+{y — x)-u 

— —{s — t) + {s — t)uo ■ cu 

= -^{s-t)\uj-uJo\'' (7.44) 

which is the desired estimate (7.38). 
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7.4.2. A lower bound for |0| when {s,y) E Am (proof of (7.39)). (7.43) yields 

u{s,y,u!) — u{t,x,u) = —{s — t) + {y — x)-u 

= — (s — t) + \y — x\u ■ Uq 

= -{s-t) + \y-x\-^\x-y\\uj-uJo\^. (7.45) 
Now, ii \x — y\ < \{s — t) we have, 

3 11 1 

u{s,y,u) -u{t,x,uj) < --{s-t) + -{s-t)-\uj-uJof <--{s-t). 

On the other hand, ii \x — y\ > \{s — t) 

~ 2 



u{s,y,u!) — u{t,x,u!) < — \x — y\\u! — u!q\'^ < —-{s — t)\u! — COqI'^. 



Thus, in both cases, 

u{s,y,u) - u(t,x,u) < -^\x -y\\LJ -cuq]^ < -^{s - t)\uj - ujo\'^ 

which is the desired estimate (7.39). 

7.4.3. A lower bound for when {s,y) e A^xt (proof of (7.40) (7.41)). (7.43) 
yields 

u{s,y,uj) — u{t,x,uj) — —{s — t) + {y — x)-uj 

= —{s — t) + \x — y\uj-uJo 

= — (s — t)\oj — ojq\'^ + mou ■ Uq. (7-46) 

Recall that we have defined the set D by 

D — {u! eEi^ / u{t, X, uj) — u{s, y, u)}. 

Also, for fixed uq, uj wc defined ui E D io lie on the same plane great circle of as uq, oj. 
Clearly, since oui e D, and in view of (7.43), (7.46) and the definition of D, we have 

001-000 = (7.47) 

\x-y\ 

Fix now ooi & D and let z — 7a,i(s — t), i.e. 

z^x + {s- t)ooi e S. (7.48) 
Note that in view of the definition of D, 

{y- z)-ui^ -{s -t) + {y-x)-uji^ u{s, y, uji) - u{t, x, Ui) = 0. 
Hence, with the notation 

Vo = y- z, 
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we obtain 

• cui = 0. (7.49) 

Now, we have 

u{s,y,uj) — u{t,x,uj) — u{s,y,uj) — u{s, z,uj) + u{s, z,uj) — u{t,x,uj) 

— vo ■ ou + u{s, z,uj) — u(t,x,uj) 

— Vo ■ {uj — uJi) + u{s, z,uj) — u{t,x,uj). (7.50) 

Note that since 2; G we can apply the estimate obtained in section 7.4.1. Since the 
maximum in m(s, z) is attained at a; = a;i, we have 

u{s, Z, Oj) — u{t, X,Oj) — —-{s — t)\0J — OJi\'^ (7-51) 

and we infer that, 

ii(s,y,a;) — w(i,x, a;) = — -(s — i)|a; — + vo • (u; — a;i). (7-52) 

Recall that we have also denoted by Q and Q\ the positive angles between wq, u; and 
respectively cuq, Wi. If ^1 < ^ < tt - which corresponds to v^- {uj — < - we have in 
view of (7.52) 

u^s^y^iJ) — uit^x^Lo) < —-{s — t)\u) — u)i\'^ (7.53) 

which is the desired estimate (7.40). 

The delicate case is when < ^ < ^1 which corresponds to 

VQ-{yj- wi) > 0. (7.54) 

In the rest of the proof, we assume (7.54), and we focus on the remaining estimate (7.41). 
In view of the definition of cjq and (7.52), we have 



-{s-t)^\x-y\ = u{s,y,ujQ) -u{t,x,ujo) 
_ 1 
~ ~2 



^{S- t)\oJo - OJif + Vq ■ {OJQ - OJi). 



Thus, 



mo^ -{s-t) + \x-y\ = -^{s -t)\ujQ - uji\^ + vq ■ {ujQ - uji). (7.55) 
Since mo > we deduce, 

Vq ■ {ujo — COi) > 0. (7.56) 

Let a, ai be the positive angles between vo and (uj — uji) and, respectively vq and {ujq — uji). 
In view of (7.54) and (7.56) we infer that 

< q;,q;i < 7r/2. (7.57) 



Figure 2. Definition of tlie angles 6i and ai 




Figure 3. Definition of tlie angles 6 and a 
Also, in view of (7.47) we have, 

< ^1 < 7r/2. 
Simple considerations on angles imply^ (see figure 2), 

01 = 2ai. (7.5^ 

Therefore, 

mo = -^{s-t)\ujo-uJi\'^ + \z-y\\ujo-uJi\cos(^^ 

and 



I^ol = ^ . 7g-. • (7.59) 

1^0 — UJi \ COS (^j 

Using the same type of argument^ as in (7.58) we also deduce (see figure 3) 

a = (7.60) 
Therefore, according to (7.52), (7.59) and (7.60), we obtain 



^Let ipi the angle defined on figure 2. Then 2ipi + 9i = tt, and (fi + ai = ^. Hence 9i — 1a\ 
^Let the angle defined on figure 3. Then lip + 16*1 — 6*1 = tt, and + a = |. Hence \d\ ~ 9\ — 2a 
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/ 

u{s,y,uj) — u{t,x,uj) — — -{s — t)\uj — uji\'^ + \vo\\uj — uji\ cos 

- + (7.61) 



|u;o-'^i| cos(|-) 2 



where A{lo) is given by 



cos { ^1^) 

A(a;) = -|a;-a;i| + ^^^^ - (7.62) 

cos (^) 

We have the following lemma (see [46] for a proof). 
Lemma 7.21. For allO <9 < 9i, we have 

A(uj) > 0. (7.63) 

Back to (7.61), we thus derive, 

, , , , Iw - Wi| cos (%^) 

u{s,y,uj) -u{t,x,uj) > r T^^^o- 

1^0 -'^i I cos (^j 

Using our angle restriction 

< ^ < ^1 < |, 

we deduce 

u{s, y, uj) - u{t, x,uj)>^ — — ^mo. (7.64) 

Z \0Oq — OO'i I 

Since 6 is the angle between u and Uq, and 6i is the angle between ui and Uq, wc have 

I = \/2v'l -cos(^i -^), - a;o| = v^Vl - cos(^i). (7.65) 
In view of (7.65), we can rewrite (7.64) in the form. 



At-f ^> /l-cos(^-ei) 

(f){t,x,s,y,uj)>moA— y—— (7.66) 

y 1 — cosi^^ij 

which is the desired estimate (7.41). This concludes the proof of Lemma 7.19 in the flat 
case. 

Remark 7.22. Let us indicate how to prove Lemma 7.19 in the general case. The 
whole point is to realize that the only estimates for which the precise regularity ofu matters 
are the ones corresponding to (7.44), (7.45), (7.46), (7.50) and (7.51). Indeed, once this 
has been achieved, the rest of the argument is then essentially the one of the flat case. 

Now, to prove the estimates corresponding to (7.44), (7.45), (7.46), (7.50) and (7.51) 
in the general case, one needs the following two additional ingredients (see [46] for the 
details): 



7.4. LOWER BOUND FOR \<f)\ 
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(1) These estimates are obtained by using the following standard identity 

u{r){l),uj)^u{r){0),uj)+ [ g{Lr,^,),r)'{a))da, (7.67) 

^0 

where rj is a curve in M., and where L denotes the space-time gradient of u. 
It turns out that one may chose suitable curves'^ rj allowing us to deduce from 
(7.67) the estimates corresponding to (7.44), (7.45), (7.46), (7.50) and (7.51) 
under our assumptions (7.16), (7.17) and (7.18). This changes the constants in 
the inequalities due to the presence of additional 0(e) terms, but does not change 
the nature of the estimates for £ > small enough. 

(2) The above mentioned curves rj start on S , and a crucial point is to check that 
such curves end up exactly at {s,y). To this end, one uses the global coordinate 
system {u{t,x,(x!Q),di^u{t,x,(x!o)) of section 7.1.4 on for a well-chosen angle 
ujq e S^, which allows us to identify (s, y) as the unique point p on such that 

u{p,Uo) = u{s,y,uJo) and d^u{p,Uo) = d^u{s,y,u}o). 



In the flat case, the corresponding curves r] are straight hnes. 
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